Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,0010AC38.1% 68.2 / 10.1 43H26P (CAC→CCC) ftnBputative ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (20/6);  new base C (0/16);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 4.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATTTAAGCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTT  >  NC_000913/1986752‑1987235
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           
gatTTAAGCATTTTCTTATACCCGTTCAGACGTTATTGTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCAc                                                                                                                                                                                                                                            <  2:91866/251‑1 (MQ=255)
                   aCCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCt                                                                                                                                                                                                                         <  1:102601/250‑1 (MQ=255)
                               gTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTg                                                                                                                                                                                                             >  1:106835/1‑250 (MQ=255)
                                     cTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGc                                                                                                                                                                                                       >  1:178984/1‑250 (MQ=255)
                                              gATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCGCCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACtt                                                                                                                                                                                             <  1:385783/251‑1 (MQ=255)
                                                 catcGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTcc                                                                                                                                                                                          >  2:352514/1‑251 (MQ=255)
                                                         gAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGcc                                                                                                                                                                                  >  1:181324/1‑251 (MQ=255)
                                                                       aTCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCCACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCCTCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCCGGCACAGAGTa                                                                                                                                                                     >  2:88204/1‑250 (MQ=255)
                                                                                cGGACGGGCAAAAAGATATTATTCTCAAAACAATAATTTTGGTGTTTTGTTTTATCTTGCCCCTCCCCACTTTTGTATATAAGAATATTAGGTATGGCACCCGCTGAAATCCTTCTCAAACCC‑AACTCTAAAAAGAACCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTACC‑TGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGAccc                                                                                                                                                            <  2:298757/250‑1 (MQ=255)
                                                                                cGGACCGGCAGAAAGGATTTATTCTGCAAACAGTAATTTTGGTGTTTTTATTTATCTTGCCCCTCTCCCCTTCTGAATATAAGGATATTGGGTATGGCACCCGCTGGAATGCTTCTCAAACTC‑ACCTCTCAAAGGAACCGCGATTTTTACGCATCCAATCTCTCCCTTCCCCGGAGTAAC‑GGGTGTTCTGAACAGAGTCTGAACGGCCCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCa                                                                                                                                                           <  2:29159/251‑1 (MQ=255)
                                                                                 ggACCGGAAAAAAGGAATTATTCTGCAAAAAGTAATTATGGTTTTTTTATTTTTCTTCCCCCTCTCCACTTCTGGATAAAAGGATTTTAGGTATGGCACCCGCGGAAATGCTTCTCAAACTC‑AACTCTCAAATGACCCGCGTTTTTAACGCACCCACTCTCTCCCTTCCCCTGAGTAAC‑TGGTGTTCTGAACAGAGTATGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCaa                                                                                                                                                          <  2:106835/251‑1 (MQ=255)
                                                                                      ggaagaaagaaTTTATTCTGCAAACAGTATTTATGGTGTTTTGATTTATCTGGCACCTCTCCCCTTTGGGAAATAAGAATATTAGTTAGGGCACCCGCTGGAATGCTTCTCAAACCC‑AACTCTCAAAGGACCCGCGATTTTTACGCATCCAATCTCTCCCTCCCCCTGAGTAAC‑TGGTGTTCGGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAatgat                                                                                                                                                     <  2:178984/248‑1 (MQ=255)
                                                                                        cAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgc                                                                                                                                                   >  2:24552/1‑251 (MQ=255)
                                                                                             aGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATg                                                                                                                                              >  1:541749/1‑251 (MQ=255)
                                                                                              ggAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGt                                                                                                                                             >  1:100519/1‑251 (MQ=255)
                                                                                                   tattCTGCAAACATTAATTATGGTTTTTTGTTTTACCTTCCCCCTCCCCCCTTCTGAATATAAGGATATTAGGTATGGCAACCGCTGAAATGCTTCTCAAACTC‑ACCTCTCAAATGAACCGCGATTTTTACGCATCCAATCTCTCCCTCCCCCTGAGTACC‑TGGTTTTCGGAACAGATTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTaa                                                                                                                                         <  1:523114/250‑1 (MQ=255)
                                                                                                       cTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCAcc                                                                                                                                                                                                    <  1:540629/187‑1 (MQ=255)
                                                                                                       cTGCAAACAGAAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCAcc                                                                                                                                                                                                    >  2:540629/1‑187 (MQ=255)
                                                                                                             aCAGTAATTATGGTTTTTTGATTTATCTTGCCCCTCTCCCCTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑ACCTCTAAAATGACCCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAaga                                                                                                                              <  1:24552/251‑1 (MQ=255)
                                                                                                                gTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTCACTTTATGAAGAgtg                                                                                                                           >  2:495772/1‑251 (MQ=255)
                                                                                                                               tGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTCCCCCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTAccccc                                                                                                             <  1:159501/250‑1 (MQ=255)
                                                                                                                                   ttATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCACCCGCTGGAATGCTTCTCAAACTC‑AACTCTAAAATGAACCGCGAGTTTTACGCATCCAATCTCTCCCTTCACCTGAGTAAC‑TGGTGTTCTAAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCg                                                                                                         <  1:176067/250‑1 (MQ=255)
                                                                                                                                                         tGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGttaacatc                                                                                  >  1:433264/1‑245 (MQ=255)
                                                                                                                                                                agtgATATTAGTTATGGAACCCGCTGGAATGCTTC‑CAAACCCCACCCCTAAAAGGACCCGTGGTTTTTACACATCCATTCTCTCCCTCCCCCGAAGTACC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAAc                                                                           <  2:412942/248‑1 (MQ=255)
                                                                                                                                                                  ggATATTAGGTATGGCAACCGCTGAAATGCTTCTAAAACCC‑ACCTCTAAAATGACCGGCGATTTTTCCGCCTCCAATCCCTCCTTTCCCCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTAAACGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTg                                                                         <  2:352200/251‑1 (MQ=255)
                                                                                                                                                                                  caCCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAAc                                                         >  2:189730/2‑251 (MQ=255)
                                                                                                                                                                                     cGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGt                                                      >  1:112550/1‑251 (MQ=255)
                                                                                                                                                                                        tGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGATTTTTACGCATCCAATCTCTCCCTCCCCCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTAAAGGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgc                                                                                                                                                   <  2:34496/155‑1 (MQ=255)
                                                                                                                                                                                        tGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgc                                                                                                                                                   >  1:34496/1‑155 (MQ=255)
                                                                                                                                                                                             tGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCaaaa                                               >  2:390135/1‑250 (MQ=255)
                                                                                                                                                                                              gCTTCTCAACCAC‑AACTCCCAAATGACCCGCGAGTTTTACGCATCCAATCTCTACCTCCCCCTGATTACC‑TGTTGTTCTGAACAGAGTCTGAACGGCTCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCaaaaa                                              <  1:405505/250‑1 (MQ=255)
                                                                                                                                                                                                      aactc‑aactcTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGga                                      >  2:103431/1‑250 (MQ=255)
                                                                                                                                                                                                          c‑aactcTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAAt                                 >  2:503629/1‑251 (MQ=255)
                                                                                                                                                                                                             gcactcAAATGAACCGCGATTTTTACGAACCAAATCTCCACCTTCCCCGGCGAAAT‑TGGTGTTCTGACCAGAGGCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATAc                               <  1:390135/248‑1 (MQ=255)
                                                                                                                                                                                                             actcgaaaaTGAACCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTAAC‑GGGTGTTCTAAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATAc                               <  1:88204/245‑1 (MQ=255)
                                                                                                                                                                                                             actcTAAAATAAACCGCGAGTTTAACGCATCCACTCCCTACCTCCACCTGAGTAAC‑GGGGTTTCGGAACAGAGTGTAAAGGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATAc                               <  2:181324/251‑1 (MQ=255)
                                                                                                                                                                                                                cccaaaaaacccGCTAGTTTGACCCAACCAATCTCTCCTTTCCCCTGATTATC‑TGGGGTTCTGAACAGAGTCAGACCGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAg                            <  2:495173/241‑1 (MQ=255)
                                                                                                                                                                                                                   aaaTGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCa                          >  1:336591/1‑250 (MQ=255)
                                                                                                                                                                                                                      tGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGt                      >  1:395797/1‑251 (MQ=255)
                                                                                                                                                                                                                              gagttttcccctcccAATCTCTACCTTCACCTGAGTACC‑TGTTTTTCTGAACAGAGTCTGACCGGCCCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTAc              <  1:44434/238‑1 (MQ=255)
                                                                                                                                                                                                                                    tACGCACCCATTCTCTCCCTTCCCCTGA‑TAACCTGGTGTTCTAAACAGAGTCTGACGGGCCCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTgg         <  2:395602/250‑1 (MQ=255)
                                                                                                                                                                                                                                     aCGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGc        >  1:517951/1‑250 (MQ=255)
                                                                                                                                                                                                                                          cccAATCTCTCCCTCCCCCTGAGAAAC‑TGGTGTTCTGAACGGAGTCTGAACGCCCCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGtt  <  1:39348/250‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           
GATTTAAGCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTC‑AACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC‑TGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTT  >  NC_000913/1986752‑1987235

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