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,2120GA100.0% 236.7 / ‑8.3 98intergenic (+838/+535)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/56);  minor base C (4/1);  total (39/59)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCTTT  >  NC_000913/2904068‑2904357
                                                                                                                                                |                                                                                                                                                 
aGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACGAAGCCGTGTAGCTGCTCGGGTTATCCCCGCTGAcgcg                                                                                                                                               >  2:931786/1‑149 (MQ=255)
 ggatCCTGTTTATCCCCGCTGGCTCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCgg                                                                                                                                              <  1:312024/146‑1 (MQ=255)
 ggCTCCGGTTTATCCCCGCTGGCGCGGGGAACCCCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGATTATCACCGCTGGCGCGGGGAACACGGGCGCACGGATTACAAAGCCGTGGATCTGCTCGGTTTATCCGCGCTGACGCgg                                                                                                                                              >  1:685972/1‑149 (MQ=255)
  gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCggg                                                                                                                                             <  1:181814/149‑1 (MQ=255)
  gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCGGTTTTCGCAAATCTATGGAGTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCGGGTCGGTTTATCCCCGCTGAGGCggg                                                                                                                                             >  1:707439/1‑149 (MQ=38)
   cTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGTGGAACACGGGCGCACGGAGTACAAAGCCGTGTATCTGCGCGGTTTATCCCCGCTGACGCgggg                                                                                                                                            >  1:683190/1‑149 (MQ=255)
    tCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGCCGCGGGGa                                                                                                                                           >  2:28557/1‑149 (MQ=255)
      cGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGCac                                                                                                                                         >  1:467459/1‑149 (MQ=255)
      cGGTTAATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTCTGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAac                                                                                                                                         <  1:47185/149‑1 (MQ=255)
       ggTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAaca                                                                                                                                        <  1:112856/149‑1 (MQ=255)
       ggTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAaca                                                                                                                                        >  1:1162916/1‑149 (MQ=255)
        gTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCCCGGAAGTCAAAGCCGTGGAGCTGCTCGGTGTTGACCGGCGGACGCGGGGAacac                                                                                                                                       >  1:1135930/1‑149 (MQ=255)
         tttATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACt                                                                                                                                      >  1:951488/1‑149 (MQ=255)
          ttATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTg                                                                                                                                     <  2:1016894/149‑1 (MQ=255)
           tATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACttg                                                                                                                                    <  2:532300/149‑3 (MQ=255)
           tATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTgg                                                                                                                                    <  1:927768/149‑1 (MQ=255)
           tATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTgg                                                                                                                                    <  2:499622/149‑1 (MQ=255)
           tATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTgg                                                                                                                                    <  1:92181/149‑1 (MQ=255)
             tCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATCGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGct                                                                                                                                  <  2:231157/149‑1 (MQ=255)
               cgcGCTGGCGCGGGGTACACCTGTTTTCACAAATCTATGGACTATTGCTCTTCGGTTTATCCCCGCTGGAGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGctct                                                                                                                                <  2:219546/147‑1 (MQ=255)
                ccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTg                                                                                                                               >  1:599754/1‑149 (MQ=255)
                ccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTAGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCACGCTGACGCGGGGAACACTTGCTCTg                                                                                                                               >  1:79345/1‑149 (MQ=255)
                ccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATAGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGTGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTg                                                                                                                               <  2:409194/149‑1 (MQ=255)
                     ggCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAca                                                                                                                          <  1:587498/149‑1 (MQ=255)
                        gcgGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAcagca                                                                                                                       <  2:951488/149‑1 (MQ=255)
                           gggAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAgcagca                                                                                                                    <  1:1129767/149‑1 (MQ=255)
                              agccccTGTTTTCGCATCTTTCTGGCCTCGTGCTAGTCGGGTGATAAAAGATGGAGCGGGGCACACGTGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCAccc                                                                                                                 <  2:1186376/145‑1 (MQ=255)
                                 acCTGTTTTCGCAAATCTATGGCCTATTGCTATGCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATCCAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATg                                                                                                              <  2:1154280/149‑1 (MQ=255)
                                  ccTGTATTCGCAAAGCGATGGCCTATTGCTATTCGGGTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCGACAGCAGCACCCATGa                                                                                                             <  2:668537/149‑1 (MQ=255)
                                   cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATAAAAAGGCGTGTATCTGCTCGGTTTATCCCCCCTGACGCCGGGAACACTGCCTCTGCTACAGCGGAACCCAAGAc                                                                                                            >  2:663188/1‑149 (MQ=255)
                                    tGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGAcc                                                                                                           <  1:170465/149‑1 (MQ=255)
                                            cAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATAAAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTCACGCGGGGAACACTGGCTCTGCAACCGCAGCACCCATGACCACGTctgg                                                                                                   >  2:934635/1‑146 (MQ=255)
                                            cAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGAGGGGAACACGGGCGCACGGAATACAAAGCCGTGTAACAGCTCGGTTTATCCCCGCTGACGCGGTGAACACTGGCTCTGCAACAGCAGCCCCCATGCCCACGTCGGt                                                                                                   >  1:992771/1‑149 (MQ=255)
                                             aaaTCTATTGACTATTGCTATTCGGTTTATCCACGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCAGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACACATGACCACGTCGGtt                                                                                                  <  2:545955/149‑1 (MQ=255)
                                              aaTCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGttt                                                                                                 <  1:335692/149‑1 (MQ=255)
                                              aaTCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGGCGGttt                                                                                                 <  2:627549/149‑1 (MQ=255)
                                                  tATTGACTATTGCTATTCGGTTTATCCCAGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGGTTATcc                                                                                             <  2:635290/149‑1 (MQ=255)
                                                      gACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCTTGTATCTTCTAGGTTTATCCCCGCTGACGCGGGGAAAACTGGCTCTGCAACAGCCGCACCCATGACCACGTCGGTTTATccactc                                                                                         >  2:872569/1‑145 (MQ=255)
                                                         tATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTgg                                                                                      >  2:84545/1‑149 (MQ=255)
                                                          aTTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGgc                                                                                     <  2:1162916/149‑1 (MQ=255)
                                                             gCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGGCGGTTTATCCCCGCTGGcgcg                                                                                  <  2:887974/149‑1 (MQ=255)
                                                              cTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCgg                                                                                 <  2:599754/149‑1 (MQ=255)
                                                              cTATTCGGTTTATCCCCGCTGGAGCGGGTAACACTGTAGCACGGAATACTAAGCCGTTTATCTGCTCGTTGTATCCCCGCTGTCGCGGTGAACACTGGCTCTGCTACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCgg                                                                                 <  2:780576/149‑1 (MQ=255)
                                                              cTATTCGGTGTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCGCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCgg                                                                                 <  2:728814/149‑1 (MQ=255)
                                                               tATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCggg                                                                                <  2:350298/149‑1 (MQ=255)
                                                                 ttCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGCCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGa                                                                              >  2:496427/1‑149 (MQ=255)
                                                                  tCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGaa                                                                             <  1:936442/149‑1 (MQ=255)
                                                                  tCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGGGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCCGCACCCATGACCACGGCGGTTTATCCCCGCTGGCGCGGGGaa                                                                             <  1:107689/149‑1 (MQ=255)
                                                                     gTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAacac                                                                          <  1:28557/149‑1 (MQ=255)
                                                                       ttctCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGa                                                                        <  2:949345/146‑1 (MQ=255)
                                                                       ttATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGa                                                                        >  2:1182636/1‑149 (MQ=255)
                                                                         aTCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGaaa                                                                      <  1:668710/149‑1 (MQ=255)
                                                                            cccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGCCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGGACCCGAAATGc                                                                   >  2:610755/1‑149 (MQ=255)
                                                                              cGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                                 >  1:105974/1‑149 (MQ=255)
                                                                              cGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                                 >  1:648922/1‑149 (MQ=255)
                                                                              cGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                                 <  2:955617/149‑1 (MQ=255)
                                                                              cGCTGGCGCGGGGAACACGGGAGCTCGGAATTCAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                                 <  1:1046685/149‑1 (MQ=255)
                                                                              cGCAGGCGCGGGTAACACGGGCGCACGGAATACTAAGCCGTGTATCTGTTCGTTTTATACCAGCTGACGCGGGGAACACTGTCTCTGCAACAGCAGCACCGATGACCACTTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTg                                                                 <  2:760619/149‑1 (MQ=255)
                                                                               gCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGCCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATTCTgg                                                                >  1:12086/1‑149 (MQ=255)
                                                                               gCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTgg                                                                <  2:355393/149‑1 (MQ=255)
                                                                                cTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGt                                                               <  2:1179866/149‑1 (MQ=255)
                                                                                  ggCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGa                                                             <  1:982397/149‑1 (MQ=255)
                                                                                        gggAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGAACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTa                                                       >  2:1046651/1‑149 (MQ=255)
                                                                                          taaCACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAAt                                                     <  2:887364/148‑1 (MQ=255)
                                                                                          gAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGTCCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAAt                                                     <  1:585545/149‑1 (MQ=255)
                                                                                          gAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTAGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACGCATGACCTCGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAAt                                                     <  1:640562/149‑1 (MQ=255)
                                                                                           aaCACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATg                                                    <  1:218825/149‑1 (MQ=255)
                                                                                                gggCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCa                                               >  1:1107596/1‑149 (MQ=255)
                                                                                                  gcgcACGGAATACAAATCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCaaa                                             <  1:682519/149‑1 (MQ=255)
                                                                                                  gcgcACGGAATACAAAGGCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCaaa                                             <  2:557631/149‑1 (MQ=255)
                                                                                                    gcACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCCTGACCACGTCGGTTTATCCCCGCGGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAaca                                           >  2:1047733/1‑149 (MQ=255)
                                                                                                    gcACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAaca                                           <  2:648922/149‑1 (MQ=255)
                                                                                                       cGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAg                                        <  1:607072/149‑1 (MQ=255)
                                                                                                       cGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAg                                        <  1:571075/149‑1 (MQ=255)
                                                                                                         taaTACAAATCCGTGTATCTGCTAGGTTTATCCCCGCTGACTCGGGGATCACTGGATCTCCAACAGCAGCAACCATGACCACGTCTGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGt                                      <  1:669380/148‑1 (MQ=255)
                                                                                                           aTACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCCACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGAGTTAATGCCGCAAACACAGGttt                                    >  2:957940/1‑149 (MQ=255)
                                                                                                               aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAGATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATcc                                >  2:1092105/1‑149 (MQ=255)
                                                                                                               aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATcc                                >  2:655313/1‑149 (MQ=255)
                                                                                                               aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCAACCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCACACACAGGTTTATcc                                >  2:607700/1‑149 (MQ=255)
                                                                                                                aaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATccc                               <  2:447706/149‑1 (MQ=255)
                                                                                                                 aGCCGTGTATCTGCTCGGTTAATCCCCGCTGCCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGCCCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATcccc                              <  2:663802/149‑1 (MQ=255)
                                                                                                                   ccGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGc                            >  1:855277/1‑149 (MQ=255)
                                                                                                                   ccGTGTATCTGCTCGGGTTATCCCCGCTGACGCGGGGAACACTTGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGc                            <  1:1074721/149‑1 (MQ=255)
                                                                                                                      tgtATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTgg                         >  1:414907/1‑149 (MQ=255)
                                                                                                                         aTCTGCGCGGTTTATCCCCGCTGACGCGGGGAACACTTGCTCTGCAACAGCAGCACGCATGACAACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGgcgc                      <  1:1143352/149‑1 (MQ=255)
                                                                                                                          tctcCTCGGTTTATCCCCTCTGTCGCCGGGAAGACTGGCTCTGAAACAGTAGCTCCGATGACCTCTTCGGTTTATCCCTGCTGTCGCGGGGAACACGACATTATTGTTAGCGTTAATGTCTTATTCACAGGTTTATCCCCGCTGGcgcg                     <  2:685972/149‑1 (MQ=255)
                                                                                                                             gCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCATCAGCAGCCCCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACAGGAAATGCTGGTGAGCGTTAATGCCGCAAACCCAGGTTTATCCCCGCTGGCGCgggg                  >  1:372577/1‑149 (MQ=255)
                                                                                                                             gCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGATATTCTTGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCgggg                  >  2:804477/1‑149 (MQ=255)
                                                                                                                              cTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCCCGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTACTGCCGCAAACACAGGTTTATCCCCGCTGGCGAGGGGa                 >  1:243857/1‑149 (MQ=255)
                                                                                                                               tCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGACGGTTTATCCCCGCTGGCGCGGGGAACACGAAAAGCTGGTGAGCGTTAAAGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGaa                >  2:393815/1‑149 (MQ=255)
                                                                                                                               tCGGTTTATCCCCGCTGACCCGGGGAACACTGGCGCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGaa                <  2:58882/149‑1 (MQ=255)
                                                                                                                                cGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAac               >  2:399743/1‑149 (MQ=255)
                                                                                                                                 cgTTTAACCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCTGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCTCAAACACAGGTTTATCCCCGCTGGCGCGGGGAaca              <  1:760053/148‑1 (MQ=255)
                                                                                                                                  gTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAacac             >  2:26725/1‑149 (MQ=255)
                                                                                                                                     tATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACAtt          <  1:1182636/149‑1 (MQ=255)
                                                                                                                                          ccGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGcc     >  2:1028512/1‑149 (MQ=255)
                                                                                                                                           cGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACAAATGACCAATTCTGTTTATCCCCGCTGGCGCGGGTCACACGAAATGCTGGGGAGCGTTAATGCCGAAAACACAGGTTTATCCACTCTGTCGCGGGGAAAACATTACGCCt    >  2:215940/1‑149 (MQ=255)
                                                                                                                                             cTGACGCGGGGAACACTGGCTCTGCATCATCAGAAACCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCttt  <  2:683190/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCTTT  >  NC_000913/2904068‑2904357

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