Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,162,150 T→G 5.5% intergenic (+20/+20) yijD → / ← trmA DUF1422 family inner membrane protein/tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,162,1500TG 5.5% 735.2 / 22.6 199intergenic (+20/+20)yijD/trmADUF1422 family inner membrane protein/tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (98/90);  new base G (3/8);  total (101/98)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 207 total aligned reads displayed.

GTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTA  >  NC_000913/4162001‑4162299
                                                                                                                                                     |                                                                                                                                                     
gTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCt                                                                                                                                                       <  2:225707/150‑1 (MQ=255)
   cAGCACGGTAGTTCGTGCGGACTATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGgc                                                                                                                                                    >  2:202346/1‑150 (MQ=255)
   cAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAAATTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTCAGCGGTATTATGCCGGATGATGgc                                                                                                                                                    >  2:1216053/1‑150 (MQ=255)
      cACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCa                                                                                                                                                 >  1:161177/1‑150 (MQ=255)
         ggTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATcc                                                                                                                                              <  1:355648/150‑1 (MQ=255)
           tAGTTCGTGCGGAATATCCGGATATCGTCTCTATCTTCTTCCCGGCAGTACTCTGCGTCATTATGGTGTTCTGGATTGTCGCGAAGATGCTTAACATTAAGAAGGAAGTTGCTGAGTAATCTTTATTATGCCGGATGCGCCAGCATCCgg                                                                                                                                            <  2:131787/150‑1 (MQ=255)
           tAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCgg                                                                                                                                            <  2:663905/150‑1 (MQ=255)
           tAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCGCCAGCATCCgg                                                                                                                                            <  2:1583264/150‑1 (MQ=255)
           tAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCGCCAGCATCCgg                                                                                                                                            <  2:1010371/150‑1 (MQ=255)
                  tGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGt                                                                                                                                     <  1:819595/150‑1 (MQ=255)
                    cGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGttt                                                                                                                                   <  1:1284695/150‑1 (MQ=255)
                      gAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTa                                                                                                                                 <  1:319907/150‑1 (MQ=255)
                      gAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTa                                                                                                                                 >  1:2450102/1‑150 (MQ=255)
                      gAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTa                                                                                                                                 <  1:441203/150‑1 (MQ=255)
                      gAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTa                                                                                                                                 <  2:533397/150‑1 (MQ=255)
                      gAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTa                                                                                                                                 >  1:609961/1‑150 (MQ=255)
                      gAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTa                                                                                                                                 >  2:1365664/1‑150 (MQ=255)
                      gAATATCCGGAGATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTa                                                                                                                                 <  2:180977/150‑1 (MQ=255)
                           tCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcg                                                                                                                            >  2:1291291/1‑150 (MQ=255)
                           tCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcg                                                                                                                            >  1:981950/1‑150 (MQ=255)
                           tCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcg                                                                                                                            >  1:2263224/1‑150 (MQ=255)
                            ccGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcgc                                                                                                                           <  2:519430/150‑1 (MQ=255)
                             cGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTcgcg                                                                                                                          <  1:1538400/150‑1 (MQ=255)
                                 tatCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCa                                                                                                                      <  2:408552/150‑1 (MQ=255)
                                 tatCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCa                                                                                                                      <  2:2003767/150‑1 (MQ=255)
                                 tatCGGCTCTAACTTCGTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCa                                                                                                                      <  1:2177186/150‑1 (MQ=255)
                                     ggCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGAGGTCAGTaa                                                                                                                  >  2:1658876/1‑150 (MQ=255)
                                       ctctAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATa                                                                                                                >  1:1528919/1‑150 (MQ=255)
                                       ctctAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATa                                                                                                                >  2:1922352/1‑150 (MQ=255)
                                        tctAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAc                                                                                                               >  1:1759344/1‑150 (MQ=255)
                                        tctAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAc                                                                                                               >  1:965710/1‑150 (MQ=255)
                                          tAACTTCGTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATAcgc                                                                                                             <  2:955691/150‑1 (MQ=255)
                                              ttcttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCa                                                                                                         <  2:863730/150‑1 (MQ=255)
                                              ttcttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGCAATACGCCGCa                                                                                                         >  2:1899509/1‑150 (MQ=255)
                                               tcttcCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCAc                                                                                                        >  1:454447/1‑150 (MQ=255)
                                                   cccGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCa                                                                                                    <  2:450882/150‑1 (MQ=255)
                                                   cccGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCa                                                                                                    >  2:2514770/1‑150 (MQ=255)
                                                   cccGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCa                                                                                                    <  1:1182942/150‑1 (MQ=255)
                                                    ccGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCat                                                                                                   >  1:1743042/1‑150 (MQ=255)
                                                           tACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGc                                                                                            >  1:2219873/1‑150 (MQ=255)
                                                                tcCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTa                                                                                       >  1:1453391/1‑150 (MQ=255)
                                                                tcCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTa                                                                                       >  2:555661/1‑150 (MQ=255)
                                                                 ccGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGCAATACGCCGCACTCCATATGGTGCGTGTAg                                                                                      <  1:1857936/150‑1 (MQ=255)
                                                                  cGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAgg                                                                                     >  2:2294268/1‑150 (MQ=255)
                                                                   gTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAggg                                                                                    >  1:1247180/1‑150 (MQ=255)
                                                                   gTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAggg                                                                                    <  2:591345/150‑1 (MQ=255)
                                                                    tCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAgggg                                                                                   >  2:430719/1‑150 (MQ=255)
                                                                    tCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAgggg                                                                                   >  2:1917220/1‑150 (MQ=255)
                                                                    tCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAgggg                                                                                   >  2:1678290/1‑150 (MQ=255)
                                                                           ggTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGAt                                                                            <  2:260851/150‑1 (MQ=255)
                                                                           ggTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGAt                                                                            <  2:47772/150‑1 (MQ=255)
                                                                             tgtTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGAtga                                                                          <  2:1062286/150‑3 (MQ=255)
                                                                               ttCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCaaa                                                                        <  1:1936058/150‑1 (MQ=255)
                                                                               ttCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCaaa                                                                        <  1:2379642/150‑1 (MQ=255)
                                                                                tCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAAc                                                                       >  1:412403/1‑150 (MQ=255)
                                                                                tCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCAGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAAc                                                                       >  1:1712923/1‑150 (MQ=255)
                                                                                 cTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACa                                                                      >  1:1023634/1‑150 (MQ=255)
                                                                                  tGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACag                                                                     >  2:2341090/1‑150 (MQ=255)
                                                                                  tGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACag                                                                     >  2:711627/1‑150 (MQ=255)
                                                                                  tGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACag                                                                     >  2:478662/1‑150 (MQ=255)
                                                                                         gcgcGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCATGTAGGGGAACTGATCAAACAGAGCCAGa                                                              <  1:2559179/150‑1 (MQ=255)
                                                                                                gATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGa                                                       <  1:1283627/150‑1 (MQ=255)
                                                                                                gATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGa                                                       <  1:1596040/150‑1 (MQ=255)
                                                                                                gATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGa                                                       <  2:1691014/150‑1 (MQ=255)
                                                                                                gATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGa                                                       <  2:1108258/150‑1 (MQ=255)
                                                                                                gATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGa                                                       <  1:1117862/150‑1 (MQ=255)
                                                                                                 aTGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGAc                                                      <  1:1915103/150‑1 (MQ=255)
                                                                                                  tGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGAcc                                                     >  2:2027145/1‑150 (MQ=255)
                                                                                                  tGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGAcc                                                     >  2:2066516/1‑150 (MQ=255)
                                                                                                   gCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCt                                                    >  2:2388717/1‑150 (MQ=255)
                                                                                                      tAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTtgt                                                 >  2:918123/1‑150 (MQ=255)
                                                                                                      tAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTtgt                                                 >  2:1159176/1‑150 (MQ=255)
                                                                                                          cGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGt                                             >  1:1043208/1‑150 (MQ=255)
                                                                                                            tAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCt                                           <  2:843893/150‑1 (MQ=255)
                                                                                                            tAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCt                                           >  2:1011571/1‑150 (MQ=255)
                                                                                                              aGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTgg                                         <  2:1456394/150‑1 (MQ=255)
                                                                                                              aGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTgg                                         <  1:2504377/150‑1 (MQ=255)
                                                                                                               gCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGc                                        <  1:1776915/150‑1 (MQ=255)
                                                                                                               gCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGc                                        <  1:1228835/150‑1 (MQ=255)
                                                                                                                 aGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCtt                                      >  2:311823/1‑150 (MQ=255)
                                                                                                                   gAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTaa                                    >  1:1886608/1‑150 (MQ=255)
                                                                                                                      gTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGt                                 <  2:1291161/150‑1 (MQ=255)
                                                                                                                      gTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGt                                 <  2:2445342/150‑1 (MQ=255)
                                                                                                                        tGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGttt                               <  1:140567/150‑1 (MQ=255)
                                                                                                                         gCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTc                              >  2:1130421/1‑150 (MQ=255)
                                                                                                                           tGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCa                            >  2:2295364/1‑150 (MQ=255)
                                                                                                                             aGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACAGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGa                          >  1:546623/1‑150 (MQ=255)
                                                                                                                                      tattatGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCAt                 <  1:1625063/150‑1 (MQ=255)
                                                                                                                                       attatGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGTGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATa                <  2:1442/150‑1 (MQ=255)
                                                                                                                                       attatGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATa                <  1:586754/150‑1 (MQ=255)
                                                                                                                                          atGCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACg             >  1:2577320/1‑150 (MQ=255)
                                                                                                                                          atGCCGGATGCGCCAGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACg             >  1:2076426/1‑150 (MQ=255)
                                                                                                                                             ccGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGttt          >  2:1162331/1‑150 (MQ=255)
                                                                                                                                             ccGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGttt          >  2:2273843/1‑150 (MQ=255)
                                                                                                                                             ccGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCAGATGGTGCGTGGAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCAGAACGgtt          >  2:2220390/1‑150 (MQ=255)
                                                                                                                                                gATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCg       >  2:2451821/1‑150 (MQ=255)
                                                                                                                                                gATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCg       >  2:2451732/1‑150 (MQ=255)
                                                                                                                                                  tGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCggg     >  1:910404/1‑150 (MQ=255)
                                                                                                                                                     tGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGtaa  >  1:1793213/1‑148 (MQ=255)
                                                                                                                                                     tGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTa  >  2:1256135/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAACCGTAAGCAGGAAGTTGCTGAGTAATCGGTATTATGCCGGATGCTGGCGCATCCGGCATGGGTTTTACTTCGCGGTCAGTAATACGCCGCACTCCATATGGTGCGTGTAGGGGAACTGATCAAACAGAGCCAGACGTTCGACCTTGTGCGTCTGGCTTAATGTTTCCAGATTCTTGCATAACGTTTCCGGGTTA  >  NC_000913/4162001‑4162299

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