Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,149 G→A 5.6% noncoding (88/171 nt)
R10* (CGA→TGA) 
rttR ←
tpr ←
rtT sRNA, processed from tyrT transcript
protamine‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,1490GA 5.6% 994.2 / 34.0 357noncoding (88/171 nt)
R10* (CGA→TGA) 
rttR
tpr
rtT sRNA, processed from tyrT transcript
protamine‑like protein
Reads supporting (aligned to +/- strand):  ref base G (156/167);  new base A (11/8);  total (175/182)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 365 total aligned reads displayed.

TTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTC  >  NC_000913/1287001‑1287298
                                                                                                                                                    |                                                                                                                                                     
ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCGCGCTCGCTGGCGCTCGa                                                                                                                                                      >  2:4933131/1‑150 (MQ=255)
  ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                                                                                                                                    <  1:1688416/150‑1 (MQ=255)
   cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACGCGCTGGCGCTCGAGTc                                                                                                                                                   >  2:8359378/1‑150 (MQ=255)
   cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCAGCGCGCTTCGCCCTTAGGGTCGTTGCCTGCGGCAACGCACTCTAGCTCGCGCTCGAGTc                                                                                                                                                   >  2:6361956/1‑150 (MQ=14)
    ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGGCGCACGAGTCg                                                                                                                                                  >  2:7148588/1‑150 (MQ=255)
     cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCAACGCTGTCTAGCTTGCGgtcgagtcga                                                                                                                                                 >  1:6897125/1‑150 (MQ=14)
     cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGa                                                                                                                                                 >  2:5218762/1‑150 (MQ=255)
        aTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCACGCTGGCGCGCGAGTCGAAcc                                                                                                                                              >  2:1778345/1‑150 (MQ=255)
         tGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCAGCCGCAACGCACTCTCGCTCGCGCTCGCGTCGCACCt                                                                                                                                             >  2:3616005/1‑150 (MQ=14)
             cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGt                                                                                                                                         >  2:1957858/1‑150 (MQ=255)
              aGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTCTCGCGGGCGCTCGAGTCGACCCAAGGTc                                                                                                                                        >  1:2788120/1‑150 (MQ=14)
                aaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCACTCTCGCATGCGCACGAGTCGAAACATGGTCGa                                                                                                                                      >  1:9365093/1‑150 (MQ=14)
                 atatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGaa                                                                                                                                     <  1:1903195/150‑1 (MQ=255)
                      tGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                                                <  1:7549826/150‑1 (MQ=255)
                         ttGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCCACGCTCTCACGCTGGCGCTCGAGTCGAACCATGGTCGAAGCTTCTCa                                                                                                                             >  1:8588869/1‑150 (MQ=255)
                           gCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCaac                                                                                                                           >  2:9092511/1‑148 (MQ=255)
                                   gCTTGAGAGTTCAGGGACTTGTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGGGGCAACGCTCACTCGCTGGCGCTCGAGTCGACACATGGTGGAAGCTTCTCATCCTTCCCCg                                                                                                                   >  2:5779332/1‑150 (MQ=14)
                                     ttGATAGTTCATGGTGTTTTGATAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                                 <  1:6733569/150‑1 (MQ=14)
                                     ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCTCCAGTCGAACCGTGGTCGAAGCTTCTCAACCATCCCCCCa                                                                                                                 >  1:8029983/1‑150 (MQ=14)
                                     gtGAGAGTTAAGGGACTTTTGAAAGTGATGTTGGTGGGGGAAGGATTACTCAGCGCTGCGCTCTTCGCCCTTCGGGACGTTGACTTAGGCAACGTTCTCTATCTATCCCTCGAATCGAACCTTAGTCGAAGCTTCTCAACCTTCCCCGct                                                                                                                 >  2:7571280/2‑149 (MQ=11)
                                             tCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                                                                                                         <  1:6000234/150‑1 (MQ=255)
                                                 ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTAGAAGCTTCTCATCCATCCCCGCAGTGGCCGAaaat                                                                                                     >  2:2327247/1‑150 (MQ=14)
                                                 ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCGAGTCGAACCATGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAatat                                                                                                     >  2:555080/1‑150 (MQ=255)
                                                 ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGCCAACGCACACTCGCTGGCGCTCGAGTCGAACCATGGTCGAAGCTTCTCAACCTTCCCCGCATGGGCAGAatat                                                                                                     >  2:548614/1‑150 (MQ=14)
                                                    cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTa                                                                                                  >  2:5901195/1‑150 (MQ=38)
                                                     ttttGACAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                                                 <  1:6598483/150‑1 (MQ=25)
                                                     ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCAACCTTCCCCGCAAGGTCAGaagagcgac                                                                                                 >  2:8712912/1‑143 (MQ=11)
                                                       ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAtt                                                                                               <  1:7047419/150‑1 (MQ=39)
                                                       ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAtt                                                                                               <  2:1610232/150‑1 (MQ=39)
                                                       atggaaGGGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTGGTCGAAGCTTCTCACCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                               <  2:4429415/146‑1 (MQ=11)
                                                         gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACACTCGCGGGCGCGCGAGTCGAAACTTTGTAGAAGCATATCATACATAACCCCAAGGGCAGAAAAATTAATTga                                                                                             >  1:8739954/1‑149 (MQ=2)
                                                               gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTAATTGCGGATTc                                                                                       <  2:9523128/150‑1 (MQ=37)
                                                                aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAgc                                                                                                                                                    <  1:8837364/88‑2 (MQ=17)
                                                                aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCg                                                                                      <  1:903491/150‑1 (MQ=37)
                                                                  ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTACCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                                                    >  2:7768058/1‑150 (MQ=37)
                                                                  ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCACTCGCAGGCGCGCGCGTCGAACCATGGTCGAAGCTTATCAAACATCCCCCCCGGGGCAGACTAATTAATAGCGGATTCggg                                                                                    >  1:5266712/1‑148 (MQ=11)
                                                                   gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTg                                                                                   <  2:8861092/150‑1 (MQ=38)
                                                                        gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAg                                                                              <  1:5400484/150‑1 (MQ=255)
                                                                           ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGaa                                                                                                                                     <  1:8863052/92‑1 (MQ=33)
                                                                           ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTGTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTc                                                                           <  1:2237380/150‑1 (MQ=255)
                                                                             aaGGATTACTGAGCGCTGCTCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                                                         <  1:7879284/150‑1 (MQ=255)
                                                                              aGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                                                        <  1:531660/150‑1 (MQ=255)
                                                                                gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                                                                                                                                    <  1:5169599/72‑1 (MQ=34)
                                                                                gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                              <  1:1251490/94‑1 (MQ=33)
                                                                                 aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGAc                                                                     <  2:3337837/150‑1 (MQ=255)
                                                                                  ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACt                                                                    <  2:531646/150‑1 (MQ=255)
                                                                                    aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACttt                                                                  <  2:3791302/150‑1 (MQ=255)
                                                                                     cacAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCACCGCTCACTCGGTGGCGCTCGAGTAGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTGGGGAAGTTCAGGGACtttt                                                                 <  1:1378808/148‑1 (MQ=14)
                                                                                        aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                                              >  2:6102010/1‑150 (MQ=255)
                                                                                         gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGTACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                                             <  1:4136385/150‑1 (MQ=255)
                                                                                         gcgcTGCGCGCCTCGCCCTTCGGGTCGTTGCCTGCGGCTACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGATGCTTCTCAACCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                                             <  2:7227157/150‑1 (MQ=255)
                                                                                          cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAGTCGTTGGGAAGTTCAGGGACGTTTGAAAg                                                            >  2:1319349/1‑150 (MQ=255)
                                                                                          cgcTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAg                                                            <  2:5659808/150‑1 (MQ=255)
                                                                                          cgcAGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGACGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAg                                                            <  2:2145251/150‑1 (MQ=255)
                                                                                           gcgcagcgcTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGt                                                           >  2:644166/6‑150 (MQ=255)
                                                                                             cgcgcgcGTCGCCCTTCGGGTCGGTGCCCGCGGCAACGCTCTCGCGCTGCCGCTCGAGTCGAACCCTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                         <  2:8296158/149‑1 (MQ=14)
                                                                                               cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtg                                                       >  2:85577/1‑150 (MQ=255)
                                                                                                 cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCCGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggt                                                     <  1:4471124/150‑1 (MQ=255)
                                                                                                       gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                               <  1:6102010/150‑1 (MQ=255)
                                                                                                        cgcctCGGGTCGTTGCCTGCTGATACGCTCTCTCGCTGTCGCTCGAGTCGCACCTTGGTAGATGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                              <  2:6849853/146‑1 (MQ=14)
                                                                                                        cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCCCATGGGCCGAAAATTTAATTGCGGAAACGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTtggg                                              >  2:6697541/1‑150 (MQ=17)
                                                                                                        cccTTCGGGTCGTTGCCTGCGGCAACGCTATCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCGTATCATCCTTCCCCCCAAGGGCAGAAAATTTAAATGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                              >  2:2903700/1‑150 (MQ=17)
                                                                                                        cccTTCGGGTCGTTGCCTGCGGCAACGCCCTCTCGCTGGCGCTCGAGTCGAACCTTGGTATAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                              <  2:4050515/150‑1 (MQ=255)
                                                                                                         ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAAAATTTAAATGCGGAACCGTTGGGAAGGTAAGGGCCTTTTGAAAGTGATTGGGGTggggg                                             >  2:8086074/1‑150 (MQ=17)
                                                                                                          cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGa                                            <  2:9482676/150‑1 (MQ=255)
                                                                                                           ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGTAGCTTCTCATCCTTCCCCGGATGTGCTGAATATTTTATTGCGGATTCGTTGGGTAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                           <  1:6926823/150‑1 (MQ=14)
                                                                                                           gtgtggTCGTGGGCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                           <  1:4963263/146‑1 (MQ=14)
                                                                                                              gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGTTGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                        <  1:6361956/150‑1 (MQ=255)
                                                                                                              gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                        <  1:8359306/150‑1 (MQ=255)
                                                                                                                gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCACATCCTTACCCGCATGGGCAGTATATTTAATTTCGGATTCGTTGGGAAGTTAAGGGACTGTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                      >  2:2249917/1‑150 (MQ=14)
                                                                                                                gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTAGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                      <  1:5779332/150‑1 (MQ=255)
                                                                                                                gTCGGTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTGCCCCGCATGGGCAGAATATTTATTTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                      <  1:8290925/150‑1 (MQ=255)
                                                                                                                 tCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc                                     <  1:7188307/150‑1 (MQ=255)
                                                                                                                   gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTg                                                                                              >  2:4832987/1‑91 (MQ=40)
                                                                                                                   gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTg                                                                                              <  1:4832987/91‑1 (MQ=40)
                                                                                                                    ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa                                  <  1:891203/150‑1 (MQ=255)
                                                                                                                    gtGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa                                  <  1:4603017/149‑1 (MQ=255)
                                                                                                                        cagCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGACCCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGTATATTTTTTTGCGGATTCGTTGGGTAGTTCAGTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCtt                              <  1:3695442/148‑1 (MQ=17)
                                                                                                                        cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                                   >  2:9373247/1‑81 (MQ=38)
                                                                                                                         tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc                             <  1:7772287/150‑1 (MQ=255)
                                                                                                                         tGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtg                                                       >  2:621105/1‑124 (MQ=255)
                                                                                                                           cGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGa                           >  2:7784157/1‑150 (MQ=255)
                                                                                                                           cGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAAACTTGGTCGAAGCTTATCATCCATCCCCCCAGGGGCCGAAAATTTAATTTCGGATTCGTTGGGAAGTTCAGGGACTTTTTAAAGTGATGGTGGTGGGGGAAGGATTCGAAccaaaga                           >  1:9526026/1‑145 (MQ=17)
                                                                                                                               aaatctctctCGCTGACGCTCGTGCCGAACCTTAGTCGATTCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTTATTGCGGATTCGTTGGTTAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTc                       <  2:2379788/146‑1 (MQ=14)
                                                                                                                                  gccggtgcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGAt                    <  2:1240073/143‑1 (MQ=255)
                                                                                                                                  gCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGAt                    >  2:4523533/1‑150 (MQ=255)
                                                                                                                                  gCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                                             >  1:9463982/1‑77 (MQ=21)
                                                                                                                                      tctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACg                >  2:3201963/1‑150 (MQ=255)
                                                                                                                                       ctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACgg               <  1:266280/150‑1 (MQ=255)
                                                                                                                                        tcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                               <  1:1461035/117‑1 (MQ=255)
                                                                                                                                        tcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTCATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                               <  2:6962552/117‑1 (MQ=255)
                                                                                                                                        tcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                               <  1:2949973/117‑1 (MQ=255)
                                                                                                                                         cGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCa             <  1:5638653/150‑1 (MQ=255)
                                                                                                                                           cTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCAACCAAACCCGCATGGGCAGAATAAATAATTGCGGATTCGGTGGGAAGTTCACGGACATTTTAAAGTGATGGTGGTGGGGGCAGGATTCGAACCATCGAAGTAGATTACAGCAGa           >  1:3809458/1‑150 (MQ=255)
                                                                                                                                             ggCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATGCTTACCCGCATGGGCAGAATATTTTATTTCGGTTTGGTTTGGTAGTTGATGGACTTTTGTAAGTGTTGGTTGTGGGGGTAGGATTGGTACCTTAGAATTAGATTTAGGCAGAtt         >  1:5042322/1‑150 (MQ=9)
                                                                                                                                              gcgcTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                                                 <  2:3514082/61‑1 (MQ=255)
                                                                                                                                                 cTCGAGTCGAACCTTGCTCGACGCTTCTCATCCTTCCACGCCTGGTCAGAATATTTAATTTCGGATTCGTTGGTCAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACa     <  1:6556355/150‑1 (MQ=21)
                                                                                                                                                    gAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTc  <  1:2429946/150‑1 (MQ=255)
                                                                                                                                                    gAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTc  <  1:1141598/150‑1 (MQ=255)
                                                                                                                                                    gAGTCGAACCTTGGTCGAAGCTTCTCAACCATCCCCCCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTAGAAGTCGATGACGGCAGATTTACAGTc  >  1:9639591/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTC  >  NC_000913/1287001‑1287298

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: