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.

GATTTAAGCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTT  >  NC_000913/1986752‑1987235
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           
gatTTAAGCATTTTCTTATACCCGTTCAGACGTTATTGTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCAc                                                                                                                                                                                                                                            <  2:91866/251‑1 (MQ=255)
                   aCCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCt                                                                                                                                                                                                                         <  1:102601/250‑1 (MQ=255)
                               gTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTg                                                                                                                                                                                                             >  1:106835/1‑250 (MQ=255)
                                     cTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGc                                                                                                                                                                                                       >  1:178984/1‑250 (MQ=255)
                                              gATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCGCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACtt                                                                                                                                                                                             <  1:385783/251‑1 (MQ=255)
                                                 catcGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTcc                                                                                                                                                                                          >  2:352514/1‑251 (MQ=255)
                                                         gAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGcc                                                                                                                                                                                  >  1:181324/1‑251 (MQ=255)
                                                                       aTCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCCACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCCTCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCCGGCACAGAGTa                                                                                                                                                                     >  2:88204/1‑250 (MQ=255)
                                                                                cGGACGGGCAAAAAGATATTATTCTCAAAACAATAATTTTGGTGTTTTGTTTTATCTTGCCCCTCCCCACTTTTGTATATAAGAATATTAGGTATGGCACCCGCTGAAATCCTTCTCAAACCCAACTCTAAAAAGAACCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTACCTGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGAccc                                                                                                                                                            <  2:298757/250‑1 (MQ=255)
                                                                                cGGACCGGCAGAAAGGATTTATTCTGCAAACAGTAATTTTGGTGTTTTTATTTATCTTGCCCCTCTCCCCTTCTGAATATAAGGATATTGGGTATGGCACCCGCTGGAATGCTTCTCAAACTCACCTCTCAAAGGAACCGCGATTTTTACGCATCCAATCTCTCCCTTCCCCGGAGTAACGGGTGTTCTGAACAGAGTCTGAACGGCCCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCa                                                                                                                                                           <  2:29159/251‑1 (MQ=255)
                                                                                 ggACCGGAAAAAAGGAATTATTCTGCAAAAAGTAATTATGGTTTTTTTATTTTTCTTCCCCCTCTCCACTTCTGGATAAAAGGATTTTAGGTATGGCACCCGCGGAAATGCTTCTCAAACTCAACTCTCAAATGACCCGCGTTTTTAACGCACCCACTCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTATGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCaa                                                                                                                                                          <  2:106835/251‑1 (MQ=255)
                                                                                      ggaagaaagaaTTTATTCTGCAAACAGTATTTATGGTGTTTTGATTTATCTGGCACCTCTCCCCTTTGGGAAATAAGAATATTAGTTAGGGCACCCGCTGGAATGCTTCTCAAACCCAACTCTCAAAGGACCCGCGATTTTTACGCATCCAATCTCTCCCTCCCCCTGAGTAACTGGTGTTCGGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAatgat                                                                                                                                                     <  2:178984/248‑1 (MQ=255)
                                                                                        cAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgc                                                                                                                                                   >  2:24552/1‑251 (MQ=255)
                                                                                             aGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATg                                                                                                                                              >  1:541749/1‑251 (MQ=255)
                                                                                              ggAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGt                                                                                                                                             >  1:100519/1‑251 (MQ=255)
                                                                                                   tattCTGCAAACATTAATTATGGTTTTTTGTTTTACCTTCCCCCTCCCCCCTTCTGAATATAAGGATATTAGGTATGGCAACCGCTGAAATGCTTCTCAAACTCACCTCTCAAATGAACCGCGATTTTTACGCATCCAATCTCTCCCTCCCCCTGAGTACCTGGTTTTCGGAACAGATTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTaa                                                                                                                                         <  1:523114/250‑1 (MQ=255)
                                                                                                       cTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCAcc                                                                                                                                                                                                    <  1:540629/187‑1 (MQ=255)
                                                                                                       cTGCAAACAGAAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCAcc                                                                                                                                                                                                    >  2:540629/1‑187 (MQ=255)
                                                                                                             aCAGTAATTATGGTTTTTTGATTTATCTTGCCCCTCTCCCCTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCACCTCTAAAATGACCCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAaga                                                                                                                              <  1:24552/251‑1 (MQ=255)
                                                                                                                gTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTCACTTTATGAAGAgtg                                                                                                                           >  2:495772/1‑251 (MQ=255)
                                                                                                                               tGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTCCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTGACCGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTAccccc                                                                                                             <  1:159501/250‑1 (MQ=255)
                                                                                                                                   ttATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCACCCGCTGGAATGCTTCTCAAACTCAACTCTAAAATGAACCGCGAGTTTTACGCATCCAATCTCTCCCTTCACCTGAGTAACTGGTGTTCTAAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCg                                                                                                         <  1:176067/250‑1 (MQ=255)
                                                                                                                                                         tGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGttaacatc                                                                                  >  1:433264/1‑245 (MQ=255)
                                                                                                                                                                agtgATATTAGTTATGGAACCCGCTGGAATGCTTCCAAACCCCACCCCTAAAAGGACCCGTGGTTTTTACACATCCATTCTCTCCCTCCCCCGAAGTACCTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAAc                                                                           <  2:412942/248‑1 (MQ=255)
                                                                                                                                                                  ggATATTAGGTATGGCAACCGCTGAAATGCTTCTAAAACCCACCTCTAAAATGACCGGCGATTTTTCCGCCTCCAATCCCTCCTTTCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTAAACGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTg                                                                         <  2:352200/251‑1 (MQ=255)
                                                                                                                                                                                  caCCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAAc                                                         >  2:189730/2‑251 (MQ=255)
                                                                                                                                                                                     cGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGt                                                      >  1:112550/1‑251 (MQ=255)
                                                                                                                                                                                        tGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGATTTTTACGCATCCAATCTCTCCCTCCCCCTGAGTAACTGGTGTTCTGAACAGAGTCTAAAGGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgc                                                                                                                                                   <  2:34496/155‑1 (MQ=255)
                                                                                                                                                                                        tGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATgc                                                                                                                                                   >  1:34496/1‑155 (MQ=255)
                                                                                                                                                                                             tGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCaaaa                                               >  2:390135/1‑250 (MQ=255)
                                                                                                                                                                                              gCTTCTCAACCACAACTCCCAAATGACCCGCGAGTTTTACGCATCCAATCTCTACCTCCCCCTGATTACCTGTTGTTCTGAACAGAGTCTGAACGGCTCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCaaaaa                                              <  1:405505/250‑1 (MQ=255)
                                                                                                                                                                                                      aactcaactcTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGga                                      >  2:103431/1‑250 (MQ=255)
                                                                                                                                                                                                          caactcTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAAt                                 >  2:503629/1‑251 (MQ=255)
                                                                                                                                                                                                             gcactcAAATGAACCGCGATTTTTACGAACCAAATCTCCACCTTCCCCGGCGAAATTGGTGTTCTGACCAGAGGCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATAc                               <  1:390135/248‑1 (MQ=255)
                                                                                                                                                                                                             actcgaaaaTGAACCGCGAGTTTTACGCATCCAATCTCTCCCTTCCCCTGAGTAACGGGTGTTCTAAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATAc                               <  1:88204/245‑1 (MQ=255)
                                                                                                                                                                                                             actcTAAAATAAACCGCGAGTTTAACGCATCCACTCCCTACCTCCACCTGAGTAACGGGGTTTCGGAACAGAGTGTAAAGGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATAc                               <  2:181324/251‑1 (MQ=255)
                                                                                                                                                                                                                cccaaaaaacccGCTAGTTTGACCCAACCAATCTCTCCTTTCCCCTGATTATCTGGGGTTCTGAACAGAGTCAGACCGCCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAg                            <  2:495173/241‑1 (MQ=255)
                                                                                                                                                                                                                   aaaTGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCa                          >  1:336591/1‑250 (MQ=255)
                                                                                                                                                                                                                      tGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGt                      >  1:395797/1‑251 (MQ=255)
                                                                                                                                                                                                                              gagttttcccctcccAATCTCTACCTTCACCTGAGTACCTGTTTTTCTGAACAGAGTCTGACCGGCCCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTAc              <  1:44434/238‑1 (MQ=255)
                                                                                                                                                                                                                                    tACGCACCCATTCTCTCCCTTCCCCTGATAACCTGGTGTTCTAAACAGAGTCTGACGGGCCCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTgg         <  2:395602/250‑1 (MQ=255)
                                                                                                                                                                                                                                     aCGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGc        >  1:517951/1‑250 (MQ=255)
                                                                                                                                                                                                                                          cccAATCTCTCCCTCCCCCTGAGAAACTGGTGTTCTGAACGGAGTCTGAACGCCCCCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGtt  <  1:39348/250‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           
GATTTAAGCATTTTCTTATACCCGTTCAGACGTTATTCTTATTTCAGATCATCGTCAGAATTGACTCCACGATCACATTTCGGACCGGCAGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGCACCGCCACTTTCCTTCGCGCCCAGGCACAGAGTAATGTGACCCAAATGATGCGCATGTTTAACTTTATGAAGAGTGTCGGCGCTACCCCCATCGTTAAAGCCATTGATGTTCCCGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTT  >  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: