Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,904,212 G→A 88.9% 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,2120GA88.9% 187.6 / 16.3 83intergenic (+838/+535)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base A (36/35);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

GGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCTTTTT  >  NC_000913/2904069‑2904359
                                                                                                                                               |                                                                                                                                                   
ggCTCCGGTGTAGCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCgg                                                                                                                                                <  2:689325/149‑1 (MQ=255)
 gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTCTGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCggg                                                                                                                                               <  2:212764/149‑1 (MQ=255)
 gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCggg                                                                                                                                               <  1:482678/149‑1 (MQ=255)
 gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCggg                                                                                                                                               <  2:238931/149‑1 (MQ=255)
 gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCggg                                                                                                                                               <  1:282004/149‑1 (MQ=255)
 gCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAAGCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCggg                                                                                                                                               <  1:202656/149‑1 (MQ=255)
 gCTCCGGTTTATCCCCGCGGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTACTCGGTTTATCCCCGCTGACGCggg                                                                                                                                               <  1:29887/149‑1 (MQ=255)
  cTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCTCGGAATACATAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCgggg                                                                                                                                              >  1:252849/1‑149 (MQ=255)
  cTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCgggg                                                                                                                                              >  1:878284/1‑149 (MQ=255)
   tCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATATATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGTGGAACACGGGCGCACGGAATACAAAGAAGTGTATCAGCTCGGTTTATACCCGCTGACGCGGGGa                                                                                                                                             >  2:368571/1‑149 (MQ=255)
     cGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCATATCTATGGAGTATTGCTATTCGGTTTATCCACGCTGGCGCGGGGAACACGGGCGCTCGGAATACAAAGCGGGGTATCTGCTCGGGTTATCCCCGCTGCCGCGGGGAac                                                                                                                                           >  2:977391/1‑149 (MQ=255)
       gTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGCCGCGGGGAacac                                                                                                                                         >  1:865270/1‑149 (MQ=255)
       gTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAacac                                                                                                                                         <  2:431961/149‑1 (MQ=255)
        tttATCCCCGCTGGCGCGGGGAACACCTGTTTTCGGAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGGTTATCCCCGCTGACGCGGGGAACACt                                                                                                                                        <  1:919381/149‑1 (MQ=255)
        tttATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACt                                                                                                                                        <  1:526501/149‑1 (MQ=255)
          tATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTgg                                                                                                                                      <  1:805038/149‑1 (MQ=255)
          tATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCGCGGTTTATCCCCGCTGACGCGGGGAACACTgg                                                                                                                                      <  1:474479/149‑1 (MQ=255)
            tCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGcg                                                                                                                                    >  1:899951/1‑148 (MQ=255)
              cccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGctct                                                                                                                                  >  1:430527/1‑149 (MQ=255)
              cccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGctct                                                                                                                                  >  1:700683/1‑149 (MQ=255)
               ccgccgGCGCCGGGAACCACTGTATTCGCAAATCAATGGACTATTGCAATTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCGCACTGAATACAAAGGCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTg                                                                                                                                 <  2:223803/149‑1 (MQ=255)
               ccGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATAGCTAATCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTTGCTCTg                                                                                                                                 <  2:227454/149‑1 (MQ=255)
                cGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGc                                                                                                                                <  1:443388/149‑1 (MQ=255)
                cGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGc                                                                                                                                <  1:787994/149‑1 (MQ=255)
                  cTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATATGCTCGGTTTATCCCCGCTGACGCGGGGCACACTGGCTCTGCat                                                                                                                              >  1:326177/1‑148 (MQ=255)
                     gcgcGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCCCACGGAATACAAAGACGTGTATCTGATCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCAGCAAcag                                                                                                                           >  2:285186/1‑149 (MQ=255)
                      gggggtcaaCACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGTAACACGGGCGCACGGAATACAAAGCCGCGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAcagc                                                                                                                          <  2:933262/142‑1 (MQ=255)
                       gcgGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTTTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTTCAAcagca                                                                                                                         >  1:222343/1‑149 (MQ=255)
                       gcgGGGAACACCTGTTTTCGCAAATATATGGACTATTGCTATTCGGCTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAAcagca                                                                                                                         <  2:559649/149‑1 (MQ=255)
                          gggAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAgcagca                                                                                                                      >  1:647447/1‑149 (MQ=255)
                             aaCACCTTTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCTGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCAccc                                                                                                                   <  2:764714/149‑1 (MQ=255)
                               cacCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCAt                                                                                                                 <  1:446090/149‑1 (MQ=255)
                                acCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCCCGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACTCGGGGAACGCTGGCTCTGCAACAGCAGCACCCagg                                                                                                                >  2:256016/1‑147 (MQ=255)
                                  cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGAc                                                                                                              >  1:14672/1‑149 (MQ=255)
                                  cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCAGCAACAGCAGCACCCATGAc                                                                                                              >  2:694215/1‑149 (MQ=255)
                                    gTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCa                                                                                                            >  1:234656/1‑149 (MQ=255)
                                    gTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCa                                                                                                            <  1:321868/149‑1 (MQ=255)
                                          gCAAATCTCTGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCgg                                                                                                      <  1:240804/149‑1 (MQ=255)
                                          gCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCgg                                                                                                      <  2:267379/149‑1 (MQ=255)
                                           cAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACCGGCGCACGTAATACAAAGCCGTGAATCTGCTCGGTTTATCCCCGCTGACGCGGGCAACACTGGCTCTTCAACAGCAGCACCCATGACCACCTCGGt                                                                                                     >  1:783607/1‑149 (MQ=255)
                                               tCTATGGACTATTGCTATTCGGTTTATGCCCGCTGGCGCGGGGAACACGGGCGCACGGACTACAAAGCCGGGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTag                                                                                                 <  2:805605/149‑2 (MQ=255)
                                                 tATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATcc                                                                                               >  2:368190/1‑149 (MQ=255)
                                                      aCTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCt                                                                                          >  1:139631/1‑149 (MQ=255)
                                                          ttGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGAGAACGGAATACAAAGAAGTGTATCTGCTCGGTTTATACCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGAACCCATGACCAAATCCGTTTATACCCGCGGgcg                                                                                      >  2:992682/1‑149 (MQ=255)
                                                           tGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCCGCACCCATGACCACGTCGGTTTATCCCCGCTGgcgg                                                                                     >  2:992563/1‑148 (MQ=255)
                                                            gCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGcgcg                                                                                    <  2:234656/149‑1 (MQ=255)
                                                             cTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGACTACAAAGCCGTGTATCTGCTCCGTTTATCCCCGCTGACTCGGGGAACACTGGCTCTGCAACAGCAGCCCCCATGCCCACGTCGTTTTATCCCCGCGGGCGCgg                                                                                   >  1:579050/1‑149 (MQ=255)
                                                                ttCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATGCAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCACCAGCAGCACCCATGCCCACATCGGTTTATCCCCGCTGGCGCggggg                                                                                >  1:62681/1‑148 (MQ=255)
                                                                   ggTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCCCCCTAGCCCCCGTCGGTTTATCCCCGCTGGCggggggtcac                                                                             >  1:93340/1‑145 (MQ=255)
                                                                    gTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCAGCCATGACCACGTCGGTTTAACCCCGCTGTCGGGGGGAacac                                                                            >  2:611522/1‑149 (MQ=255)
                                                                        aTCCCCGCTGGCGCGGGGAACACGGGCGCACCGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCTCGGGGAACACGaaa                                                                        <  2:706850/149‑1 (MQ=255)
                                                                           cccGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGGTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGc                                                                     <  1:915749/149‑1 (MQ=255)
                                                                            ccGCTGGCGCGGGGAACCCGGGCGCACGGAATACAACGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCGGGCGCGGGGTGCACGAATTgat                                                                    >  2:51906/1‑147 (MQ=255)
                                                                               cTGGCGCGGGGAACACGGGCGCACGGAATACATAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTGTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGt                                                                 <  1:198893/149‑1 (MQ=255)
                                                                                    gcgGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCg                                                            <  2:574143/149‑1 (MQ=255)
                                                                                     cgGGGAACACGGGCGCACGGACAACAAAGCCGGGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCTCCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGt                                                           <  1:51595/149‑1 (MQ=255)
                                                                                       gggAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACCGCAGCACCCATGACCCCGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTa                                                         >  2:488662/1‑149 (MQ=255)
                                                                                             acGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCg                                                   <  2:14672/149‑1 (MQ=255)
                                                                                               gggCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGTCTCTGAAACAGCAGCACCCATGAACACGTCTGTTTATCCCCGCACGCGCGGGGCACAAGAAATGCTGGTTCGCGTTAATGCCGCa                                                 >  2:744841/1‑149 (MQ=255)
                                                                                                 gcgcACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCaaa                                               <  2:878284/149‑1 (MQ=255)
                                                                                                 gcgcACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGCCCCCGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCaaa                                               >  1:974289/1‑149 (MQ=255)
                                                                                                       ggAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAGCACAgg                                         >  2:52916/1‑149 (MQ=255)
                                                                                                       ggAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCCTGACCACGTCGGTTTATCCCCGCGGGCGCGGGGCACACGAACTGCGGGTGAGCGTTAATGCCGCAAAcacacg                                         >  1:117139/1‑147 (MQ=255)
                                                                                                       ggAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAgg                                         <  1:722354/149‑1 (MQ=255)
                                                                                                         aaTACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGtt                                       <  1:1010429/149‑1 (MQ=255)
                                                                                                         aaTACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGtt                                       >  2:133671/1‑149 (MQ=255)
                                                                                                          aTACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCCCCAGCCGCACCCATGCCCCCGGCGGTTTATCCCCGCGGGCGCGGGGAGCACGAAATGCTGGTGAGCGTTAATGCCGCAAACCCAGGttt                                      >  2:719422/1‑149 (MQ=255)
                                                                                                             cAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCAGCAAACACAGGTTTATc                                   >  2:415545/1‑149 (MQ=255)
                                                                                                              aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCCCACATGACCAAGTCGGTTTATCACCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGATAATGCCGGAAACACAGGTTTATcc                                  >  1:773589/1‑149 (MQ=255)
                                                                                                              aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTACTGCCGCAAACACAGGTTTATcc                                  >  1:755435/1‑149 (MQ=255)
                                                                                                              aaaGCCGTGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATcc                                  >  1:276046/1‑149 (MQ=255)
                                                                                                               aaGCCGGGTATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATccc                                 <  2:600322/149‑1 (MQ=255)
                                                                                                                     tgtATCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGTAGCACCCATGACCAAGTCGGTTTATCCCCGCTGGCGCGGGGAACAAGAAATGCTGGTGAGCGTAAATGCCGCAATCACAGGTATATCCCCGCTgg                           >  2:756331/1‑149 (MQ=255)
                                                                                                                         tCTGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGcgcg                       >  1:363982/1‑149 (MQ=255)
                                                                                                                           tGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGATGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCggg                     <  2:647447/149‑1 (MQ=255)
                                                                                                                           tGCTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGTCCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCggg                     <  2:429439/149‑1 (MQ=255)
                                                                                                                             cTCGGTTTATCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCATGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGa                   >  2:598834/1‑149 (MQ=255)
                                                                                                                                ggTTTATCCCCGCTGACGCGGGGAACACTGGATCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAaca                <  2:689947/149‑1 (MQ=255)
                                                                                                                                     atcacagCTGACGCGGGGAAAACTGGCTCTGCAACTGCAGCACCTATGTCTAAGTAGGTTTTTCCGCGCTGGCGCTGGTAACACGAATTGCTGGTGATTGTTAATGACGCATACACTGGTTTATCCCCGCTGGCGCGGGGAACACATTa           <  2:700683/143‑1 (MQ=255)
                                                                                                                                     aTCCCCGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTa           <  1:694215/149‑1 (MQ=255)
                                                                                                                                       ccccGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCCAACACAGGTTTATCCCCGCTGGCGCGGGGAACACAGTACg         >  2:120631/1‑149 (MQ=255)
                                                                                                                                        cccGCTGCCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGCCCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGCGGAGCGTTAATGCCGCAACCACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGc        <  1:190843/149‑1 (MQ=255)
                                                                                                                                         ccGCTGACGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGGTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGcc       >  2:268363/1‑149 (MQ=255)
                                                                                                                                              ggCGCGGGGAACACTGTCTCTGCAACAGCAGCACCCAAGAGCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCttttt  <  2:1001810/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAACAGCAGCACCCATGACCACGTCGGTTTATCCCCGCTGGCGCGGGGAACACGAAATGCTGGTGAGCGTTAATGCCGCAAACACAGGTTTATCCCCGCTGGCGCGGGGAACACATTACGCCTTTTT  >  NC_000913/2904069‑2904359

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: