Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 989,945 A→G 21.4% S138P (TCA→CCA)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,9450AG21.4% 80.8 / 2.8 48S138P (TCA→CCA) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (14/21);  new base G (3/8);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCA  >  NC_000913/989711‑990175
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
ccgtCGGCGCGAGGGAAAGCCGACGCGCCAGTTCGTAGTTGCTTGTGCCCCGAAACCAGGATTCCTGTTGCGGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAACGCGGCGACGGGTGCCAATATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                       <  2:1051987/236‑1 (MQ=255)
ccacgggCGGGAGGGAAAGCCGACGCGCCAGAACGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGGCGACTAATTTGCGGTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                       <  2:336952/234‑1 (MQ=255)
ccATCGGCGTGACGGGAAGCCGGCGCGCCAGAGCGCAGTTGCTGGTGCTCACAAACCAGGATTCCTGTCGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                       <  2:168200/239‑1 (MQ=255)
     ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGGCGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑CCCCGGAAGACGGCTAATTTGCCTTTCCGCGGGGTGGGGGGTTG‑CGCCACGTCGGCCCCCGGGGagcgataacg                                                                                                                                                                                                                                  >  1:438034/1‑239 (MQ=255)
     ggCGTGAAGGAAAGCCGACTCGCCAGAACGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                  <  1:616190/239‑1 (MQ=255)
     ggCGGGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCCAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTAA‑GACAGCAAGGCGCCGAACTTGTTTTTCATAGCGTGGAGGGGTTG‑CGCCACGGCGCGCGGCGGTGggcgctgacg                                                                                                                                                                                                                                  >  2:1176973/1‑233 (MQ=255)
       cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcgatagcg‑a                                                                                                                                                                                                                                <  1:976188/238‑1 (MQ=255)
        gggAGGGAAAGACGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATg                                                                                                                                                                                                                              <  1:673649/237‑1 (MQ=255)
          gAGGGAAAGCCGACGCGACAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAAGCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGtt                                                                                                                                                                                                                            <  1:496940/239‑2 (MQ=255)
           aGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGAGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTa                                                                                                                                                                                                                            <  1:75419/238‑1 (MQ=255)
                    cGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTACAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCAc                                                                                                                                                                                                                 <  2:831923/238‑1 (MQ=255)
                             agagCGTAGTTGCGGGTGCCCACAAACCAGGATTCCTGTGGCGAGCGCTTAACATTGGTTTCTTGTACTTC‑GGGAGAAAAACGGGGACGGGTGCCAAAATCC‑ATAAGAGGGAAGCGCGACATATAAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGTTTG‑CGAAACGTCGGCCGGCGGGGAGCGATAGCG‑ATGTC‑CCA‑TTTCACTGATAACCg                                                                                                                                                                                                        <  2:95167/238‑1 (MQ=255)
                                           ggTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAA‑ACCGAGAAGGCGACTAATTTGCTTTCCAGCGGGTCGGGGGGTTG‑CGTGACGGCGGCCCGCCGTGGTCGGTAGGG‑GCGCA‑CCA‑TTTTCATGATAAACGCCAGCACGGCAACtt                                                                                                                                                                                         >  2:509553/1‑239 (MQ=255)
                                                      aCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcgat                                                                                                                                                                                                                                      <  1:1129982/186‑1 (MQ=255)
                                                           gATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAA‑CCCCAGAAAACGACTAATTTGCTTTTCAGCGCGGCCGGGGCTTG‑CGCAACGTCGGCCCGCGGTAAGCTGTAGCG‑GTGTACCCA‑TTCCAC‑GATGACCTCCTCCAACAGGACTCACCAGAGGAGGACtt                                                                                                                                                                          >  2:406889/1‑238 (MQ=255)
                                                                cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGGCGAGGGCTTG‑CGCAACGTCGGCCCGCGGTGAGCGATAGCG‑ATGGA‑CCA‑TTTCACTGATAACCCCCAGCAAAGGAACTTCCCAGAGGGAGACTTCACGc                                                                                                                                                                     >  1:934156/1‑238 (MQ=255)
                                                                  gTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGCCTAATTTTCTTTTCAGCGTGTCGGGGGGTTG‑CGCAAAGTCGGCCTGGGGGGGGCGATAGCG‑GTGTA‑CCA‑TTTCACTGATAACCGCCACCAAAGGAACTTCCCAAGGGATGACTTCCCGcccc                                                                                                                                                                  >  1:90289/1‑237 (MQ=255)
                                                                     gCAGACGCTAAACGATGGTTTCTGGTACTGC‑GCGAGAAAACCGGGGACGGGGGCAAAAATCC‑ATCAGAGGGAGGCGCGACATATCAGGACGGGGGGTAAAC‑GCGGAGAGGTCGTCTAATTTGCTTTCCGGCGTGTCGAGGGCTTG‑CCAACCGCCGGCCGGAGGGGGGAGATAGCG‑AGGTA‑CAT‑TTTCACGGAGAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACggg                                                                                                                                                               <  2:438034/239‑1 (MQ=255)
                                                                       aGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCTGTGagcgata                                                                                                                                                                                                                                     <  1:207612/170‑1 (MQ=255)
                                                                       aGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAA‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCTGTGagcgata                                                                                                                                                                                                                                     >  2:207612/1‑170 (MQ=255)
                                                                             ttAAGGTTGGTTTCTGGTACTAG‑GCGAAAAAAAGGGGGA‑GGGTGGCAAAAACCAATCAGAGGGAGGGGCTAAATAACAAGACGGGGGGTGACC‑CCGTAGAAGTCGACTATGTTGATTTCCGGGGGGCGGGGGGTTGG‑CGCAACGTCGGCCTGCGGTGGAGGTTAGCG‑TGGTA‑CCA‑TTTCACTAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTa                                                                                                                                                        <  2:341389/238‑1 (MQ=255)
                                                                             taaaggATGGTTTCTTG‑ACATCGGCGAGAAAAACGGGGCGGGGTGCAAAAATCC‑AACAGGTGGAAGAGCGAAATATCAAGACCGGGGGTTAAC‑GCGGAAAAGGCGACTAATTTGTTTCCCAGCGTGTGGAGGGTTTG‑CGAAACGCCGGCCGGCGGTGGGCGATAGCG‑AGGTA‑CCA‑TTTCCCGAAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTa                                                                                                                                                        <  2:383272/233‑1 (MQ=255)
                                                                                     ggTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGGCGACTCAATTTCTTTCCAGAGGGGCGGGGGGCTG‑CGCAAAGCCGGCCTGCGGGGAGCGGAAGCG‑AGGTA‑ACA‑ATTTTCTCATAAACCCCAGCAAAGGGACTTCCCACAGGAGGGATTTACTCCACGGGCCCCTTTAAAGaaaa                                                                                                                                                >  2:981972/1‑236 (MQ=255)
                                                                                       tttCTTGTCCTCC‑GGAAAAAAACCGGCGCGGGGGGCAAAAACCC‑ATCAGATGGAAGCGGGAAATACAACGACCGGCGTTTACC‑GCCGAGAAGCCGACTAATTTGCTTCCCAGGGTGTCGAGGGTTTG‑CGAACCGCCGGCCTGCGGTGGGCGTTAGCG‑ATTTA‑CCA‑TTTCCCGGAAACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATc                                                                                                                                              <  2:1059599/238‑1 (MQ=255)
                                                                                         tgttgcacggc‑gcgAAAAAACCGGCGAGGGGTGCAAAAACCC‑ATAAGATGGAAGCGCGACATACCAAGACGGGGGTTTAAC‑GCCGAGAAGGCGACTAATTTGCTTCCCGGCGTGTCGGGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAg                                                                                                                                           <  2:1030287/229‑1 (MQ=255)
                                                                                               aCTTC‑GCGAGAAAAACGGCGACGGGTCCCAAAACCC‑ATCAGATGGAAGCGCGACATATCAAGCCCGGCGTTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTGG‑CGAAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGc                                                                                                                                      <  1:838957/238‑1 (MQ=255)
                                                                                                   c‑gcgAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGGCGACTAATTTGCTTTTCAGCGGGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGAGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGATTGCCATTa                                                                                                                                 >  1:858265/1‑239 (MQ=255)
                                                                                                   c‑gcgAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCCAGAAAGGGGATAAATTGCTTTCCAGAGGGTCGGGGGGTTGCCCCAACG‑CGGCCCGCCGGGGGGCATAGCG‑GTGGA‑CCATTTTCCCTATCAACCCCCCGAAAAGAACCTTCCC‑GGGGATTGATTCCCGCCCCCGGCCCCCTAAACGAATATCCaccctggccttt                                                                                                                                 >  2:129666/1‑228 (MQ=255)
                                                                                                   c‑gcgAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGCTAAC‑GCCCAGAAGTCGACTAATTTGCTTTCCAGCGCGTCGGGGGCTTG‑CGCAAAGTCCGCCTGCCGTGAGCGATAGCG‑GTTGA‑CCA‑TTTCACTGATAACCGCCCGCAAAGGAACTTCCCCGAGGATGACTTCACACCCCCGGCCCCCTTAACGAATATCCCACCTTCCATTa                                                                                                                                 >  2:176576/1‑239 (MQ=255)
                                                                                                       gataaagagcgcggcgGGTCCACAAATCC‑ATGCTATGGAAGGTCGACATAGCAAGCCCGGTGGTTAAC‑GCAGGAAAGCCGTCTTCTTCTCTTTCCAGGGTGTGGGGGTTGTG‑GCCACCGTCGGGTCGGGGGGGCCGATAACT‑ATGTC‑CCA‑TTACACGGATATCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTGGCCATTATCg                                                                                                                              <  2:742607/229‑1 (MQ=255)
                                                                                                                 ggggcGGGTGAAAAAACCC‑ATAATAGGGAAGCGCGACATACAAAGCCTGGCGGTTAGC‑GCCGAGAGTCCGATTAGTGTGCTTTCCGGCGTGCCGAGGGGTGG‑CGCAACTCCGGCCGGGGGGGAGCGAT‑GCGAATGTC‑CCT‑TTTCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACg                                                                                                                     <  2:779208/234‑1 (MQ=255)
                                                                                                                                            ggAAGCGCGACATATCAAGACCGGCGGATAAA‑GCCGCGAAATCGAATAATTTTCTTTCCAGCGTGTCGGGGGCCTG‑CGCAACGGCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTTACTTATAACCGGCAGCAAAGGGACTTTCCAGAGGAAGACTTCAAGCCACGGCCCGACTAAACGAATATCCAACTTGCCCATATAGTTGGACCCTGTGGCTTTTTCCGGGTTAAAGCggc                                                                                            >  2:897926/1‑236 (MQ=255)
                                                                                                                                               aGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGGCGACTAATTTGCTTTCCAGCGGGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAAACGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAAAATCCAGCTTGCCCTTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTc                                                                                        >  1:1072813/1‑238 (MQ=255)
                                                                                                                                               aGCGCGACATATCAAGACCGGCGGTTAAA‑GCCCGGAAGTCGGCTAATTTTCTTTCCCGCGTGTCGAGGGCTTT‑CGCCACGTCGGCCTGCGGGGTGCGCTAGCG‑GTGTA‑CCA‑TTTCCCTGATAACCGCCCGCAAAGGAACTTCCCCGAGGAGGACTTTACCCCACCGGCCGCCTTAACGAAAACTCCGCCTGCCCTTATCTTTGTGCACGGGGACTTGTTCCAGGGTAAAGCgggaacc                                                                                        >  2:943205/1‑233 (MQ=255)
                                                                                                                                                     ccATATAAGGAGCGGTGGTAAAC‑GCGGAAAATGCGCTTAATGTGCTTCCCGGGGTTTGGGGGCCTGG‑CGCAGCGTCGGCCGGCGGTGGGCAATGGCG‑ATGTA‑CCA‑TTTCCCGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                  <  2:1087425/237‑1 (MQ=255)
                                                                                                                                                     caaagaacaagccGGGGGGTAACGGCCGAGATGCCGAATAATTTGCTTTCGAGGGGTTCGAGGTCTTC‑GGCAACGTCGGCCGGGGGTGAGCATTAGGT‑GTGTA‑CCT‑TTCCCCGGATAGCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                  <  2:576046/228‑1 (MQ=255)
                                                                                                                                                     aCATATAAATCCCGGCGGTTACG‑GCCGAGAAGTCGACTAATTTGCTTTCCGGGGTGTCGGGGGCTTG‑CGCAACGGCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                  <  1:34652/238‑1 (MQ=255)
                                                                                                                                                        tatCAAGCCGGGCGGTTACC‑GCCGAGAAGTCGCTTAATTTGTTTTCCGGCGTGTCGAGGGTTTG‑CAAAACGTCGGCCTGCGGTGAGCAATAGCA‑AGGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAg                                                                              <  2:1072813/239‑1 (MQ=255)
                                                                                                                                                                   gcggTAAAC‑GCCGAGAAGTCGACTAATTTGCGTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGTTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                    <  1:347951/238‑1 (MQ=255)
                                                                                                                                                                                    agtgtgatagtttGATGTCCAGGGTTTCGAGGGTTTG‑CGCAACGCCGGCCGGCGTGGAGCAAAAACA‑ATGTA‑CCA‑TTTCCCTAAAAGCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggc                                                   <  2:1060354/229‑1 (MQ=255)
                                                                                                                                                                                                tGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAg                                        >  1:490804/1‑238 (MQ=255)
                                                                                                                                                                                                ggCTTTCCCGCTTGTCGAGGGCTTG‑CCGACCAGCGGCGTGGGCAGCGAGATAGAG‑ATGTC‑CCA‑TTTCTATGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                                       <  1:329344/238‑1 (MQ=255)
                                                                                                                                                                                                   ttcccggCGTGTGGTGGGCTTG‑CCCAACGTCGCCCGCGGGTGAGCGATAGCA‑AGTTA‑CAT‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCa                                     <  1:689224/232‑1 (MQ=255)
                                                                                                                                                                                                   ttccaagggtgcCGAGGTTTTG‑CACAACGTCGGCCTGCGGTGGGCGATACCA‑TTTCA‑CCA‑TTGCACTTATAGCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCAc                                    <  1:1031769/231‑1 (MQ=255)
                                                                                                                                                                                                    ttccacgtggtCGAGGGTTTG‑CACAACGTCGGCCTGGGGTGAGCGATAGCA‑ATTAT‑CCA‑TTACACTTAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACt                                   <  1:553835/230‑1 (MQ=255)
                                                                                                                                                                                                    ttCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGcc                                      >  1:898767/1‑236 (MQ=255)
                                                                                                                                                                                                       cggcggttCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGat                                <  1:934470/232‑1 (MQ=255)
                                                                                                                                                                                                             gtCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAACG‑ATGTA‑CCA‑TTTCACTGCTAACCGCCAGCAAAGGCACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGCATATCCAGCATGCCATTATCGTTGGACACGGTGACTTGTTACGGGTTAAAGCGGAACTCGAGTAACCAGTTAAGATAGTCGGCCTTACAGACAGGCAGGGCGGAAAGGCACTGAACTtcat                          >  2:666170/1‑239 (MQ=255)
                                                                                                                                                                                                                     cTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGc                  >  1:1067453/1‑239 (MQ=255)
                                                                                                                                                                                                                      ttG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCa                 >  2:353335/1‑239 (MQ=255)
                                                                                                                                                                                                                        g‑cgcAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                >  1:904207/1‑238 (MQ=255)
                                                                                                                                                                                                                               cGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGCAAGCCACTGATATTCATCATCCTGTAGGCGCAgg         >  1:1061032/1‑239 (MQ=255)
                                                                                                                                                                                                                                      cggCGGTGAGCGATTGCT‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  1:1079509/237‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGGGGTGGGCGCTAGCG‑ATGTA‑CCT‑TTTCACTGCTACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  2:257671/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGGGGTGAGGGATAGCG‑ATGTC‑CAA‑TTTCGCGAATAACCGCTAGCAAAGGAACTTCCAAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  1:323069/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGCGGTGAGCATTAGCA‑ATATC‑CCA‑TTTCACTTATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  1:174618/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGCGGTGAACGATAGCG‑ATGTA‑CCA‑TTTCACTGATAGCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  2:858265/239‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCGGTGAGAGAGAGCG‑AGGTC‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa  <  1:589659/237‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
CCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCA  >  NC_000913/989711‑990175

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