Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,751,5370AC43.8% 48.4 / 12.4 48V111G (GTT→GGT) yiaTouter membrane protein YiaT
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base C (0/21);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 6.44e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATCCTGCGCGGAATAACTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGCGATTCATTATAAGC  >  NC_000913/3751295‑3751778
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
cATCCTGCGCGGAATAACTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGcc                                                                                                                                                                                                                                           <  2:453723/251‑1 (MQ=255)
              aaaCTTGCCAGACCGCTACGACGGGACTCCCTCTCTGAAATCCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATCCCGATAGCTCCCCCCCCCAGCCGTTGTTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATc                                                                                                                                                                                                                             <  1:328866/250‑1 (MQ=255)
              aaaCTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTAAAATTCTCGTCATAATAGAGAACGCCCAGCGCGGGTGTCAGCGACAGACGACCTATCGGCATTTTGTGGAATACCGATAGCTCCCCCCCCCAGCCGTTGCTGTTATCCAGACCGTCCGCAGCTGCAGAGGCTTTCACGTTTC‑‑CCAAACGCTCGTGGTAATACCAGTAACCCCCCGCCATAGCCGTACTATc                                                                                                                                                                                                                             <  1:114842/250‑1 (MQ=255)
                 cTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGCTACCAGGCAACCCCCGCCATAGCCGTACTATCACg                                                                                                                                                                                                                          >  1:251427/1‑251 (MQ=255)
                    gCCAGACCGCTCCGAGGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCGGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGAAATACCGATAGCTCCCCCCCCCAGCCTTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCt                                                                                                                                                                                                                        <  1:509097/250‑1 (MQ=255)
                                          cactcTCTGAAAGGCCATAGTAATAGCCCCTGAATTTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGCCCTATCTGCATTTTGTGGAATCCCGAAAGCCCCCCCCCCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATa                                                                                                                                                                                                  <  1:206305/250‑1 (MQ=255)
                                                    aaTGCCATATTAATAGCCACTGAATCTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGCCGCCCTATCTTCATTTGTGGGAATCCCGATAGCCCCCCCACCCAGCCGTTGCTTTTACCCAGAACGTCCGCAGCTGCAGAGGCTCTCACGTTTC‑‑CCCAACGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTAt                                                                                                                                                                                       <  2:481241/251‑1 (MQ=255)
                                                     aTGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGCCGGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCCCCCAGCCGTTGCTTTTATCCAGAACGTCCGCAGCTGCAGAGGTTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATc                                                                                                                                                                                      <  1:507478/251‑1 (MQ=255)
                                                         cATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                                                                                                                   >  2:112500/1‑250 (MQ=255)
                                                         cATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                                                                                                                   >  2:192100/1‑250 (MQ=255)
                                                         cATAGTAATAGCAACTGAAATTCTCGTCAAAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATCCCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTa                                                                                                                                                                                  <  1:348397/251‑1 (MQ=255)
                                                           tAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAgg                                                                                                                                                                                >  1:298311/1‑251 (MQ=255)
                                                                 tAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGCCAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGATCaa                                                                                                                                                                           >  1:421604/1‑250 (MQ=255)
                                                                 tAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCAGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACACTCGTGGTGATACCATGCAACCCCCGCCATAGCCGTCCTATCACGCTTATAAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCaa                                                                                                                                                                           >  1:561099/1‑250 (MQ=255)
                                                                       cTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTcc                                                                                                                                                                      >  1:162482/1‑249 (MQ=255)
                                                                                          aaagaccGCCCGGCGCAGGTGTCGCGCACAGACGCCCTATCTGCTTTTTGGGGAATCCCGAAACCCCCCCCACCCAGCCGTTGCTTTAATCCAAACCGCCCGCAGCTGCAGAGGCTTTCACGGTCC‑‑CCCACCGCCCTTGGTAATACCGGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGt                                                                                                                                                  <  2:467547/244‑1 (MQ=255)
                                                                                            agaccGCCCAGCGCAGGTTTCAGCGACAGCCGACCTATCTGCATTTTGTGGAATCCCGATGCCTCCCCCCCCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCCCGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGTTCAg                                                                                                                                               <  1:189156/247‑1 (MQ=255)
                                                                                              aaCGCCCACCCCAGGTTTCAGCGACAGACGACCTACCTGATTTTTGTGAAATACCGATAGCCCCCCCCCCCAGCCGTTGCGTCAACCCACACCGTCCGAAGCGGCGGAGGTTTTCTCGTTCC‑‑CCCAACGCTCGTGGGAATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAAGTCCAGCGGCATCCAGGATGCGGTCAGGc                                                                                                                                             <  2:122506/251‑1 (MQ=255)
                                                                                                   ccAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGCCCGCAGCTGCAGCGGCTTTCACGGTTC‑‑CCCAACGCTCGGGGTGATACCAGGCAACCCCCGCCATCGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTGGGGTCAAATTCCAGCGGCATCCAGGATGCgacacggctaagt                                                                                                                                        >  2:232635/1‑239 (MQ=255)
                                                                                                                cggCGACAGACGCCCTATCGTCATTTGGGGGAATACCGATAGCTCCCCCCCCCAGCCGTTGCTGTTACCCAGACCGTCCGCAGCTGCAGAGCCTTCCACGTTTC‑‑CCCACCGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCAc                                                                                                                           <  1:112500/249‑1 (MQ=255)
                                                                                                                cAGCGACAGACGACCTATCTGCATTTTGTGGACTCCCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCACATTCCAGCGGCATCCAGGCTGCGGTCAGGCTAAGTTCATTTTTTTCAc                                                                                                                           >  2:550183/1‑251 (MQ=255)
                                                                                                                      acaacccaccttccttctttttGGGGAATACCGATAGCCCCCCCCCCCAGCCGTTGCTGTTATCCAGAACTTCCGCAGCTGCAGAGCCTTCCACGGTCC‑‑CCCACCGCTCGTGGTGATCCCGGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGc                                                                                                                      <  1:232635/233‑1 (MQ=255)
                                                                                                                         gacgacCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTACGTTCATTTTTTTCACTTTGCGAc                                                                                                                   >  2:207942/1‑251 (MQ=255)
                                                                                                                            gacCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCTAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGAc                                                                                                                   >  1:405573/1‑248 (MQ=255)
                                                                                                                                         ttttgtagaacccagacaGCTCCCCCCCCCACCGGTCGCTTTTATCCAGAACCCCCGCAGCTTAAGAGGTTTCCACGGTCC‑‑CCCAACGCTCGTGTTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGTCAGAATAAAACCTaa                                                                                                     <  2:298311/233‑1 (MQ=255)
                                                                                                                                         tttGTGGAATACCGATATCCCCCCCACCCACCCGTTGCGGTTATCCAGACCGCCCGCAGCTGCAGAGGTTTCCCCGGTCC‑‑CCCAACGCTCGTGGTGATTCCAGGAACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAAc                                                                                                    <  1:579596/250‑1 (MQ=255)
                                                                                                                                            tgtgaataccggatgtccccccccccCAGCCGTGGCTGTAACCCGAAACGTCCGCACCTGAAGAGGTTTACACGGTCC‑‑CCCAACGCTCGTGGTGTTACCAGGCACCCCCCGCCATAGCCGTACGATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGt                                                                                                  <  2:12277/232‑1 (MQ=255)
                                                                                                                                              ggAATACCGATAGCTCCCCCCCCCAGCCGTTGCTGTTATCCAGACCGTCCGCAGCTGCAGAGGTTTTCACGGTCC‑‑CCCAACGCTCGTGTTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGt                                                                                               <  1:148575/250‑1 (MQ=255)
                                                                                                                                                  atcccGATACCTCCCCCCCCCAGCCTTTGTTGTTATCCAGAACGTCCTCAGCTGCAGAGGCTTCCGCGTTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCTCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGa                                                                                           <  1:43463/248‑1 (MQ=255)
                                                                                                                                                   agcctaTAGCCCCCCCGCCCACCCGTTGTTGTAATCCAGAACGTCCGCATCTGCAGAGGCTTTCACGGTCC‑‑CCCAACGCTCGGGGGGAGACGAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCGTTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGAc                                                                                          <  1:246169/246‑1 (MQ=255)
                                                                                                                                                          cgccccccccccCGGCCGTTCCGGTTATCCAGAGCGTCCGCAGCTGCAGAGGCTTTCCCGGTTC‑‑CCCAACGCTCGTGGTTATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATa                                                                                   <  1:311753/247‑1 (MQ=255)
                                                                                                                                                                cccccccaGCCTTGGCTGTTATCCAGAACGCCCGCAGCTGCAGAGGCTTCCCCGTTTC‑‑CCCACCGCCCGTGGGGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGt                                                                             <  1:410955/247‑1 (MQ=255)
                                                                                                                                                                        gCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCtt                                                                     >  1:237338/1‑251 (MQ=255)
                                                                                                                                                                           gTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAAAGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCAt                                                                  >  1:551018/1‑251 (MQ=255)
                                                                                                                                                                              gCTGTTATCCAGAACGTCCCAAGCGGCAGAGGCTTTCCCGGTTC‑‑CCCAACGCTGGGGGTGATACCCGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAAc                                                               <  2:550629/251‑1 (MQ=255)
                                                                                                                                                                               cTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACt                                                              >  1:84060/1‑251 (MQ=255)
                                                                                                                                                                                 gTTATCCAGACCTCCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGTTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTaa                                                            <  1:93376/251‑1 (MQ=255)
                                                                                                                                                                                 gTTATCCGAACCGTCCGCAGCTGCAGAGGCTTTCACGTTCC‑‑CCCAACGCTCGTGTTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACt                                                              <  1:232568/249‑1 (MQ=255)
                                                                                                                                                                                      cccAGAACGTCCGCAGCTGCAGAGGCTTTCCCGGTCC‑‑CCCAACGCTCGTGGTGATACCAGGAAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAg                                                        <  1:123694/250‑1 (MQ=255)
                                                                                                                                                                                                  cccAGCTGCAGAGGTTTTCACGGTTC‑‑CCCAACGCTCGGGTTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCt                                             <  1:207942/249‑1 (MQ=255)
                                                                                                                                                                                                    gagctccattggcttccccGGTCC‑‑CCCAACGCCCGTGGGTATACCAGTCACCCCCCGCCATAGCCGTACTATCACGCGTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGTAATTGCCTTc                                           <  2:251427/233‑1 (MQ=255)
                                                                                                                                                                                                      gCTGCAGAGACTTTCACGGCTCAACCCA‑‑GCTCGTGGTGATACCAGGCAACCACCGCCCTAGTCGTCCTATCCCGCTTATCAAGCTGTTGCATGGCATAATCGTAATTATCGGTAGGGTCAAATTCCAGCGGCATCCCGGATACGGTCCTGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACATAACGTGGTCTGACGAACATACAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGt                                         >  2:534880/1‑251 (MQ=255)
                                                                                                                                                                                                                  ttCACGTTCC‑‑CCCAACGCTCGTGGTGTTCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATATCTAATCAGCGGAATTGCCTTCGTATTTTCATta                               <  1:347103/249‑1 (MQ=255)
                                                                                                                                                                                                                      cacgtttcccacaCGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGTCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATGTGCATTATaa                            <  1:175019/238‑1 (MQ=255)
                                                                                                                                                                                                                           c‑‑cccACCGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATa                    <  2:535935/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aaCGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGc                >  2:472985/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGCGATTc           <  1:313822/250‑1 (MQ=255)
                                                                                                                                                                                                                                               acccggtcAACACCCGCCATATCCGTACTATCACGCTTATCAAGCTGGTGCATGGCATAATCGTCATTATCTGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATACAACCTAACGTGGTTTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGCGATTCATTATAAGc  <  1:9964/244‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
CATCCTGCGCGGAATAACTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCGACAGACGACCTATCTGCATTTTGTGGAATACCGATAGCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGCGATTCATTATAAGC  >  NC_000913/3751295‑3751778

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