Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,007,816 G→A 5.1% intergenic (‑216/‑28) ycbX ← / → rlmL 6‑N‑hydroxylaminopurine detoxification oxidoreductase/23S rRNA m(2)G2445 and m(7)G2069 methyltransferases, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,007,8160GA 5.1% 667.5 / 16.1 196intergenic (‑216/‑28)ycbX/rlmL6‑N‑hydroxylaminopurine detoxification oxidoreductase/23S rRNA m(2)G2445 and m(7)G2069 methyltransferases, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (80/106);  new base A (5/5);  total (85/111)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 198 total aligned reads displayed.

AACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCAAGGGC  >  NC_000913/1007669‑1007960
                                                                                                                                                   |                                                                                                                                                
aacaaTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCgg                                                                                                                                                 >  4:1174659/1‑149 (MQ=255)
        tCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGctct                                                                                                                                         >  3:688023/1‑149 (MQ=255)
         cTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCtt                                                                                                                                        <  1:171814/149‑1 (MQ=255)
         cTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACCGCAGTAGGcgtaa                                                                                                                                        <  3:766873/149‑5 (MQ=255)
          ttAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTa                                                                                                                                       <  2:1016349/149‑1 (MQ=255)
          ttAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTa                                                                                                                                       <  4:544048/149‑1 (MQ=255)
          ttAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTa                                                                                                                                       <  1:575393/149‑1 (MQ=255)
          ttAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTa                                                                                                                                       <  4:688023/149‑1 (MQ=255)
          ttAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTa                                                                                                                                       <  3:172349/149‑1 (MQ=255)
          ttAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACCGCAGTAGGcgtaat                                                                                                                                       <  3:1046045/149‑6 (MQ=255)
              cgcCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGt                                                                                                                                   >  2:1364661/1‑149 (MQ=255)
                 cTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTaaa                                                                                                                                <  1:465579/149‑1 (MQ=255)
                 cTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTaaa                                                                                                                                <  2:812337/149‑1 (MQ=255)
                  tGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTaaaa                                                                                                                               <  1:843522/149‑1 (MQ=255)
                  tGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTaaaa                                                                                                                               <  2:485269/149‑1 (MQ=255)
                      ccATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGa                                                                                                                           >  3:323589/1‑149 (MQ=255)
                      ccATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGa                                                                                                                           >  2:602833/1‑149 (MQ=255)
                       cATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGAc                                                                                                                          >  4:350618/1‑149 (MQ=255)
                       cATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGAc                                                                                                                          <  2:1441113/149‑1 (MQ=255)
                             gCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATg                                                                                                                    >  4:415956/1‑149 (MQ=255)
                               tattttatttGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGaa                                                                                                                  <  2:773531/149‑1 (MQ=255)
                                    tatttGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATtctc                                                                                                             <  3:1115276/149‑1 (MQ=255)
                                     atttGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTctct                                                                                                            <  2:1019420/149‑1 (MQ=255)
                                     atttGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTctct                                                                                                            <  2:33007/149‑1 (MQ=255)
                                       ttGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGt                                                                                                          <  1:1153185/149‑1 (MQ=255)
                                        tGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGtt                                                                                                         <  4:493352/149‑1 (MQ=255)
                                         gTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGttt                                                                                                        >  2:729494/1‑149 (MQ=255)
                                         gTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGttt                                                                                                        <  4:784801/149‑1 (MQ=255)
                                          tCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTg                                                                                                       >  2:249810/1‑149 (MQ=255)
                                             gTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTaga                                                                                                                                     <  3:120409/116‑1 (MQ=255)
                                             gTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCa                                                                                                    <  2:1291645/149‑1 (MQ=255)
                                             gTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCa                                                                                                    <  4:627605/149‑1 (MQ=255)
                                                tGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTa                                                                                                 <  4:1258977/149‑1 (MQ=255)
                                                                tCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGt                                                                                                          <  2:154091/124‑1 (MQ=255)
                                                                tCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCa                                                                                                    >  1:988604/1‑130 (MQ=255)
                                                                tCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGa                                                                                 >  1:1257904/1‑149 (MQ=255)
                                                                tCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGa                                                                                 >  1:245529/1‑149 (MQ=255)
                                                                tCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGa                                                                                 <  4:206574/149‑1 (MQ=255)
                                                                tCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGa                                                                                 <  3:454708/149‑1 (MQ=255)
                                                                tCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGa                                                                                 >  3:410790/1‑149 (MQ=255)
                                                                  gAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAag                                                                               >  1:210862/1‑149 (MQ=255)
                                                                  gAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAag                                                                               <  3:7864/149‑1 (MQ=255)
                                                                   aaaTCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAaga                                                                              <  3:159002/149‑1 (MQ=255)
                                                                     aTCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGc                                                                            <  4:1026139/149‑1 (MQ=255)
                                                                      tCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCt                                                                           <  2:797438/149‑1 (MQ=255)
                                                                       ccTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTg                                                                          >  4:889122/1‑149 (MQ=255)
                                                                        cTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGt                                                                         <  2:425523/149‑1 (MQ=255)
                                                                        cTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGt                                                                         <  4:153973/149‑1 (MQ=255)
                                                                         tCACGTAGTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGtt                                                                        >  4:1236147/1‑149 (MQ=255)
                                                                         tCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGtt                                                                        >  3:483147/1‑149 (MQ=255)
                                                                              tACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGatat                                                                                                                     <  4:748898/99‑1 (MQ=255)
                                                                              tACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGatat                                                                                                                     >  3:748898/1‑99 (MQ=255)
                                                                                 ttGTGTACGCTCCGGTTTCTCCGCGCTGTCCGTGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTAAGAGTAAAATTGACGATATGAATTCTCTGTATGCCAGTACGGACCGTGGGTTGGAAGAGCTGTTAAAaactg                                                                >  3:899643/1‑149 (MQ=37)
                                                                                   gtgtACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAaactgaa                                                              <  1:602833/149‑1 (MQ=255)
                                                                                    tgtACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGATGAGCTGTTAAAaactgaac                                                             >  4:282795/1‑149 (MQ=255)
                                                                                    tgtACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAaactgaac                                                             <  3:706676/149‑1 (MQ=255)
                                                                                           tCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGaaaa                                                      >  2:350106/1‑149 (MQ=255)
                                                                                              ggTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTTCTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCt                                                   <  3:884632/149‑1 (MQ=255)
                                                                                               gTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTg                                                  <  1:249810/149‑1 (MQ=255)
                                                                                                tttCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTgg                                                 >  2:1259925/1‑149 (MQ=255)
                                                                                                 ttCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTggg                                                <  3:658917/149‑1 (MQ=255)
                                                                                                    tcCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAGGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGc                                             >  4:1317589/1‑149 (MQ=255)
                                                                                                    tcCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAaact                                                                 >  2:1333746/1‑129 (MQ=255)
                                                                                                    tcCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTg                                                  >  2:391738/1‑144 (MQ=255)
                                                                                                    tcCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTg                                                  <  1:391738/144‑1 (MQ=255)
                                                                                                    tcCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGc                                             >  4:1146625/1‑149 (MQ=255)
                                                                                                    tcCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGc                                             >  3:1039517/1‑149 (MQ=255)
                                                                                                    tcCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGc                                             >  4:268393/1‑149 (MQ=255)
                                                                                                         gcTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTg                                        <  2:219767/149‑1 (MQ=255)
                                                                                                         gcTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTg                                        <  4:1105788/149‑1 (MQ=255)
                                                                                                           tgtccatgtccaGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGaa                                      <  2:105023/149‑1 (MQ=255)
                                                                                                           tgtccatgtccaGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGaa                                      <  2:530057/149‑1 (MQ=255)
                                                                                                            gtccatgtccaGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAAt                                     <  3:1023884/149‑1 (MQ=255)
                                                                                                             tccatgtccaGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAACCTGAACTGGAAAACCTGGGGGCCGTTGAATg                                    <  2:1121024/149‑1 (MQ=255)
                                                                                                                    ccaGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAggtg                             >  1:677417/1‑149 (MQ=255)
                                                                                                                    ccaGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTCAATGCCAggtg                             <  4:1070367/149‑1 (MQ=255)
                                                                                                                       gactgactGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTGGAATGCCAggtggct                          >  1:1453534/1‑147 (MQ=255)
                                                                                                                         ctgactgCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCa                        >  1:299987/1‑149 (MQ=255)
                                                                                                                         ctgactgCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAACGCCAGGTGGTTCa                        <  4:352458/149‑1 (MQ=255)
                                                                                                                         ctgactgCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATACGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCa                        <  1:101427/149‑1 (MQ=255)
                                                                                                                          tgactgCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAg                       >  1:1271956/1‑149 (MQ=255)
                                                                                                                             ctgCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGt                    <  1:808448/149‑1 (MQ=255)
                                                                                                                              tgCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGTATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTg                   >  4:921300/1‑149 (MQ=255)
                                                                                                                                caacaaTTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTggg                 <  3:1423745/149‑1 (MQ=255)
                                                                                                                                 aacaaTTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTgggg                <  2:69469/149‑1 (MQ=255)
                                                                                                                                 aacaaTTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTgggg                >  2:1202312/1‑149 (MQ=255)
                                                                                                                                     taagagCCTACCGCAGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCa            >  3:355445/6‑149 (MQ=255)
                                                                                                                                     taagagCCTACCGCAGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCa            >  1:1214581/6‑149 (MQ=255)
                                                                                                                                      ttACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCAt           <  1:759060/149‑1 (MQ=255)
                                                                                                                                          gCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCa       >  3:1373164/1‑149 (MQ=255)
                                                                                                                                           ccTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCaa      <  1:1017960/149‑1 (MQ=255)
                                                                                                                                           ccTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCaa      <  4:390283/149‑1 (MQ=255)
                                                                                                                                           ccTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCaa      <  2:660497/149‑1 (MQ=255)
                                                                                                                                            cTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGATGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCAAg     <  3:282795/149‑1 (MQ=255)
                                                                                                                                            cTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCAAg     <  3:257735/149‑1 (MQ=255)
                                                                                                                                            cTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCAAg     >  1:931847/1‑149 (MQ=255)
                                                                                                                                             tACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAAt                                     <  3:800168/116‑1 (MQ=255)
                                                                                                                                             tACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAAt                                     >  4:800168/1‑116 (MQ=255)
                                                                                                                                              aCTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCAAggg   <  1:1364661/149‑1 (MQ=255)
                                                                                                                                               cTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCAAGGGc  <  4:958977/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCGCTGTCCATGTCCAGACTGACTGCAACAATTACGCCTACTGCGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCAAGGGC  >  NC_000913/1007669‑1007960

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: