Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,451,602 C→G 18.7% E1279Q (GAA→CAA)  rpoB ← RNA polymerase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,451,6020C→G18.7% 180.5 ‑0.4 ‑0.9283E1279Q (GAA→CAA) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  new base (27/26):  ref base (108/122):  total (135/148)
Fisher's exact test strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.19e-01

Only 100 of 291 total aligned reads displayed.

TTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACG  >  NC_007779/3451502‑3451701
                                                                                                    |                                                                                                   
ttGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTc                                                                                                     >  1:6228208/1‑101 (MQ=255)
ttGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTc                                                                                                     >  1:6122973/1‑101 (MQ=255)
  gCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCa                                                                                                   <  1:9694574/101‑1 (MQ=255)
  gCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCa                                                                                                   >  1:8259442/1‑101 (MQ=255)
  gCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCa                                                                                                   <  1:3579643/101‑1 (MQ=255)
  gCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCa                                                                                                   <  1:1780456/101‑1 (MQ=255)
  gCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCa                                                                                                   <  1:12400068/101‑1 (MQ=255)
    cGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAgc                                                                                                 >  1:4480438/1‑101 (MQ=255)
       ccACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGcc                                                                                              >  1:5529979/1‑101 (MQ=255)
        cACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGccc                                                                                             >  1:6956173/1‑101 (MQ=255)
        cACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGccc                                                                                             <  1:12998717/101‑1 (MQ=255)
          cGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCcac                                                                                           <  1:8940758/101‑1 (MQ=255)
          cGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCcac                                                                                           >  1:12811254/1‑101 (MQ=255)
             tGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACt                                                                                        <  1:9421013/101‑1 (MQ=255)
                   tataCATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCAtc                                                                                  <  1:10678317/101‑1 (MQ=255)
                   tataCATCTTAGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCAtc                                                                                  >  1:9197732/1‑101 (MQ=255)
                     taCATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCAtctc                                                                                <  1:663187/101‑1 (MQ=255)
                      aCATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTcc                                                                               <  1:10463338/101‑1 (MQ=255)
                        aTCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTcccc                                                                             <  1:6255800/101‑1 (MQ=255)
                        aTCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTcccc                                                                             >  1:6834294/1‑101 (MQ=255)
                        aTCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTcccc                                                                             >  1:14186603/1‑101 (MQ=255)
                        aTCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTcccc                                                                             <  1:13049856/101‑1 (MQ=255)
                         tCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCg                                                                            <  1:5929016/101‑1 (MQ=255)
                           ttGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGaa                                                                          <  1:9377906/101‑1 (MQ=255)
                            tGGTACGACCGTTCACGTCATCAGACTTAGCGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGaaa                                                                         >  1:7175444/1‑101 (MQ=255)
                             ggTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAAc                                                                        <  1:11482843/101‑1 (MQ=255)
                                 cGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTg                                                                    <  1:860525/101‑1 (MQ=255)
                                 cGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTg                                                                    <  1:581175/101‑1 (MQ=255)
                                  gACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGa                                                                   >  1:11331576/1‑101 (MQ=255)
                                    ccGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGacc                                                                 <  1:4737236/101‑1 (MQ=255)
                                     cGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGacca                                                                <  1:8022789/101‑1 (MQ=255)
                                     cGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGacca                                                                <  1:12667464/101‑1 (MQ=255)
                                      gTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGaccac                                                               >  1:10295287/1‑101 (MQ=255)
                                       ttCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGaccacc                                                              >  1:3409716/1‑101 (MQ=255)
                                       ttCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGaccacc                                                              >  1:2114296/1‑101 (MQ=255)
                                       ttCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGaccacc                                                              <  1:2802301/101‑1 (MQ=255)
                                       ttCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGaccacc                                                              >  1:5218806/1‑101 (MQ=255)
                                         cACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGa                                                            >  1:1106538/1‑101 (MQ=255)
                                           cGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAAc                                                          >  1:4203476/1‑101 (MQ=255)
                                           cGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAAc                                                          >  1:9123319/1‑101 (MQ=255)
                                           cGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAAc                                                          <  1:4533439/101‑1 (MQ=255)
                                            gTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACt                                                         >  1:10561546/1‑101 (MQ=255)
                                             tcatcgGACTTGATGGTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACtg                                                        <  1:12918595/101‑1 (MQ=255)
                                             tcatcaGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACtg                                                        <  1:4866694/101‑1 (MQ=255)
                                               atcaGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACtgtg                                                      >  1:7774656/1‑101 (MQ=255)
                                                tcaGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGc                                                     <  1:10356686/101‑1 (MQ=255)
                                                tcaGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGc                                                     <  1:7512857/101‑1 (MQ=255)
                                                  ggACTTGATGGTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCt                                                   <  1:3310051/100‑1 (MQ=255)
                                                    aCTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTa                                                 >  1:1966938/1‑101 (MQ=255)
                                                       tAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTacca                                              <  1:9479115/101‑1 (MQ=255)
                                                        aaCGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTaccac                                             <  1:13369897/101‑1 (MQ=255)
                                                         atggTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTaccacg                                            <  1:2324968/99‑2 (MQ=255)
                                                         aCGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTaccacc                                            >  1:608533/1‑101 (MQ=255)
                                                           ggTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCa                                          <  1:5660813/101‑1 (MQ=255)
                                                           ggTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACCCCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCa                                          <  1:13470428/101‑1 (MQ=255)
                                                            gTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAg                                         >  1:2193434/1‑101 (MQ=255)
                                                             tGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGc                                        >  1:6576956/1‑101 (MQ=255)
                                                             tGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGc                                        >  1:7057087/1‑101 (MQ=255)
                                                              gAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCg                                       >  1:1143307/1‑101 (MQ=255)
                                                              gAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCg                                       <  1:10922588/101‑1 (MQ=255)
                                                              gAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCg                                       <  1:13807084/101‑1 (MQ=255)
                                                                 cagctCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGgct                                    <  1:2645028/97‑1 (MQ=255)
                                                                 cATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgct                                    <  1:7134941/101‑1 (MQ=255)
                                                                 cATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgct                                    <  1:14262544/101‑1 (MQ=255)
                                                                  agctCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGgctg                                   <  1:2528316/98‑1 (MQ=255)
                                                                  aTTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgctg                                   <  1:7347439/101‑1 (MQ=255)
                                                                  aTTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgctg                                   <  1:6511984/101‑1 (MQ=255)
                                                                    ttCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgctgct                                 <  1:7992006/101‑1 (MQ=255)
                                                                      ccTGCAGGGTGTAGGCAGCGCTGTATGCCTGCATTGCCCACACCTCCATCTCTCCGAAGCTCTGACCACCGAACTGAGCCTTACCACCCAGCGGctgctgg                               >  1:3109068/1‑100 (MQ=255)
                                                                       cTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAg                              <  1:5369643/101‑1 (MQ=255)
                                                                        tGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGt                             <  1:4012802/101‑1 (MQ=255)
                                                                        tGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGt                             <  1:2986715/101‑1 (MQ=255)
                                                                         gCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTa                            >  1:10051738/1‑101 (MQ=255)
                                                                         gCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTa                            <  1:1852493/101‑1 (MQ=255)
                                                                           aGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAAc                          >  1:10751202/1‑101 (MQ=255)
                                                                           aGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAAc                          >  1:12218251/1‑101 (MQ=255)
                                                                           aGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAAc                          >  1:2832184/1‑101 (MQ=255)
                                                                            gggTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAAcc                         >  1:5893773/1‑101 (MQ=255)
                                                                             ggTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCa                        <  1:11212573/101‑1 (MQ=255)
                                                                             ggTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCa                        <  1:4260256/101‑1 (MQ=255)
                                                                                gtAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCatgg                     >  1:1361482/1‑98 (MQ=255)
                                                                                 tATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCt                    >  1:12513851/1‑101 (MQ=255)
                                                                                  aggcagcgcCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCatggag                   >  1:10626238/6‑96 (MQ=255)
                                                                                   tgctgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGt                  <  1:8457017/101‑1 (MQ=255)
                                                                                    gctgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTa                 >  1:13892874/1‑101 (MQ=255)
                                                                                      tgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACTAGGCTGTAAg               >  1:493969/1‑101 (MQ=255)
                                                                                      tgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAg               >  1:4317031/1‑101 (MQ=255)
                                                                                       gcgcCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGa              >  1:11889021/1‑101 (MQ=255)
                                                                                       gcgcCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCatggagtacgg              >  1:10084120/1‑91 (MQ=255)
                                                                                        cgcCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGaa             >  1:3714459/1‑101 (MQ=255)
                                                                                            gTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCgg         <  1:9022672/101‑1 (MQ=255)
                                                                                            gTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCgg         >  1:11268666/1‑101 (MQ=255)
                                                                                             tATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGt        <  1:7892306/101‑1 (MQ=255)
                                                                                              aTGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTg       <  1:2595052/101‑1 (MQ=255)
                                                                                              aTGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTg       <  1:2697894/101‑1 (MQ=255)
                                                                                              aTGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTg       <  1:11384985/101‑1 (MQ=255)
                                                                                                 cTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGaa    <  1:4069545/101‑1 (MQ=255)
                                                                                                   tCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAAcg  <  1:2839492/101‑1 (MQ=255)
                                                                                                   tCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAAcg  <  1:640496/101‑1 (MQ=255)
                                                                                                   tCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGGGGTAcg  >  1:8962467/1‑101 (MQ=255)
                                                                                                    |                                                                                                   
TTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACG  >  NC_007779/3451502‑3451701

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.