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

GCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGC  >  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)
                            attatCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑ttt                                                                                                                  >  1:480406/1‑151 (MQ=255)
                            attatCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑tt                                                                                                                   <  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)
                                 gagtacacgagTGACGGGGCTGGAGAGGCAGGGGCCGTTCTTGTGACGTCACTCGGGGGCGACGGCTGCGAAAAAAAAGCATGTAAAATCTGATATTGAGATGACGGATGCGGCGTAAACGCCTTATACAGCCTACGGCTCGG‑TTTGtg                                                                                                               >  1:2314417/11‑148 (MQ=255)
                                 cAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAgg                                                                                                             <  2:2282063/151‑1 (MQ=255)
                                 cAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAgg                                                                                                             >  1:729715/1‑151 (MQ=255)
                                    aaCACCCGTGACGGCGCTGAAGACGCACAAGACGCTCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGc                                                                                                            >  1:478543/1‑149 (MQ=255)
                                     acacCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTg                                                                                                         >  2:1825558/1‑151 (MQ=255)
                                     aatgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCt                                                                                                          >  2:2375416‑M1/98‑151 (MQ=255)
                                      atgcgtgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTg                                                                                                         <  2:3308986‑M1/55‑1 (MQ=255)
                                         ccGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATaa                                                                                                     <  1:89463/151‑1 (MQ=255)
                                           tgtcggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATa                                                                                                      <  2:2826025‑M1/58‑1 (MQ=255)
                                           gTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGa                                                                                                   <  1:1848841/151‑1 (MQ=255)
                                             gACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGAc                                                                                                  <  1:3137822/150‑1 (MQ=255)
                                              aCGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGAcgc                                                                                                <  1:115500/151‑1 (MQ=255)
                                               ggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgc                                                                                                <  2:3184358‑M1/64‑1 (MQ=255)
                                               ggatgcggcgtaaacgccttatccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgc                                                                                                >  2:3957953‑M1/88‑151 (MQ=255)
                                                 gcgcTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGcc                                                                                             <  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)
                                                               caAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAgg                                                                               <  1:3038062/151‑1 (MQ=255)
                                                                     ccgacctacggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAggcac                                                                            <  2:3426268‑M1/84‑3 (MQ=255)
                                                                         tGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGc                                                                     >  2:1438957/1‑151 (MQ=255)
                                                                          ggTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGt                                                                       <  1:3263484/148‑1 (MQ=255)
                                                                            acggttctacccctgcgtaggcctgataagacgcgccagcgtcgcatcaggctccgggTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTcgcatcaggcacaggatgc                                                                     >  2:3605588‑M1/59‑140 (MQ=255)
                                                                               cGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCCCCAGCGTCGCATCAGGCTCCGGGTGCCGGATg                                                               <  1:145278/151‑1 (MQ=255)
                                                                                         ggACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAAc                                                     <  2:3079205/151‑1 (MQ=255)
                                                                                             gATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGACAAGACGCGCCAGCGCCGCATCAGTCTCCGGGAGCCGGATGCAGCGTGAACGcc                                                  >  1:3727611/1‑150 (MQ=255)
                                                                                                ggCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGAc                                                                                                  >  2:3304471/1‑99 (MQ=255)
                                                                                                ggCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGAc                                                                                                  <  1:3304264/99‑1 (MQ=255)
                                                                                                  cTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATcc                                            >  1:2429954/1‑151 (MQ=255)
                                                                                                        aaGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGcc                                        >  1:2066654/1‑149 (MQ=255)
                                                                                                          gcggcgCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACg                                    <  2:106252/151‑1 (MQ=255)
                                                                                                                       gATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCCGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGGTTTGt                        >  1:311077/1‑150 (MQ=255)
                                                                                                                            gATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAgg                     <  1:3732929/148‑1 (MQ=255)
                                                                                                                             atatTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTg                 >  2:1633318/1‑151 (MQ=255)
                                                                                                                                ttGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGAt               >  2:2460914/1‑150 (MQ=255)
                                                                                                                                 tGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATa              >  2:802347/1‑150 (MQ=255)
                                                                                                                                   agaTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGAt               >  2:1563972/1‑147 (MQ=255)
                                                                                                                                       gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGAcgcg       <  1:248895/151‑1 (MQ=255)
                                                                                                                                        ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACgcgc      <  2:2122763/151‑1 (MQ=255)
                                                                                                                                            aTGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGc  >  2:2639664/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCCTACAAATTCGGGGTCAAAACACTGTATTATCAGAACACCCGTGACGGCGCTGAAGACGCACAAGACGATCTGGTGCCGTCAATCCAGGACGATGGCTGCGAAAGCGGCGCATGTAAGATCTGATATTGAGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCTCCGGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGC  >  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: ‑