Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,184,280 A→C 7.1% E1012D (GAA→GAC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,2800AC 7.1% 334.9 / 9.6 125E1012D (GAA→GAC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (61/51);  new base C (3/7);  total (64/61)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAG  >  NC_000913/4184131‑4184427
                                                                                                                                                     |                                                                                                                                                   
gaaCTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTTAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGACACAAATTCAGCTGTaa                                                                                                                                                     >  2:738801/1‑150 (MQ=255)
gaaCTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGaa                                                                                                                                                     >  1:937126/1‑150 (MQ=255)
 aaCTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAAc                                                                                                                                                    <  1:574881/150‑1 (MQ=255)
 aaCTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAAc                                                                                                                                                    <  2:63425/150‑1 (MQ=255)
  aCTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACa                                                                                                                                                   >  2:1410916/1‑150 (MQ=255)
     gCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAgct                                                                                                                                                >  2:1403706/1‑150 (MQ=255)
      cAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAgctg                                                                                                                                               <  2:810343/150‑1 (MQ=255)
        taTCCTCGAAGCGGGTCTGGTCTGCCCTCTCCGGGCTGGTCTGGACGCCGGTGGCGTTGAAGCGGATAAGATCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGTCAGACGAAGAGAAACAAAATCAGCTGGAACAgctggc                                                                                                                                             <  2:959200/149‑1 (MQ=255)
        taTCCGCGAAGCGGGTCTGTTCAGCCGGATCCGTGCTGTGCTGGGAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAgctggc                                                                                                                                             <  1:1245022/149‑1 (MQ=255)
         aTCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgatt                                                                                                                                            <  1:1687377/150‑4 (MQ=255)
         aTCCTCGAAGCGGGTCGGTGCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTTAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgatt                                                                                                                                            <  1:1321045/150‑4 (MQ=255)
          tCCGCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgattt                                                                                                                                           <  1:1705306/150‑5 (MQ=255)
           ccTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgatttt                                                                                                                                          <  2:594659/150‑6 (MQ=255)
           ccTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTTCCAgctgatttt                                                                                                                                          <  2:954770/150‑6 (MQ=255)
            cTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGTGAAACAAAATCAGCTGGAACAGCTGGCTGAg                                                                                                                                         >  2:376976/1‑150 (MQ=255)
            cTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAg                                                                                                                                         <  1:1310818/150‑1 (MQ=255)
            cTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAg                                                                                                                                         >  2:1109581/1‑150 (MQ=255)
             cccaaGCGCGTCCGATCCGCCGTAGCCGTGCGGGGCGGGAAGCCGGTGGCGTTGACGCTGCGAAGCTCGACACACTGCCGCGCGACCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGCACAGCTGGCTGAGc                                                                                                                                        <  2:1642123/147‑1 (MQ=255)
                   cGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTATCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATg                                                                                                                                  <  2:106044/150‑1 (MQ=255)
                     ggTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGAc                                                                                                                                >  1:101169/1‑150 (MQ=255)
                     ggTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTCGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGAc                                                                                                                                <  2:50078/150‑1 (MQ=255)
                       tCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGa                                                                                                                              <  1:1078706/150‑1 (MQ=255)
                          gTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGTAACAGCTGGCTGAGCAGTATGACGAACt                                                                                                                           <  2:789440/150‑1 (MQ=255)
                              aGCCGTATCCGTGATGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGaaa                                                                                                                       <  1:336377/150‑1 (MQ=255)
                               gCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAac                                                                                                                      <  2:463673/150‑1 (MQ=255)
                               gCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAac                                                                                                                      <  2:268511/150‑1 (MQ=255)
                                ccGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAaca                                                                                                                     >  2:1713307/1‑150 (MQ=255)
                                  gTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGAAAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACg                                                                                                                   <  1:841543/150‑1 (MQ=255)
                                    aTCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCAGGAGCTGTGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGATGGCTGAGCAGAATGACGAACTGAAACACGAg                                                                                                                 >  2:704781/1‑150 (MQ=255)
                                        gtgctgtgctGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCg                                                                                                             <  1:1052790/150‑1 (MQ=255)
                                        gggccgtgctgGTAGCCGGTGGCGTTGTCGCTGAGAAGCAGGACAAACTGCCGCTCGATCGCTGGCTGGTGCTGGGCCTGAGAGACGAAGAGAAACAAAATCAGCTGGATCAGCGGGCTGAGCAGTATGACGAACTGAAACACGAGTTCg                                                                                                             <  1:98023/145‑1 (MQ=255)
                                          gctgtgctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCgag                                                                                                           <  2:67027/150‑1 (MQ=255)
                                          gctgtgctgGTAGACGGTGGCATTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGATGGAGCTGGGACTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACTAACTGACACACGATATCgag                                                                                                           >  1:58023/1‑150 (MQ=255)
                                             gtgctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGagaag                                                                                                        >  1:1437280/1‑150 (MQ=255)
                                              tgctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGATGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGagaaga                                                                                                       >  2:629037/1‑150 (MQ=255)
                                                ctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGaaa                                                                                                     >  2:1474953/1‑150 (MQ=255)
                                                ctgGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGaaa                                                                                                     >  1:41076/1‑150 (MQ=255)
                                                  ggTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACt                                                                                                   <  2:151388/150‑1 (MQ=255)
                                                  ggTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACt                                                                                                   >  1:1691125/1‑150 (MQ=255)
                                                       ccGGTGTCGTTGAATCTGAGAAGCTCGACAAACTTCCGCGCGATCTCTGGCTGGAGCTGGGCCTGACAGACGAAGAGATACAAAATCAGCTGGATCAGCTGGCTGAGCAGTATGACTAACTGAAACACGAGTTCGAGAAGAAACTCGAAg                                                                                              <  2:506985/150‑1 (MQ=255)
                                                       ccGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAg                                                                                              <  1:1144741/150‑1 (MQ=255)
                                                         ggTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACATACTAAGAGAAACAAAATCAGCTGGACCAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGct                                                                                            >  2:841141/1‑149 (MQ=255)
                                                            ggCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGaaa                                                                                         >  1:725415/1‑150 (MQ=255)
                                                              cGTTGTAGCTGAGAAGCTCGACAAACTGCAGCGCGATCGATGGCTGGAGCTGGGCCTTACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTCAGCAGTATGACGCACTGAAACACGATTTCGAGCAGAAACTCGAAGCGAAAcg                                                                                       >  1:1565559/1‑150 (MQ=255)
                                                              cGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAAcg                                                                                       >  2:856588/1‑150 (MQ=255)
                                                               gTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGtt                                                                                                               <  1:277215/125‑1 (MQ=255)
                                                               gTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGtt                                                                                                               >  2:277215/1‑125 (MQ=255)
                                                                  caattaaagaAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAAcgccgc                                                                                   <  1:1038639/143‑1 (MQ=255)
                                                                     gCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAg                                                                                                                 <  1:288967/117‑1 (MQ=255)
                                                                     gCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAg                                                                                                                 >  2:288967/1‑117 (MQ=255)
                                                                      cTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGTAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGAAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCTAAACGCCGAaaaa                                                                               >  2:684270/1‑150 (MQ=255)
                                                                       tGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAAt                                                                              >  2:1695546/1‑150 (MQ=255)
                                                                          taaGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCAc                                                                           <  1:647218/149‑1 (MQ=255)
                                                                            aGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACTGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGAAAAATCAccc                                                                         >  1:162011/1‑150 (MQ=255)
                                                                            aGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCAccc                                                                         <  2:391757/150‑1 (MQ=255)
                                                                              cTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAg                                                                       <  1:608583/150‑1 (MQ=255)
                                                                              cTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGAAAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAg                                                                       >  1:1491126/1‑150 (MQ=255)
                                                                                cGACAAACTGCCGCGCGATCGCTGGCTGGCGCTGGGCCTGACAGCCGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTACGACGAACTGAAACACGAGTTCGAGAAGAAACTCGCAGCGAAACGCCGCAAAATCACCCAggg                                                                     >  2:277928/1‑150 (MQ=255)
                                                                                  aCAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAACTCAGCTGGCACAGCTGGCTGAGCAGTATGACGATCTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGcg                                                                   >  1:83848/1‑150 (MQ=255)
                                                                                     aaCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGcgacg                                                                >  1:1214669/1‑150 (MQ=255)
                                                                                     aaCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGcgacg                                                                >  1:502100/1‑150 (MQ=255)
                                                                                     aaCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGcgacg                                                                <  2:458349/150‑1 (MQ=255)
                                                                                        tGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATc                                                             <  1:1403706/150‑1 (MQ=255)
                                                                                         gCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACCCGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGCCGATCt                                                            >  2:904418/1‑150 (MQ=255)
                                                                                          ccGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTg                                                           <  1:1109581/150‑1 (MQ=255)
                                                                                                gATCGCTGGCTGGAGCTAGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCg                                                     >  2:1471640/1‑150 (MQ=255)
                                                                                                 aTCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCgg                                                    >  2:365142/1‑150 (MQ=255)
                                                                                                  tCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCggg                                                   <  2:976891/150‑1 (MQ=255)
                                                                                                       ggctggAGCTGGGCCTGACAGACGAAGAGAAACAAATTCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGc                                              <  2:824020/150‑1 (MQ=255)
                                                                                                        gctggAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCt                                             <  1:628513/150‑1 (MQ=255)
                                                                                                         ctggAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTg                                            <  2:1437280/150‑1 (MQ=255)
                                                                                                         ctggAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTg                                            >  1:1428762/1‑150 (MQ=255)
                                                                                                            gAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAg                                         <  2:1261496/150‑1 (MQ=255)
                                                                                                              gCTGGGCCTTACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCATGGCGACGATCTGGCACCGGGCGTGCTGAAGAt                                       <  1:1124846/150‑1 (MQ=255)
                                                                                                                tGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTg                                     >  2:474875/1‑150 (MQ=255)
                                                                                                                 gggCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGCGCGCCGATCTGGCACCCGGAGTGCTGAAGATTGt                                    >  1:979984/1‑150 (MQ=255)
                                                                                                                   gCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAACCTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTa                                  >  1:1293876/1‑150 (MQ=255)
                                                                                                                     gtGACAGTCGAAGAGAAACAAAATCAGCTGGACCAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGTAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAg                                <  1:1063125/149‑1 (MQ=255)
                                                                                                                      tGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAgg                               <  2:544516/150‑1 (MQ=255)
                                                                                                                           gACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTatat                          >  2:1283141/1‑150 (MQ=255)
                                                                                                                            aCGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTAAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATc                         >  1:45945/1‑150 (MQ=255)
                                                                                                                               aaGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACGGAAACACGAGTTCGAGAAGAAACTCGAAGCGAACCGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTgg                      >  2:102418/1‑150 (MQ=255)
                                                                                                                                agagAAACAAAATCAGCTGGCACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGAGCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGc                     <  1:563527/150‑1 (MQ=255)
                                                                                                                                  agaAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGATCTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCgg                   <  1:230188/150‑1 (MQ=255)
                                                                                                                                   gaAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGCAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGt                  >  2:1288709/1‑150 (MQ=255)
                                                                                                                                   gaAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGt                  >  1:1429950/1‑150 (MQ=255)
                                                                                                                                   gaAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGt                  >  1:1586095/1‑150 (MQ=255)
                                                                                                                                   gaAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGt                  >  2:1295512/1‑150 (MQ=255)
                                                                                                                                     aaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGACCCTCGAAGCGAAACGCCGCAAAATCACACAGGGCGACGAGCTGGCACCGGGAGTGCTGAACATTGTTAACGTATATCTGGCGGTTa                >  1:1221743/1‑150 (MQ=255)
                                                                                                                                     aaCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGCCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTa                >  1:191903/1‑150 (MQ=255)
                                                                                                                                      aCAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGCCGAACAGACACACGAGTTCGAGAAGAAACTCGAAGCTAAACGACGCAAAATCACCCAGGGCGCCGAGCTGGCACCGGGCGGGCTGTAGTTGGTTAAGGTATATCTGTCTGttca               >  2:1082925/1‑148 (MQ=255)
                                                                                                                                        aaaaTCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAAc             >  1:998581/1‑150 (MQ=255)
                                                                                                                                            tctgCTGGAACATCTGGCTGATCAGTATGAAGAACTGAAACACGAGTTCGATAAGAAACTCGAAGCGAAACGCCGCAATATAACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCg         <  1:684270/147‑1 (MQ=255)
                                                                                                                                            tCAGCTGGAACAGCGGGCTGAGCATTATGACTAACTGAAACCCGAGTGCGAGAAGAAACTCGAAGCGCAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCGGAAGATTGTTAAGGTATAGCTGTCGGGTAAACGCCg         <  1:260892/150‑1 (MQ=255)
                                                                                                                                             cAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGt        <  2:725415/150‑1 (MQ=255)
                                                                                                                                              ggCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGTCGACTATCTGGCAACGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACCCCGTa       >  1:229840/2‑150 (MQ=255)
                                                                                                                                              aGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCCCCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTa       >  2:369211/1‑150 (MQ=255)
                                                                                                                                              aGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCgaa       >  1:229837/1‑148 (MQ=255)
                                                                                                                                                cTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATc     >  2:206055/1‑150 (MQ=255)
                                                                                                                                                   gAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAg  >  2:807644/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAACTGCAGATCCTCGAAGCGGGTCTGTTCAGCCGTATCCGTGCTGTGCTGGTAGCCGGTGGCGTTGAAGCTGAGAAGCTCGACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACAGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAG  >  NC_000913/4184131‑4184427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑