Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,347,1690GA43.3% 13.0 / 97.0 67intergenic (+19/‑215)nrdA/nrdBribonucleoside‑diphosphate reductase 1, alpha subunit/ribonucleoside‑diphosphate reductase 1, beta subunit, ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (19/19);  new base A (9/20);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-02

GCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGC  >  NC_000913/2347025‑2347315
                                                                                                                                                |                                                                                                                                                   
aacgccatcgactttttatgcctttgcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCgtg                                                                                                                                                <  2:1255399‑M1/16‑2 (MQ=255)
 ccTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTaa                                                                                                                                               >  2:3403054/1‑150 (MQ=255)
 acgccatcgactttttatgcctttgcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAgcgtgaa                                                                                                                                              >  1:2477932‑M1/134‑148 (MQ=255)
  tgccatcgactttttatgcctttgcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAAc                                                                                                                                             <  1:3506154‑M1/19‑1 (MQ=255)
  cgccatcgactttttatgcctttgcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAAc                                                                                                                                             <  1:3283678‑M1/19‑1 (MQ=255)
   gccatcgactttttatgcctgtgcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAAc                                                                                                                                             >  1:2599033‑M1/132‑150 (MQ=255)
    ccatcgactttttatgcctttgcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGc                                                                                                                                           <  2:1951649‑M1/21‑1 (MQ=255)
    aCAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAAc                                                                                                                                             <  2:2334006/149‑1 (MQ=255)
     catcgactttttatgcctttgcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGcc                                                                                                                                          <  2:3789037‑M1/22‑1 (MQ=255)
          tCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                     >  2:3813538/1‑151 (MQ=255)
           ctttttatgcctttgcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATcc                                                                                                                                    <  1:3223818‑M1/28‑1 (MQ=255)
              tgTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                 >  1:1430324/2‑151 (MQ=255)
                     ctttgcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGc                                                                                                                          <  2:1481823‑M1/38‑1 (MQ=255)
                     ctttgcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGc                                                                                                                          <  2:409983‑M1/38‑1 (MQ=255)
                       acTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTc                                                                                                                        <  1:3858874/151‑1 (MQ=255)
                         gcggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCgg                                                                                                                      <  2:3033112‑M1/42‑1 (MQ=255)
                          gTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCgg                                                                                                                      <  2:2112065/150‑1 (MQ=255)
                          cggcatcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCgg                                                                                                                      <  2:2038646‑M1/42‑1 (MQ=255)
                            attatCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGttt                                                                                                                  >  1:480406/1‑151 (MQ=255)
                            attatCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGtt                                                                                                                   <  2:1153121/150‑1 (MQ=255)
                             catcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGgctcggattt                                                                                                                  <  2:169878‑M1/46‑5 (MQ=255)
                             catcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGgctcggattt                                                                                                                  >  2:2998939‑M1/106‑147 (MQ=255)
                              atcgggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACggctcggatttg                                                                                                                 <  2:50653‑M1/47‑6 (MQ=255)
                                 gggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTAcggctcggatttgt                                                                                                                <  2:2651724‑M1/48‑7 (MQ=255)
                                 gggcaatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGTTCGGATTTGTAg                                                                                                              >  2:1252793‑M1/102‑151 (MQ=255)
                                 gagtacacgagTGACGGGGCTGGAGAGGCAGGGGCCGTTCTTGTGACGTCACTCGGGGGCGACGGCTGCGAAAAAAAAGCATGTAAAATCTGATATTGAGATGACGGATGCGGCGTAAACGCCTTATACAGCCTACGGCTCGGTTTGtg                                                                                                               >  1:2314417/11‑148 (MQ=255)
                                 cAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAgg                                                                                                             <  2:2282063/151‑1 (MQ=255)
                                 cAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAgg                                                                                                             >  1:729715/1‑151 (MQ=255)
                                    aaCACCCGTGACGGCGCTGAAGACGCACAAGACGCTCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGc                                                                                                            >  1:478543/1‑149 (MQ=255)
                                     acacCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTg                                                                                                         >  2:1825558/1‑151 (MQ=255)
                                     aatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCt                                                                                                          >  2:2375416‑M1/98‑151 (MQ=255)
                                      atgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTg                                                                                                         <  2:3308986‑M1/55‑1 (MQ=255)
                                         ccGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATaa                                                                                                     <  1:89463/151‑1 (MQ=255)
                                           tgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATa                                                                                                      <  2:2826025‑M1/58‑1 (MQ=255)
                                           gTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGa                                                                                                   <  1:1848841/151‑1 (MQ=255)
                                             gACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGAc                                                                                                  <  1:3137822/150‑1 (MQ=255)
                                              aCGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGAcgc                                                                                                <  1:115500/151‑1 (MQ=255)
                                               ggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgc                                                                                                <  2:3184358‑M1/64‑1 (MQ=255)
                                               ggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgc                                                                                                >  2:3957953‑M1/88‑151 (MQ=255)
                                                 gcgcTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGcc                                                                                             <  2:3052398/151‑1 (MQ=255)
                                                     ggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAg                                                                                           <  1:1994790‑M1/69‑1 (MQ=255)
                                                     ggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGc                                                                                          <  1:831601‑M1/70‑1 (MQ=255)
                                                      gcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCg                                                                                         >  1:2358907‑M1/81‑151 (MQ=255)
                                                       cgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGt                                                                                        >  1:1842636‑M1/80‑151 (MQ=255)
                                                          aaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGc                                                                                     <  2:3946512‑M1/75‑1 (MQ=255)
                                                               caAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAgg                                                                               <  1:3038062/151‑1 (MQ=255)
                                                                     ccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAggcac                                                                            <  2:3426268‑M1/84‑3 (MQ=255)
                                                                         tGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGc                                                                     >  2:1438957/1‑151 (MQ=255)
                                                                          ggTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGt                                                                       <  1:3263484/148‑1 (MQ=255)
                                                                            acggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTcgcatcaggcacaggatgc                                                                     >  2:3605588‑M1/59‑140 (MQ=255)
                                                                               cGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCCCCAGCGTCGCATCAGGCTCCGGGTGCCGGATg                                                               <  1:145278/151‑1 (MQ=255)
                                                                                         ggACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAAc                                                     <  2:3079205/151‑1 (MQ=255)
                                                                                             gATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGACAAGACGCGCCAGCGCCGCATCAGTCTCCGGGAGCCGGATGCAGCGTGAACGcc                                                  >  1:3727611/1‑150 (MQ=255)
                                                                                                ggCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGAc                                                                                                  >  2:3304471/1‑99 (MQ=255)
                                                                                                ggCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGAc                                                                                                  <  1:3304264/99‑1 (MQ=255)
                                                                                                  cTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATcc                                            >  1:2429954/1‑151 (MQ=255)
                                                                                                        aaGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGcc                                        >  1:2066654/1‑149 (MQ=255)
                                                                                                          gcggcgCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACg                                    <  2:106252/151‑1 (MQ=255)
                                                                                                                       gATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCCGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGGTTTGt                        >  1:311077/1‑150 (MQ=255)
                                                                                                                            gATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAgg                     <  1:3732929/148‑1 (MQ=255)
                                                                                                                             atatTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTg                 >  2:1633318/1‑151 (MQ=255)
                                                                                                                                ttGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGAt               >  2:2460914/1‑150 (MQ=255)
                                                                                                                                 tGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATa              >  2:802347/1‑150 (MQ=255)
                                                                                                                                   agaTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGAt               >  2:1563972/1‑147 (MQ=255)
                                                                                                                                       gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgcg       <  1:248895/151‑1 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACgcgc      <  2:2122763/151‑1 (MQ=255)
                                                                                                                                            aTGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGc  >  2:2639664/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGGTTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGC  >  NC_000913/2347025‑2347315

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