Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913815,8400TC26.6% 116.5 / 10.3 66I270V (ATT→GTT) ybhKputative transferase YbhK
Reads supporting (aligned to +/- strand):  ref base T (34/13);  new base C (0/17);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 8.68e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGTATAAAAACTGCCAGG  >  NC_000913/815598‑816080
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
cccttaaccggCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTGGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACGGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGctct                                                                                                                                                                                                                                           <  2:172915/241‑1 (MQ=255)
                   ttCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCAc                                                                                                                                                                                                                       >  1:307798/1‑251 (MQ=255)
                                ttggAATTTGCGGGGTTTGTGGCCCGGATAAGGGGTTCCCGCCGCACCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTCCCTTCGGATTCCCGTTACCCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCGGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCcgac                                                                                                                                                                                                           <  1:214904/247‑1 (MQ=255)
                                  tGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGACTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGAt                                                                                                                                                                                                        >  1:569534/1‑251 (MQ=255)
                                           gcgGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGCCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCCCCTGG‑ATCACAATCCGCCCTTTCACCGCCGCGACATCCACTTTTGGCCCGACGATGAcaccaac                                                                                                                                                                                               >  2:128507/1‑245 (MQ=255)
                                                  gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGAcc                                                                                                                                                                                        >  2:530167/1‑251 (MQ=255)
                                                   taagCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCt                                                                                                                                                                                       >  1:456027/4‑251 (MQ=255)
                                                         ggATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTAc                                                                                                                                                                                 >  1:237952/1‑251 (MQ=255)
                                                          gATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTAcc                                                                                                                                                                                >  1:117184/1‑251 (MQ=255)
                                                                 ggtccccgccgcATCCGGCATAAACACCGCCCCCTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTCCCGTTAACCTAAAGCCGGTAACGCC‑TTTCCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCGGGCCTCCAT‑TACCCCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATAc                                                                                                                                                                         <  1:8129/245‑1 (MQ=255)
                                                                  gTTCACGCCGCATCCGGCATAAACAACGCGACCTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTACCCTAAAGCCTGAAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCAGGACG‑ATACGGAATATCGCGGGCCTCCAG‑TACCTCCTGG‑ATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACt                                                                                                                                                                        <  1:187996/251‑1 (MQ=255)
                                                                        gcctcATCCAGCATAAACAACGCGCACTTTG‑TATTCAATAGGCACCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTCAC‑AACTGGCGGTCATGACG‑ATACGGTAAATCGCTGGCCTCCAG‑TACCTCCTGG‑TTCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGAAGTCGACAGCATCCATGACCTCTTTACCAACATACttcccc                                                                                                                                                                   >  2:48509/5‑245 (MQ=255)
                                                                         ccgcaccgcgaaTAAACACCGCTCACTTTC‑TCTTCAAACGGACCCCGAAATTTTCCTTCGATTTTCCGTAACCCAAAAGCCT‑TAACCCCCTTTACC‑AGCGCTTTAGGTAAC‑ACCGGCCGTTCATGACG‑ATCCGGATTATCGCTG‑CCCCCAGGTACC‑CCTGGGACCCCATTCCCCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCAt                                                                                                                                                                 <  2:139383/240‑1 (MQ=255)
                                                                         ccccccccgtaaTAAACAACGCG‑TCTTTGGTCTCCAATGGGAAGCTGGATTTTTCCTCCGGATTCCCTTAAACCAAACGCCTGTAACGCT‑TTTCCC‑CCCGCGTAATGAAAC‑AACGGGCGGTCAGGAGC‑ATCCGAAGTATCGCGGGCCCCCAT‑TCCCTCCTGG‑ACCCCAACCCCCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCAt                                                                                                                                                                 <  2:237952/240‑1 (MQ=255)
                                                                            cATCCGGCATAAACACCGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTACCCTAAAGCTTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TCCCTCCTGG‑ATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAAt                                                                                                                                                              <  1:530167/251‑1 (MQ=255)
                                                                              cccGGCATAAACAACGCGCTCTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCTTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TCCCTCCTGG‑ATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTg                                                                                                                                                            <  1:168271/250‑1 (MQ=255)
                                                                              cccGGCATAAACAACGCGCCCTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTACCCTAAACCCTGTAACGCC‑TTTTCC‑AGCGCGTTATTTAAC‑ACCTGGCGGTCAGGACG‑ATACGGAATATCGCTGGCCTCCAG‑TCCCTCCTGG‑TTCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTg                                                                                                                                                            <  1:25057/250‑1 (MQ=255)
                                                                                cGGCATAATCAACGCGCACTTTG‑TTTTAAATAGGACGCAGGAATTTTCCTTCGGTTTTCCGTTAACCTAAACCCTGTCACGCT‑TTTTCC‑AGCGCGTTATGTACC‑AACTGGCGGTCATGACT‑AGTCGGAATATCGCTGCCCCCCAG‑TCCCCCCTTG‑TTCCCATCCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGcc                                                                                                                                                          <  1:85050/251‑1 (MQ=255)
                                                                                    aTAAACAACGCGCACTTTG‑TCTCCAATAGGAAGCCGGAATTTCCCTTCGGATTCCCGTTAACCTAAAGCCTGTAACCCC‑TTTTCC‑ACCGCTTTATGTAAC‑AACTGGCGGTCATGCCG‑ATACGGAATATCGCTGGCCTCCAG‑TACCCCCTGG‑ATCACAACCCGCTCTTCCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAgct                                                                                                                                                      <  1:23329/251‑1 (MQ=255)
                                                                                      aaaCATCGCGCACTTTG‑TCTTCAATAGGAAGCCGGATTTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTCCC‑AGCCCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TCCCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAgctt                                                                                                                                                     <  1:363121/250‑1 (MQ=255)
                                                                                       aacaacGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGATGGCCTCCAG‑TACCTCCTGG‑TTCCCAATCCGCTCTTTCCCCGCCGGGACATCCACTTTTGGCCCGACGATGACCGCATCAATTCCCTTTTTACCCACCTACTGCTCCATAATTGCCAgcttgc                                                                                                                                                   >  2:222751/1‑251 (MQ=255)
                                                                                          aacGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCAACAATGACCTTTTTACCAACATCCTGCTCCATAATTGCCAGCTTGCttt                                                                                                                                                >  1:461370/1‑251 (MQ=255)
                                                                                           ccgcgCACTTTG‑TCTTCAATAGGAAGCCGAAATTTTCCTTCGGATTTCCTTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TCCCTCCTGG‑ATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTc                                                                                                                                               <  1:438620/250‑1 (MQ=255)
                                                                                           acGCGCACTTAG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTTGCCCGACGATGACCGCATCACTGACCTTTTTACCCACATACTGCTCCATAATTGCCAGCTTGCTTTc                                                                                                                                               >  2:519305/1‑251 (MQ=255)
                                                                                              cgctttttG‑TCTTCAATGGGAGGCCGGAATTTTCCTTGGGTGTTCCGTGACCCAAAAGCCTTTAACGCT‑TTTGCC‑AGCGTGTTATTTAAC‑AACTGGCGGTCATGACG‑ATACGGAATACCGCTGGCCTCCAG‑TACCCCCTGG‑ATCAAAACCCGCCCTTCCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAg                                                                                                                                            <  2:572595/246‑1 (MQ=255)
                                                                                               cccactttgatttaaaTAGGAAGCAGGATTTTTCCTGCGGATTTCCTTTAACCAAAAGCCTTAAACGCC‑TTTTCC‑AGCTCGTT‑TGGAACCAGCGCGCGTTAAGGGCG‑ATACGAACGATCGCTGGCCCCCAG‑TCCCCCTTGG‑ATCAAAATCCGTTCTTGCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATGCTGCTCCATAATTGCCAGCTTGCTTTCCAg                                                                                                                                            <  2:550242/237‑1 (MQ=255)
                                                                                                cccTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTCCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TCCCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCt                                                                                                                                          <  1:579970/249‑1 (MQ=255)
                                                                                                cACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGAc                                                                                                                                                                                                                             >  1:10887/1‑168 (MQ=255)
                                                                                                cACTTTG‑TCATCAATAGGAATCCGGAATTTTCCTTCGGATTTCCTTTAACCAAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGTTCATGACG‑ATACGTAATATCGCTGGCCTCCAG‑TCCCCCCTGG‑ATCACACCCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCt                                                                                                                                          <  1:520379/251‑1 (MQ=255)
                                                                                                 aCTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAACGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACAt                                                                                                                                                                           >  2:36048/1‑217 (MQ=255)
                                                                                                 aCTTTG‑TCTTAAATAGGAAGCCGGATTTTTCCTTCGGATTTCCGTTACCCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TCCCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACAt                                                                                                                                                                           <  1:36048/217‑1 (MQ=255)
                                                                                                           tCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCCTGACG‑ATACGGAATATCGCTGGCCTCCAG‑AAACTCCTGG‑CTCACAAACCGCTCATTCACCGCAGAGCCATCCACTTTTTGACCGACGATGACCCCACCAATTCACTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCCGCTCCAAAtttt                                                                                                                                 >  2:337793/1‑248 (MQ=255)
                                                                                                            aaaTAGGAGGCCGAAATTTCCCTTCGGTTTCCCGTTACCCAAAACCCTTTACCGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGGGGTCATGACG‑ATACGGAATATCGCTGGCCCCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGcc                                                                                                                               <  1:533136/250‑1 (MQ=255)
                                                                                                              aTAGGAAGCCGGCATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑ACCGCGTTCTGTAAC‑AACTGGCGGTCATAACG‑ATAAGGAATATCGCTGGCCTCCAG‑CACCAACTGG‑CTCACAATCCGCTCTTTCACCCCCGGGACATCCCCTTGTGGCCCGACGATAGCCGCATCAATGACCTTTTTCACAACATCCTGCTCCATAATGCCCATCTTGCTTTCCCGCTTCAAATTAGcagc                                                                                                                             >  1:550535/1‑248 (MQ=255)
                                                                                                                   aaGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGt                                                                                                                         >  1:23083/1‑250 (MQ=255)
                                                                                                                    aGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGCTGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTa                                                                                                                        >  1:373510/1‑250 (MQ=255)
                                                                                                                    aGCCGAAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTACCGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCCCCAG‑TACCTCCTGG‑ATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTaa                                                                                                                       <  1:574082/251‑1 (MQ=255)
                                                                                                                        gaaaTTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACCCC‑TTTTCC‑AGCGCGTTATGTACC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TCCCCCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAaactc                                                                                                                   <  1:357950/249‑1 (MQ=255)
                                                                                                                           aTTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑AACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATAATGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGGAcactcaac                                                                                                                >  2:359495/1‑251 (MQ=255)
                                                                                                                              ttCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTACAGCTTCAAATTAGCCGACGGTAaactcaacgc                                                                                                              >  1:128347/1‑248 (MQ=255)
                                                                                                                                    cGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCCCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCCGGAAAACTCAACTCAcccccc                                                                                                       >  2:554565/1‑251 (MQ=255)
                                                                                                                                       aTTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGCCGATGACCGCATCACTGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGcccc                                                                                                       >  1:127071/1‑248 (MQ=255)
                                                                                                                                        tttCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGAGCGTTATGTAAC‑AACTGGCAGTCATGACG‑ATACGGAATATCGCTGTCCTCCAG‑TAGATCCTGG‑ATCACAATCCGCTCTTTAACCGCCGGGACAACCAATTTTGGCCCGACGGTTACCGCAAAAATGACCTTTTTACCGACAGACTGCACAAAAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAAGTCACGCCACAgct                                                                                                   >  1:540498/1‑249 (MQ=255)
                                                                                                                                          ccctttAACCTACAGCCGGTACCGCC‑TTTCCC‑AGCGCTTTATTTAAC‑AACTGGCGGTCATGACT‑ATACGGAATATCGCTGGCCCCCAG‑TACCCCCTGG‑ACCCCAACCCCCCTTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGAtt                                                                                                  <  1:159434/246‑1 (MQ=255)
                                                                                                                                           ccGTTACCCTAAAGCCTGTAACGCT‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCAGCAGGTAAACTCAACTCACGCCCCAGATTGc                                                                                                <  1:178909/251‑1 (MQ=255)
                                                                                                                                             gTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCg                                                                                              >  1:73048/1‑251 (MQ=255)
                                                                                                                                                aaCCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGACTATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGATTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGat                                                                                            >  1:71111/1‑250 (MQ=255)
                                                                                                                                                            gTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTNCCGCCCCAGATTGCCGATATAAACCATCGgc                                                                               >  2:287146/1‑251 (MQ=255)
                                                                                                                                                             gtaaccccttgccc‑AGCGCGTTATTAAAC‑ACCGGGCGGTCAT‑ACGGATACGGAAAACCGCTGCCCCCCAT‑TCCCTCTGGG‑ATCACATTCCCCTTTTACACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCGTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTGCAAATTACCCGCAGGTAAACGCAACTCACGCCCCAGATTGCCGATATAAACCATCgg                                                                                <  1:481761/238‑1 (MQ=255)
                                                                                                                                                              aaCGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCAGCAGGTAAACTCAACTCACGCCCCAGCTTGCCGATATAAACCATCGgcgc                                                                             >  1:381685/1‑251 (MQ=255)
                                                                                                                                                               aCGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATc                                                                                  >  1:243623/1‑245 (MQ=255)
                                                                                                                                                                actc‑TTTTCCCAGCGGGTTATTTAAC‑AACGGGGTGTCATGACG‑AGACGGAATCTCGCTGGCCCCCAT‑GCCCTCCTGG‑ATCACAACCCGCCTTTTCACCTCCGAGACAGCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTGTTACCAACATACTTCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATc                                                                                  <  2:243623/242‑1 (MQ=255)
                                                                                                                                                                 gCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGcc                                                                                                                                                                                                               >  1:49127/1‑118 (MQ=255)
                                                                                                                                                                      tttCC‑AGCGCGTTATGTAAC‑AACTGGCGTTAATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCCCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGt                                                                       <  1:173562/250‑1 (MQ=255)
                                                                                                                                                                                    aTGTAAA‑AACGGGCGGTCATGACG‑ATACGGAATATCCCTGGCCTCCAG‑TACCCCCTGG‑ATCCCAACCCGCTCTTCCCCCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCt                                                         <  1:207223/251‑1 (MQ=255)
                                                                                                                                                                                                         gACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTcagcagc                                     >  1:550427/1‑251 (MQ=255)
                                                                                                                                                                                                            g‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCg                                                                                              >  1:451604/1‑191 (MQ=255)
                                                                                                                                                                                                            g‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTcagcagc                                     >  2:92126/1‑248 (MQ=255)
                                                                                                                                                                                                            ccaTACGTAACATCCCTGGCCCCCAC‑TACCTCCGGG‑ACCACAACCCGTCCTTTCCCCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTcagcagcag                                   <  1:20166/249‑1 (MQ=255)
                                                                                                                                                                                                                  ggAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTgg                             >  1:293033/1‑251 (MQ=255)
                                                                                                                                                                                                                       atCGCTGGCCTCCCG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAAt                                                             >  1:399448/1‑214 (MQ=255)
                                                                                                                                                                                                                             ggCCTCCAG‑TACGTCCTGG‑ATCACAAACCGCTCTTTCACAGCCGAGACATCCACTTTTGGACCGACGATGACCGCATCAATGACCTTTTTACCAACATCCTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGCAAACTCAACTCAAGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGCCGAAATGCCTGGTCGATTTCCTTCAGCAGCAGCATTGGCATCCGGcggg                  >  2:63831/1‑248 (MQ=255)
                                                                                                                                                                                                                               ccTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGt                 >  2:505659/1‑250 (MQ=255)
                                                                                                                                                                                                                                 cccAG‑TACCCCCTGG‑ATCACAACCCGCTTCTTCACCGCCGAGACATCCACTTTTGGCCCTACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGtata              <  1:45282/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cctcctGG‑ATCACAACCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGTATAAAAACTg       <  1:296340/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGTATAAAAACTGCCAgg  >  1:238161/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
CGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTG‑TCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCC‑TTTTCC‑AGCGCGTTATGTAAC‑AACTGGCGGTCATGACG‑ATACGGAATATCGCTGGCCTCCAG‑TACCTCCTGG‑ATCACAATCCGCTCTTTCACCGCCGAGACATCCACTTTTGGCCCGACGATGACCGCATCAATGACCTTTTTACCAACATACTGCTCCATAATTGCCAGCTTGCTTTCCAGCTTCAAATTAGCCGCAGGTAAACTCAACTCACGCCCCAGATTGCCGATATAAACCATCGGCGCTGGCGTGCGGCGTAATGCCTGGGCGATTTCCTTCAGCAGCAGAATTGGCATCAGGCTGGTATAAAAACTGCCAGG  >  NC_000913/815598‑816080

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