Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,904,212 G→A 100% intergenic (+838/+535) ygcE → / ← queE putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,904,2120GA97.2% 1271.1 / 13.8 371intergenic (+838/+535)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (7/3);  new base A (151/203);  total (162/209)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 378 total aligned reads displayed.

GTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCTTTT  >  NC_000913/2904065‑2904358
                                                                                                                                                   |                                                                                                                                                  
gTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGACTACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGacg                                                                                                                                                  >  1:3706952/1‑150 (MQ=255)
 tCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGCATACAAAGCCGTGTATCTGCTCGGTTGATCCCCGCTGacgc                                                                                                                                                 >  1:935479/1‑150 (MQ=255)
   aGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCgg                                                                                                                                               <  2:427404/150‑1 (MQ=255)
     gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCgggg                                                                                                                                             <  2:1889682/150‑1 (MQ=255)
       tCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGaa                                                                                                                                           >  2:3843776/1‑150 (MQ=255)
       tCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCGGGCGCGGGGAACACGGGCGCACGGAAAACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGaa                                                                                                                                           >  1:2052250/1‑150 (MQ=255)
        ccGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAac                                                                                                                                          <  1:690988/150‑1 (MQ=255)
         cGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAaca                                                                                                                                         <  2:982452/150‑1 (MQ=255)
         cGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATAGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGAGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAaca                                                                                                                                         >  2:2016161/1‑150 (MQ=255)
          ggTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAacac                                                                                                                                        >  1:4062806/1‑150 (MQ=255)
           gTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCACATCTATGGACTATTGGTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACt                                                                                                                                       <  2:612688/150‑1 (MQ=255)
               aTCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAAGATATGGACTAGTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCTGTTTATCCCCGCTGACGCGGGGAACACTGGcg                                                                                                                                   >  2:672440/1‑149 (MQ=255)
                tCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTTCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGctc                                                                                                                                  <  1:3043055/150‑1 (MQ=255)
                tCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGctc                                                                                                                                  >  1:2102772/1‑150 (MQ=255)
                 ccccGCTGGCGCGGGGAACACATGTAGACCCAAATCAATGGCGTAGTGGTATTCGGTTTATACCCGCTGGCGCGGGGAACACGGGCGCGCGGAATCAAAAGCAGTGTATCGGCGCGGTGTGTCCCCGGTGTCGCGGGGGGCGCTGGctct                                                                                                                                 >  2:257308/1‑150 (MQ=255)
                    cGCTGGCGCGGGGAACACCTGTTTTCGCCCATCTATGGACTATTGCGAGTCGGGGTATCCCCGCTGGCGCGGGGACCCCGGGCGCCCGGAATACCAAGCCGTGGAGCTGCTCGGTTTCTCCCCGCTGACGCGGGGCACACTGGCGCTGCa                                                                                                                              >  2:1186223/1‑150 (MQ=255)
                     gCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCaa                                                                                                                             <  1:1431651/150‑1 (MQ=255)
                     gCTGGCGCGGGGAACACCTGTTTTCGCAAATCGATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCaa                                                                                                                             <  1:2325530/150‑1 (MQ=255)
                         gcgcGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAcagc                                                                                                                         <  1:1504307/150‑1 (MQ=255)
                          cgcgGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAcagca                                                                                                                        <  1:375041/150‑1 (MQ=255)
                              gggAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCAc                                                                                                                    <  2:1662097/150‑1 (MQ=255)
                                gAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCTccc                                                                                                                  >  2:1323979/1‑150 (MQ=255)
                                gAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCAccc                                                                                                                  <  2:2480209/150‑1 (MQ=255)
                                    acCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGa                                                                                                              >  1:2531436/1‑150 (MQ=255)
                                    acCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGAGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGTCTCTGCAACAGCAGCACCCATGa                                                                                                              >  2:4014170/1‑150 (MQ=255)
                                      cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGAcc                                                                                                            >  1:1048674/1‑150 (MQ=255)
                                         ttttCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACg                                                                                                         >  2:727713/1‑150 (MQ=255)
                                          tttCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGt                                                                                                        <  2:3148874/150‑1 (MQ=255)
                                            tCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAATCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACAGGCACGTAGACAGAAGCACCCAAGACAacgacg                                                                                                      >  1:1996477/1‑150 (MQ=255)
                                              gCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGt                                                                                                    <  2:318801/150‑1 (MQ=255)
                                              gCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGt                                                                                                    <  1:3625372/150‑1 (MQ=255)
                                                aaaTCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGttt                                                                                                  <  1:623652/150‑1 (MQ=255)
                                                   tCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATc                                                                                               <  2:839648/150‑1 (MQ=255)
                                                     tATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATccc                                                                                             <  2:1992279/150‑1 (MQ=255)
                                                       tGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCg                                                                                           >  2:3087186/1‑150 (MQ=255)
                                                         gACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCt                                                                                         >  2:1234054/1‑150 (MQ=255)
                                                             aTTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGgcg                                                                                     >  1:4006480/1‑150 (MQ=255)
                                                              ttGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGgcgc                                                                                    >  1:3921422/1‑150 (MQ=255)
                                                                gCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCgg                                                                                  <  2:1864141/150‑1 (MQ=255)
                                                                 cTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCggg                                                                                 >  1:2463586/1‑150 (MQ=255)
                                                                 cTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCggg                                                                                 <  2:1001205/150‑1 (MQ=255)
                                                                 cTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCggg                                                                                 <  2:2276777/150‑1 (MQ=255)
                                                                 cTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCACTCCAACAGCAGCACCCCTGACCACGTCGGTTTATCCCCGCTGGCGCggg                                                                                 >  2:3316568/1‑150 (MQ=255)
                                                                      cGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAaca                                                                            <  1:2115126/150‑1 (MQ=255)
                                                                      cGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAaca                                                                            <  1:3804614/150‑1 (MQ=255)
                                                                             tCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATg                                                                     >  1:2902782/1‑150 (MQ=255)
                                                                               cccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCt                                                                   >  1:2632366/1‑150 (MQ=255)
                                                                               cccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCAGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCAGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCt                                                                   >  1:1103103/1‑150 (MQ=255)
                                                                                ccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                                  >  1:3111413/1‑150 (MQ=255)
                                                                                ccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                                  >  1:1232806/1‑150 (MQ=255)
                                                                                ccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                                  <  1:413925/150‑1 (MQ=255)
                                                                                ccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                                  >  1:3494905/1‑150 (MQ=255)
                                                                                  gCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGt                                                                <  1:71390/150‑1 (MQ=255)
                                                                                  gCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGt                                                                <  2:3550068/150‑1 (MQ=255)
                                                                                  gCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGt                                                                >  1:1745659/1‑150 (MQ=255)
                                                                                       cgcgGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCg                                                           >  2:2743107/1‑150 (MQ=255)
                                                                                        gcgGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGt                                                          >  2:3735046/1‑150 (MQ=255)
                                                                                          ggggAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTa                                                        >  2:3503302/1‑150 (MQ=255)
                                                                                            ggAACACGTGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAAt                                                      <  2:190635/150‑1 (MQ=255)
                                                                                            ggAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAAt                                                      <  2:3408806/150‑1 (MQ=255)
                                                                                              aaCACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGc                                                    >  2:1781070/1‑150 (MQ=255)
                                                                                               acacGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCTTGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGcc                                                   <  2:698696/150‑1 (MQ=255)
                                                                                               acacGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGcc                                                   <  2:1048674/150‑1 (MQ=255)
                                                                                                cacGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCg                                                  <  1:1851993/150‑1 (MQ=255)
                                                                                                 acGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGc                                                 <  2:1528080/150‑1 (MQ=255)
                                                                                                    ggCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCaaa                                              <  1:2405533/150‑1 (MQ=255)
                                                                                                    cgcgcACGGAATACAAAGCCGAGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACACATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCaaa                                              <  2:763032/149‑1 (MQ=255)
                                                                                                        cACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAAcaca                                          >  2:3255040/1‑150 (MQ=255)
                                                                                                         aCGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAg                                         <  1:1561011/150‑1 (MQ=255)
                                                                                                         aCGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAg                                         >  1:1143030/1‑150 (MQ=255)
                                                                                                          cGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCACGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAgg                                        <  1:3230430/150‑1 (MQ=255)
                                                                                                            gAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGtt                                      >  1:1714902/1‑150 (MQ=255)
                                                                                                             aaTACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGttt                                     <  1:3055903/150‑1 (MQ=255)
                                                                                                             aaTACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGttt                                     <  2:2225238/150‑1 (MQ=255)
                                                                                                               tACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTAt                                   >  1:46546/1‑150 (MQ=255)
                                                                                                                aCAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCTGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATc                                  <  2:3710828/150‑1 (MQ=255)
                                                                                                                 cAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATcc                                 <  1:990210/150‑1 (MQ=255)
                                                                                                                 cAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATcc                                 <  2:539095/150‑1 (MQ=255)
                                                                                                                  aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATccc                                >  2:2078263/1‑150 (MQ=255)
                                                                                                                  aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATccc                                >  2:3614821/1‑150 (MQ=255)
                                                                                                                       cGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTg                           >  2:3931667/1‑150 (MQ=255)
                                                                                                                       cGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTg                           <  2:3439021/150‑1 (MQ=255)
                                                                                                                        gtgtATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCAgg                          >  2:1337896/1‑150 (MQ=255)
                                                                                                                            aTCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACATCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGcgcg                      <  2:2102772/150‑1 (MQ=255)
                                                                                                                            aTCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGcgcg                      <  2:3715181/150‑1 (MQ=255)
                                                                                                                            aTCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGcgcg                      <  2:2671284/150‑1 (MQ=255)
                                                                                                                             tCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCgg                     <  2:2307171/150‑1 (MQ=255)
                                                                                                                             tCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCgg                     <  2:3961368/150‑1 (MQ=255)
                                                                                                                              cTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCggg                    <  1:1878734/150‑1 (MQ=255)
                                                                                                                               tGCTCGGTTCATCCCCGCTGACGAGGGGAACTCTGTCTCTGCTACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCgggg                   <  1:1669248/150‑1 (MQ=255)
                                                                                                                                  tCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAac                <  1:1859577/150‑1 (MQ=255)
                                                                                                                                   cGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAaca               <  1:3770958/150‑1 (MQ=255)
                                                                                                                                     gTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACTGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACAGGAAATGCTGGTGTGCGTTAATGCCGCAATCACAGGTTTATCGCCGCTGGCGCGGGGGacagg             >  2:2400951/1‑148 (MQ=255)
                                                                                                                                        tATCCCCGCTGTCGCGGGGAACTCTGGCTGTGCTACAGCCGGACCGTTGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTTATGCGGCACTCACAGGTTTATCCCCGCTGGCGCGGGGAACACATTa          <  1:2719838/150‑1 (MQ=255)
                                                                                                                                         aTCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTAc         >  1:3540203/1‑150 (MQ=255)
                                                                                                                                             ccGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCt     >  1:2596383/1‑150 (MQ=255)
                                                                                                                                             ccGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCt     >  2:29710/1‑150 (MQ=255)
                                                                                                                                             ccGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCt     >  1:130011/1‑150 (MQ=255)
                                                                                                                                             ccGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCt     >  1:2993152/1‑150 (MQ=255)
                                                                                                                                                cTGACGCGGGGAACACTGGCTCTTCAACAGCAGCACCCATGCCCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCTCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCtttt  <  1:3846529/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCTTTT  >  NC_000913/2904065‑2904358

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