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.

CATCCTGCGCGGAATAACTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTAATCAGCGGAATTGCCTTCGTATTTTCATTATAACCGCGATAAGGCGATTCATTATAAGC  >  NC_000913/3751295‑3751778
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
cATCCTGCGCGGAATAACTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGcc                                                                                                                                                                                                                                           <  2:453723/251‑1 (MQ=255)
              aaaCTTGCCAGACCGCTACGACGGGACTCCCTCTCTGAAATCCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATCCCGATAGCTCCCCCCCCCAGCCGTTGTTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATc                                                                                                                                                                                                                             <  1:328866/250‑1 (MQ=255)
              aaaCTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTAAAATTCTCGTCATAATAGAGAACGCCCAGCGCGGGTGTCAGCG‑ACAGACGACCTATCGGCATTTTGTGG‑AATACCGATAGCTCCCCCCCCCAGCCGTTGCTGT‑TATCCAGACC‑GTCCGCAGCTGCAGAGGCTTTCACGTTTC‑‑CCAAACGCTCGTGGTAATACCAGTAACCCCCCGCCATAGCCGTACTATc                                                                                                                                                                                                                             <  1:114842/250‑1 (MQ=255)
                 cTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGCTACCAGGCAACCCCCGCCATAGCCGTACTATCACg                                                                                                                                                                                                                          >  1:251427/1‑251 (MQ=255)
                    gCCAGACCGCTCCGAGGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCGGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGA‑AATACCGATAGCTCCCCCCCCCAGCCTTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCt                                                                                                                                                                                                                        <  1:509097/250‑1 (MQ=255)
                                          cactcTCTGAAAGGCCATAGTAATAGCCCCTGAATTTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGCCCTATCTGCATTTTGTGG‑AATCCCGAAAGCCCCCCCCCCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATa                                                                                                                                                                                                  <  1:206305/250‑1 (MQ=255)
                                                    aaTGCCATATTAATAGCCACTGAATCTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGCCGCCCTATCTTCATTT‑GTGGGAATCCCGATAGCCCCCCCACCCAGCCGTTGCTTT‑TACCCAGAAC‑GTCCGCAGCTGCAGAGGCTCTCACGTTTC‑‑CCCAACGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTAt                                                                                                                                                                                       <  2:481241/251‑1 (MQ=255)
                                                     aTGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑CCGGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCCCCCAGCCGTTGCTTT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGTTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATc                                                                                                                                                                                      <  1:507478/251‑1 (MQ=255)
                                                         cATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                                                                                                                   >  2:112500/1‑250 (MQ=255)
                                                         cATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGt                                                                                                                                                                                   >  2:192100/1‑250 (MQ=255)
                                                         cATAGTAATAGCAACTGAAATTCTCGTCAAAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATCCCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTa                                                                                                                                                                                  <  1:348397/251‑1 (MQ=255)
                                                           tAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAgg                                                                                                                                                                                >  1:298311/1‑251 (MQ=255)
                                                                 tAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑CCAGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGATCaa                                                                                                                                                                           >  1:421604/1‑250 (MQ=255)
                                                                 tAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCAGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACACTCGTGGTGATACCATGCAACCCCCGCCATAGCCGTCCTATCACGCTTATAAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCaa                                                                                                                                                                           >  1:561099/1‑250 (MQ=255)
                                                                       cTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTcc                                                                                                                                                                      >  1:162482/1‑249 (MQ=255)
                                                                                          aaagaccGCCCGGCGCAGGTGTC‑GCGCACAGACGCCCTATCTGCTTTTTGGGG‑AATCCCGAAACCCCCCCCACCCAGCCGTTGCTTT‑AATCCAAACC‑GCCCGCAGCTGCAGAGGCTTTCACGGTCC‑‑CCCACCGCCCTTGGTAATACCGGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGt                                                                                                                                                  <  2:467547/244‑1 (MQ=255)
                                                                                            agaccGCCCAGCGCAGGTTTCAGCG‑ACAGCCGACCTATCTGCATTTTGTGG‑AATCCCGATGCCTCCCCCCCCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCCCGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGTTCAg                                                                                                                                               <  1:189156/247‑1 (MQ=255)
                                                                                              aaCGCCCACCCCAGGTTTCAGCG‑ACAGACGACCTACCTGATTTTTGTGA‑AATACCGATAGCCCCCCCCCCCAGCCGTTGC‑GTCAACCCACACC‑GTCCGAAGCGGCGGAGGTTTTCTCGTTCC‑‑CCCAACGCTCGTGGGAATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAAGTCCAGCGGCATCCAGGATGCGGTCAGGc                                                                                                                                             <  2:122506/251‑1 (MQ=255)
                                                                                                   ccAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GCCCGCAGCTGCAGCGGCTTTCACGGTTC‑‑CCCAACGCTCGGGGTGATACCAGGCAACCCCCGCCATCGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTGGGGTCAAATTCCAGCGGCATCCAGGATGCgacacggctaagt                                                                                                                                        >  2:232635/1‑239 (MQ=255)
                                                                                                                cggCG‑ACAGACGCCCTATCGTCATTTGGGGG‑AATACCGATAGCTCCCCCCCCCAGCCGTTGCTGT‑TACCCAGACC‑GTCCGCAGCTGCAGAGCCTTCCACGTTTC‑‑CCCACCGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCAc                                                                                                                           <  1:112500/249‑1 (MQ=255)
                                                                                                                cAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑ACTCCCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCACATTCCAGCGGCATCCAGGCTGCGGTCAGGCTAAGTTCATTTTTTTCAc                                                                                                                           >  2:550183/1‑251 (MQ=255)
                                                                                                                      acaacccaccttccttctttttGGGG‑AATACCGATAGCCCCCCCCCCCAGCCGTTGCTGT‑TATCCAGAAC‑TTCCGCAGCTGCAGAGCCTTCCACGGTCC‑‑CCCACCGCTCGTGGTGATCCCGGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGc                                                                                                                      <  1:232635/233‑1 (MQ=255)
                                                                                                                         gacgacCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTACGTTCATTTTTTTCACTTTGCGAc                                                                                                                   >  2:207942/1‑251 (MQ=255)
                                                                                                                            gacCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCTAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGAc                                                                                                                   >  1:405573/1‑248 (MQ=255)
                                                                                                                                         ttttgtagaacccagacaGCTCCCCCCCCCACCGGTCGCTTT‑TATCCAGAAC‑CCCCGCAGCTTAAGAGGTTTCCACGGTCC‑‑CCCAACGCTCGTGTTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGTCAGAATAAAACCTaa                                                                                                     <  2:298311/233‑1 (MQ=255)
                                                                                                                                         tttGTGG‑AATACCGATATCCCCCCCACCCACCCGTTGCGGT‑TATCCAGACC‑GCCCGCAGCTGCAGAGGTTTCCCCGGTCC‑‑CCCAACGCTCGTGGTGATTCCAGGAACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAAc                                                                                                    <  1:579596/250‑1 (MQ=255)
                                                                                                                                            tgtgaataccggatgtccccccccccCAGCCGTGGCTGT‑AACCCGAAAC‑GTCCGCACCTGAAGAGGTTTACACGGTCC‑‑CCCAACGCTCGTGGTGTTACCAGGCACCCCCCGCCATAGCCGTACGATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGt                                                                                                  <  2:12277/232‑1 (MQ=255)
                                                                                                                                              gg‑AATACCGATAGCTCCCCCCCCCAGCCGTTGCTGT‑TATCCAGACC‑GTCCGCAGCTGCAGAGGTTTTCACGGTCC‑‑CCCAACGCTCGTGTTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGt                                                                                               <  1:148575/250‑1 (MQ=255)
                                                                                                                                                  atcccGATACCTCCCCCCCCCAGCCTTTGTTGT‑TATCCAGAAC‑GTCCTCAGCTGCAGAGGCTTCCGCGTTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCTCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGa                                                                                           <  1:43463/248‑1 (MQ=255)
                                                                                                                                                   agcctaTAGCCCCCCCGCCCACCCGTTGTTGT‑AATCCAGAAC‑GTCCGCATCTGCAGAGGCTTTCACGGTCC‑‑CCCAACGCTCGGGGGGAGACGAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCGTTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGAc                                                                                          <  1:246169/246‑1 (MQ=255)
                                                                                                                                                          cgccccccccccCGGCCGTTCCGGT‑TATCCAGAGC‑GTCCGCAGCTGCAGAGGCTTTCCCGGTTC‑‑CCCAACGCTCGTGGTTATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATa                                                                                   <  1:311753/247‑1 (MQ=255)
                                                                                                                                                                cccccccaGCCTTGGCTGT‑TATCCAGAAC‑GCCCGCAGCTGCAGAGGCTTCCCCGTTTC‑‑CCCACCGCCCGTGGGGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGt                                                                             <  1:410955/247‑1 (MQ=255)
                                                                                                                                                                        gCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCtt                                                                     >  1:237338/1‑251 (MQ=255)
                                                                                                                                                                           gTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAAAGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCAt                                                                  >  1:551018/1‑251 (MQ=255)
                                                                                                                                                                              gCTGT‑TATCCAGAAC‑GTCCCAAGCGGCAGAGGCTTTCCCGGTTC‑‑CCCAACGCTGGGGGTGATACCCGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAAc                                                               <  2:550629/251‑1 (MQ=255)
                                                                                                                                                                               cTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACt                                                              >  1:84060/1‑251 (MQ=255)
                                                                                                                                                                                 gT‑TATCCAGACC‑TCCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑CCCAACGCTCGTGTTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACTaa                                                            <  1:93376/251‑1 (MQ=255)
                                                                                                                                                                                 gT‑TATCC‑GAACCGTCCGCAGCTGCAGAGGCTTTCACGTTCC‑‑CCCAACGCTCGTGTTGATACCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAGGATGCGGTCAGGCTAAGTTCATTTTTTTCACTTTGCGACAGAATAAAACCTAACGTGGTCTGACGAACATAAAAAGTATCACCTTCATAACt                                                              <  1:232568/249‑1 (MQ=255)
                                                                                                                                                                                      cccAGAAC‑GTCCGCAGCTGCAGAGGCTTTCCCGGTCC‑‑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)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
CATCCTGCGCGGAATAACTTGCCAGACCGCTACGACGGGACTCACTCTCTGAAATGCCATAGTAATAGTCACTGAAATTCTCGTCATAATAGAGAACGCCCAGCGCAGGTGTCAGCG‑ACAGACGACCTATCTGCATTTTGTGG‑AATACCGATAGCTCCCCCACCCAGCCGTTGCTGT‑TATCCAGAAC‑GTCCGCAGCTGCAGAGGCTTTCACGGTTC‑‑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: ‑