Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913825,2350TA32.3% 81.5 / 19.0 64M292L (ATG→TTG) ybhRputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (34/8);  new base A (0/20);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 3.38e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCC  >  NC_000913/825010‑825448
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                         
ccTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCCTA‑AAGACAAaaac                                                                                                                                                                                                              >  1:335071/1‑238 (MQ=255)
          aaaCATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGGAAAGTGGCGAATAGGGGT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGGACCACA‑ACGACAAACACGCCCATaaaa                                                                                                                                                                                                    >  2:994431/1‑237 (MQ=255)
            aCATCGCGTGCGCCGATGCCCCTGTCGTGGCCGTTATCCCCAGTAGCGGCCACAAACTTTCCCACACAATATCCAGACTCGCTTCCTTCAAATAAATCTGCTGGGTAAGGCCCGAAAAGTGGCAAATAGGGTT‑AACCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTCCCAC‑CGAAGAAACG‑TAACCGGAAAGGAGAATGGGGGGCATCAAA‑AAGACAAACACGCCGATAAACGcc                                                                                                                                                                                                 <  2:1513112/239‑1 (MQ=255)
                cgcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAAAGCccgcc                                                                                                                                                                                             >  1:298078/1‑238 (MQ=255)
                 gcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAAACAAAAACGCCGAAAAACACct                                                                                                                                                                                                >  1:1424218/1‑235 (MQ=255)
                   gTACGCCGCTGACCATGCCGGGGCCTTTATCACTGGTGGCGCCAACAAACTATCCCAAACAATACCCAGGCACGCATCCTTCAAAAAATTCTGCTTGTTAAGGTCAGTAAAGTGGCAAATAGGGTT‑AACCCACGCCAGTTTTT‑GCAGCCAAACCGGCATTTTTTCCAC‑CGGAGAAACG‑TAACGGAAAGGGAAAAGGGCGGGCATCAAA‑AAGACAAACACGCCGATAAACGCCTGCtgtt                                                                                                                                                                                          <  2:619837/239‑1 (MQ=255)
                      cgccgcTGCCCCTGTGTTGGCCGTTCTCACCAGTAGCGGCCACAAATTATCCCAAACAATATCCAGACTCGCATCCTTCAAAAAAATCTGTTGGGTATGGCCCGAAAAGTGGCGAAGGGGGTA‑AATCCACGTCGGGTTTT‑GCAGCAATACCGGCATGTTTCCCAC‑CGGAGAAACG‑TAACCGGAAGGGGGAGGGGCGGGCATCATA‑AAGAAAAACACGCCGATAAACGCCTGCtgttgtg                                                                                                                                                                                       <  2:1178858/239‑1 (MQ=255)
                                   gTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAAAAAAACACGCCGATAAAAGCCTGCTTTTGGGTTTGACAAAGtga                                                                                                                                                                          >  1:34693/1‑239 (MQ=255)
                                        ggCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAAAAAACCGATAAAAGCCTGCTGTTGTGTTGAAAAGAGTGAGGaaa                                                                                                                                                                     >  1:1244594/1‑239 (MQ=255)
                                                cccccAGAAGCGACCAAAAACTATCCCACAAAATACCCAGACTCGCATCCTAAAAATAAATCTGCTTGGTAATGTCAGAAAAGGGGCAAATAGGGTT‑AAACCACGTCAGGTTTT‑GCAGCCATACCGGCATTTTTCCCAC‑CGGAGAAACG‑AAACCGAAAAGGAAAAGGGCGGGCATCAAA‑AAGACAAACACGCCGATAAACGCctgctg                                                                                                                                                                                            <  2:764513/205‑1 (MQ=255)
                                                       tAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAAAACCCCGATAAACGCCCGCTGTTGTGGTGAACAGAGTGGTGAAAACAACAGaaccaaa                                                                                                                                                       >  2:312357/1‑232 (MQ=255)
                                                          cGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCCATAAACGCCTGCTGTTTTGTTGaa                                                                                                                                                                                  >  1:764513/1‑208 (MQ=255)
                                                            gCCACAAACTATTCCACACAATATCCAGACTCGCCTCCTTAAAATAAATCTGCCTGGTAATGTCCGTAAAGTGGCGAATAGTGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCCGCCTGTTGTCCAC‑CGGAGAAACG‑CAACCGGTAAGGCGAATGGCGGGGTTCAGACAAGCCAAAAAGCCCAGAAAACCCCGGCTGTTGGGTGGAACAGCG‑GATGAAATAAAAAGaaccaaatcaca                                                                                                                                                  >  2:991239/1‑227 (MQ=255)
                                                                        tCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAAAACGCCGAAAAACGCCTGCTGTTGGGTTGAACAGAGGGATGAAATCAAAAGACCGAATCCCACCAACGAT‑AAACCa                                                                                                                                    >  1:723068/1‑239 (MQ=255)
                                                                                 tatCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAACGAT‑AAACCATAAATCAcc                                                                                                                           >  1:980471/1‑239 (MQ=255)
                                                                                 tatCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAAAACGCCGATAAAAGCCTGCTGGTGGGTTGAAAAAAGGGATGAAAACAACAAACCGAATCCCACCAACGGT‑AAAACATAAAACAcc                                                                                                                           >  1:295339/1‑239 (MQ=255)
                                                                                         cTCCCATCCTTAAAATAAATCTGTTGGGTAATTTCCGAAAAGTGGCAAATAGGTTT‑AATCCACGTCAGGTTTT‑GCAGCCATCCCGGCATTTTTTCCAC‑CGGAAAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCAAA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCg                                                                                                                       <  2:295339/235‑1 (MQ=255)
                                                                                         atgcaaTCCTTAAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCCCGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TACCCGGAAAGGGGAATGGCGGGCATAATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCg                                                                                                                       <  2:1384347/230‑1 (MQ=255)
                                                                                           cGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAAAGCCTGCTGTTGTGTTGAAAAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATAACCATCGTAAAGt                                                                                                                 >  1:1603637/1‑239 (MQ=255)
                                                                                           cGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAAAGCCCGCTGTTGTGTTGAACAGAGTGATGAAATCAAAAGACCGAATCCCACCAGCGAT‑AAACCATAAAACACCATCGGaaaat                                                                                                                 >  1:82084/1‑237 (MQ=255)
                                                                                             cATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGCTAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCCTCGTAAAGTAg                                                                                                               >  1:266504/1‑239 (MQ=255)
                                                                                                 cTTAAAAAAAATCGGCTTGTAAATGCCCG‑AAAGGGGGCAAAAGGGTTAAATCCCGGTCAGGTTTT‑GCAGCCATACCGGCATTTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCAAA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAAt                                                                                                                                                       <  2:1433664/196‑1 (MQ=255)
                                                                                                     aaataaatCTGTTGGTTAATGCCCGAAAAGTGGCAAAGAGGTTA‑AACCCACGT‑AGGTTTTGGCAGCCATACCGCATTGTTTTCCAC‑CGGAAAAACG‑TACCGGGAAAGGAAAGGGGGGGGCATCAAA‑AAAACAAACACGCCGATAAACGCctg                                                                                                                                                                                               <  2:683735/152‑1 (MQ=255)
                                                                                                       ataaatCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAAACAAACACGCCCATAAAAGCCCGCtgttgtg                                                                                                                                                                                       >  1:748194/1‑158 (MQ=255)
                                                                                                           atCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACAAGCCGATAAAAGCCCGCTGGTGTGTTGAACAGAGGGATTAAATCAACAGACCGAATCCCACCAGCg                                                                                                                                             >  1:1433664/1‑196 (MQ=255)
                                                                                                            tCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGGATCATA‑AAGACAAAAACGCCGATAAAAGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATACCACCAGCGAT‑AAAACATAAAACACCAACGTAAAATAGAAAAACAC‑CCGCGa                                                                                                >  2:379544/1‑238 (MQ=255)
                                                                                                               gCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAAAGCCTGCtgttgtgtt                                                                                                                                                                                     >  1:683735/1‑152 (MQ=255)
                                                                                                                cTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAAAGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑aaa                                                                                                                                       >  1:474189/1‑196 (MQ=255)
                                                                                                                         gTCCGTAAAGTGGCAAATAGGGTT‑AATCCACGTCGGGTTTT‑GCAGCCATACCGGCATGTTTGCCAC‑CGGAGAAACG‑TAACCGGAAAGGAAAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCCCCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGa                                                                                  <  2:337994/239‑1 (MQ=255)
                                                                                                                          tCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAAAGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAAACATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGAt                                                                                 >  1:1500213/1‑239 (MQ=255)
                                                                                                                            cgaaaaGTGGCAAATAGGGTA‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAAAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGAAAAACACGCCGATAAACGCctg                                                                                                                                                                                               <  1:31675/126‑1 (MQ=255)
                                                                                                                                agtgggCGAAGAGGTTA‑AATCCAGTTCAGTTGTT‑GCACCCATCCCGGCATGTTTTCCAC‑GGGAGAACCG‑TAACCGGAAGGGAGAATGGCGGGCATCAAA‑AAAACAACCACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATa                                                                           <  2:1244594/235‑1 (MQ=255)
                                                                                                                                   tGGCAAATAGGGTT‑AATCCCCGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAAAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCAAA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCccac                                                                                                                                                  <  1:668436/167‑1 (MQ=255)
                                                                                                                                     gagAATAGGGTA‑ATACCTCGTAGGGTTTG‑GCAGCCAAACCGGCATTTTTTCCAC‑CGGAGAAACG‑AAACGGGAAAGGAGAATGGGGGGCATCATA‑AAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGcc                                                                       <  2:1424218/236‑1 (MQ=255)
                                                                                                                                      cGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAAAACCCCGATAAACGCCTGCtgttgtgtt                                                                                                                                                                                     >  2:31675/1‑129 (MQ=255)
                                                                                                                                             ggTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGGGATGAAAACAACAGACCGAATCCCCCCAGCGAT‑aaa                                                                                                                                       >  2:668436/1‑167 (MQ=255)
                                                                                                                                                   acccccGTCGGGTTTT‑GCAGCCATCCCGGCATGTTTCCCAC‑CGGAAAAACG‑TAACCGGAAAGGAAAAGGGCGGGCATCAAA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGc                                                        <  2:249888/234‑1 (MQ=255)
                                                                                                                                                   aTCCACTGCAGTTTTA‑GCAGCGATGCCGACATGTTTTCCAC‑CGGAAAAACG‑TAACCGAAAGGGAAAATGGCGGCCATCAAA‑AACACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGc                                                        <  2:324807/239‑1 (MQ=255)
                                                                                                                                                   aTCCACGTCAGGTTTT‑GCTGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGAAAAACACGCCGAAAAACGCCCGGTGTTTTGTTGAACAGAGGGGTGAAATCAACACACCGAAACCCACCAGCGCT‑AAAACATAAAACAACATCGCAAAGGAGAATAGCGC‑CAGCGATCCGGCGAAGGGGAATTGATACCACCC‑GATACCAAACg                                                         >  2:1555472/1‑238 (MQ=255)
                                                                                                                                                      cACGTCAGGTTTT‑GCAGCCATCCCGGCATTTTTTCCAC‑CGGAAAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCAAA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGCCAg                                                     <  1:312357/239‑1 (MQ=255)
                                                                                                                                                        cGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGGGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATAAGCCCA‑GATACCCATCGCCCGc                                                    >  1:249888/1‑238 (MQ=255)
                                                                                                                                                               ttgt‑GCAGCAATCCCGGCATTTTTCCCAC‑CGGAGAAACG‑TAACCGAAAAGGAGAATGGCGGGCATCAAA‑AAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGCCAGCACAAtgg                                             <  2:298078/235‑1 (MQ=255)
                                                                                                                                                                  t‑GCAGCCATACCGGCATGTTTCCCAC‑CGGAAAAACG‑TACCCGGAAAGGAGAATGGCGGGAATCAAA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACag                                                                                                                                                                               <  1:1616016/108‑1 (MQ=255)
                                                                                                                                                                        ccATCCCGGCATGTTTCCCAC‑GGGAAAAACG‑TACCCGAAAAGGAGAATGGCGGGAATAAAA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGCCAGCACAATGGTGGCCTGGa                                    <  1:1437289/239‑1 (MQ=255)
                                                                                                                                                                           tcccGGCATGTTTCCCAC‑CGGAAAAACG‑TAACCGGAAAGGAGAAGGGCGGGCATCAAA‑AAAAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGtg                                                                                                                                                                           <  2:1024705/102‑1 (MQ=255)
                                                                                                                                                                            cccGGCATGTTTCCCAC‑GGGAAAAACG‑TAACCGGAAAGGAGAAGGGGGGGCATAAAA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCccacca                                                                                                                                                <  2:1466855/129‑1 (MQ=255)
                                                                                                                                                                             ccGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGGTGTGTTGAACAGAGTGATGaaa                                                                                                                                                                     >  2:1616016/1‑108 (MQ=255)
                                                                                                                                                                                   agTTTTCCACGCGAAAAAACG‑GAACCGAAAAGGAGAAGGGGGGGAATCAAA‑AAAACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGa                           <  2:266504/237‑1 (MQ=255)
                                                                                                                                                                                     ttttCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCaa                                                                                                                                                                 >  1:1024705/1‑104 (MQ=255)
                                                                                                                                                                                      tttCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATT‑AAGACAAACACCCCGATTAACGCCTTTCGTTGTGTTGTAAAGGGAGATGAAAAAAAAAGAACGAATACCACCAGCGATAAAACCAT‑AATCAACCCCGTAAAATAAACCACCGCACAGGGA‑CCGGCGAAGGGGGTTTGTATCTCCCAGGATACC‑ATAGCCCGCACAAAAGAGGCCTTGAAGGTggggctcatt                      >  2:880797/1‑231 (MQ=255)
                                                                                                                                                                                      tttCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTTAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAAcg                                                                                                                                     >  1:1466855/1‑129 (MQ=255)
                                                                                                                                                                                           ccaccgaaaaaCGTTAACGGGAAAGGGAAATGGAGGGCATCAAA‑AAGAAAACCACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAAc                   <  2:723068/231‑1 (MQ=255)
                                                                                                                                                                                                gagaAACG‑TAACCGGAAAGGAGAATGGCGGGCCTCATA‑AAGACAAAAACCCCGAGAAACACCTGCCGTTGTGTTGAACAGAGTGATGAAATCAACAGACCCAATCCCACCAGGCAT‑AAACCATAAAACACCAAACGAAAGGAGAACAAGGC‑CAAGGATTCCGCGAAGGGGATTTGATACCCCCC‑GAAACCAAACGACAGCACACTAGGGGCCTGGAAGGGGGCGACAAATAAAACCGGt             >  2:206566/1‑239 (MQ=255)
                                                                                                                                                                                                     aCG‑TAACCGGAAAGGAGAATGGCGGGCATCCTA‑AAGGCAAACACCCCGATAAACGCCCGCTGTTGTGTTTAACAAAGTGATTAAATCAACCGACCGGATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAAACGCCCGCCCCATGGTGGCCTGGAGAGTGGCGACAATTAAACCCGGGAcc          >  2:793193/1‑236 (MQ=255)
                                                                                                                                                                                                      cG‑TACCCGAAAAGGAGAAGGGCGGGCTTAAAA‑AAGACAAACACGCCGATAAACGCCTGCtgtt                                                                                                                                                                                          <  2:1149615/63‑1 (MQ=255)
                                                                                                                                                                                                          aaCCGGAAAGGAAAATGGCGGGCATCAAA‑AAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCttt     <  1:691290/238‑1 (MQ=255)
                                                                                                                                                                                                              ggAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAAAAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACCGCGC‑CCGCGAACCGGCGAAGGGGGATTGAAAAGCCCC‑GATACCAAACGCCCGCACAATGGTGGCCTGGAAGGTGGGGACAATTAAAACCGGTCCGGCCTTGcc  >  2:928279/1‑237 (MQ=255)
                                                                                                                                                                                                               gAAAGGAGAATGGCGGGCATCATA‑AAGCCAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCa                                                             <  1:1264219/177‑1 (MQ=255)
                                                                                                                                                                                                                 aaGGAGAATGGCGGGCATCATA‑AAGAAAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACa                                                                                                                                                                                >  1:1149615/1‑63 (MQ=255)
                                                                                                                                                                                                                         tGGCGGGCATCATA‑AAAACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAAAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAAAAGCGC‑CAGCCAACCGGCGAAAGGGAATTGATACCCCCC‑GAAAACAAACGCCCGca                                                   >  2:1264219/1‑177 (MQ=255)
                                                                                                                                                                                                                          ggCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGaa                                                                                                                                                                                  >  2:1290816/1‑52 (MQ=255)
                                                                                                                                                                                                                          ggCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGaa                                                                                                                                                                                  <  1:1290816/52‑1 (MQ=255)
                                                                                                                                                                                                                              ggCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCg                                                                                                                       <  2:1469851/106‑1 (MQ=255)
                                                                                                                                                                                                                              ggCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCg                                                                                                                       >  1:1469851/1‑106 (MQ=255)
                                                                                                                                                                                                                                     tA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGc                                                        >  1:32332/1‑160 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                         
CCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTT‑AATCCACGTCAGGTTTT‑GCAGCCATACCGGCATGTTTTCCAC‑CGGAGAAACG‑TAACCGGAAAGGAGAATGGCGGGCATCATA‑AAGACAAACACGCCGATAAACGCCTGCTGTTGTGTTGAACAGAGTGATGAAATCAACAGACCGAATCCCACCAGCGAT‑AAACCATAAATCACCATCGTAAAGTAGAACAGCGC‑CAGCGATCCGGCGAAGGGGATTTGATACGCCCA‑GATACCAATCGCCAGCACAATGGTGGCCTGGAAGGTGGCGACAATTAACGCCGGTACGGCTTTGCC  >  NC_000913/825010‑825448

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