Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 989,890 T→G 22.2% D156A (GAC→GCC)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,8900TG22.2% 84.7 / 3.0 56D156A (GAC→GCC) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (11/25);  new base G (9/4);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 2.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGGCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGT  >  NC_000913/989662‑990109
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                
tCCGGGCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAAt                                                                                                                                                                                                                          <  2:1172103/237‑1 (MQ=255)
   gggCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCt                                                                                                                                                                                                                     <  1:234480/239‑1 (MQ=255)
         aTTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTAA‑GCCGGGGAATCGGCTAATTTTCttttcc                                                                                                                                                                                                                >  1:1030287/1‑235 (MQ=255)
         aTTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGCTAAC‑GCCACGAAGTCGGCTAATTTTCTTTccc                                                                                                                                                                                                                >  1:59029/1‑237 (MQ=255)
           ttGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑g                                                                                                                                                                                                            <  2:158107/239‑1 (MQ=255)
                     gCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGc                                                                                                                                                                                                  <  1:91808/239‑1 (MQ=255)
                     gCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGc                                                                                                                                                                                                  <  1:813202/239‑1 (MQ=255)
                      ccTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGACAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGc                                                                                                                                                                                                  <  1:92137/238‑1 (MQ=255)
                                ttCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGAAATATCAAGACCGGCGGTTAAA‑GCCGAGAAAGCGGCTAATTTGGTTTTCAGCGGTGGCGAGGG‑TTA‑CGCAAccg                                                                                                                                                                                      >  1:496772/1‑237 (MQ=255)
                                ttCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGGTTCTGGTGCCCACAAACCAGGATTCCTGTTGCACACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCACGACCGGCGGTTCAA‑ACCCAGACATACACTACTTTTCTTTCCAGG‑GTGTCGCGGCCTTt‑cgccccgt                                                                                                                                                                                      >  2:166748/1‑234 (MQ=255)
                                           gtgtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTAA‑ACCGCCAAGTGGGAAAAATTGTTTTCCCAC‑GCGTCGGGGGCATg‑cggcaccgccgcccgtggg                                                                                                                                                                           >  2:838957/1‑222 (MQ=255)
                                            tgtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAATCGACTAATTTGCTTTCCCGC‑GTGTCGGGGGGTTG‑GGGAACGCCGGccccggggg                                                                                                                                                                          >  1:257671/1‑232 (MQ=255)
                                             gtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTg                                                                                                                                                                          <  2:774788/238‑1 (MQ=255)
                                              ttCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAAt                                                                                                                                                                                                                          <  1:637580/191‑1 (MQ=255)
                                                 ccgtCGGCGCGAGGGAAAGCCGACGCGCCAGTTCGTAGTTGCTTGTGCCCCGAAACCAGGATTCCTGTTGCGGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAACGCGGCGACGGGTGCCAATATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                     <  2:1051987/236‑1 (MQ=255)
                                                 ccacgggCGGGAGGGAAAGCCGACGCGCCAGAACGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGGCGACTAATTTGCGGTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                     <  2:336952/234‑1 (MQ=255)
                                                 ccATCGGCGTGACGGGAAGCCGGCGCGCCAGAGCGCAGTTGCTGGTGCTCACAAACCAGGATTCCTGTCGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                     <  2:168200/239‑1 (MQ=255)
                                                   aTCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAAA‑GCCCAGAAAGCGACAAATTTTCTTTTCAGC‑GTGGCGGGGGGTTG‑GGCCAAGTCgggcctcggggaggggta                                                                                                                                                                   >  1:383272/1‑223 (MQ=255)
                                                      ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGGCGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑CCCCGGAAGACGGCTAATTTGCCTTTCCGC‑GGGGTGGGGGGTTG‑CGCCACGTCGGCCCCCGGGGagcgataacg                                                                                                                                                                >  1:438034/1‑239 (MQ=255)
                                                      ggCGTGAAGGAAAGCCGACTCGCCAGAACGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                <  1:616190/239‑1 (MQ=255)
                                                      ggCGGGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCCAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTAA‑GACAGCAAGGCGCCGAACTTGTTTTTCATA‑GCGTGGAGGGGTTG‑CGCCACGGCGCGCGGCGGTGggcgctgacg                                                                                                                                                                >  2:1176973/1‑233 (MQ=255)
                                                        cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcgatagcg‑a                                                                                                                                                              <  1:976188/238‑1 (MQ=255)
                                                         gggAGGGAAAGACGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATg                                                                                                                                                            <  1:673649/237‑1 (MQ=255)
                                                           gAGGGAAAGCCGACGCGACAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAAGCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGtt                                                                                                                                                          <  1:496940/239‑2 (MQ=255)
                                                            aGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGA‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTa                                                                                                                                                          <  1:75419/238‑1 (MQ=255)
                                                                     cGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTACAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCAc                                                                                                                                               <  2:831923/238‑1 (MQ=255)
                                                                              agagCGTAGTTGCGGGTGCCCACAAACCAGGATTCCTGTGGCGAGCGCTTAACATTGGTTTCTTGTACTTC‑GGGAGAAAAACGGGGACGGGTGCCAAAATCC‑ATAAGAGGGAAGCGCGACATATAAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGTTTG‑CGAAACGTCGGCCGGCGGGGAGCGATAGCG‑ATGTC‑CCA‑TTTCACTGATAACCg                                                                                                                                      <  2:95167/238‑1 (MQ=255)
                                                                                            ggTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAA‑ACCGAGAAGGCGACTAATTTGCTTTCCAGC‑GGGTCGGGGGGTTG‑CGTGACGGCGGCCCGCCGTGGTCGGTAGGG‑GCGCA‑CCA‑TTTTCATGATAAACGCCAGCACGGCAACtt                                                                                                                       >  2:509553/1‑239 (MQ=255)
                                                                                                       aCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGagcgat                                                                                                                                                                    <  1:1129982/186‑1 (MQ=255)
                                                                                                            gATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAA‑CCCCAGAAAACGACTAATTTGCTTTTCAGC‑GCGGCCGGGGCTTG‑CGCAACGTCGGCCCGCGGTAAGCTGTAGCG‑GTGTACCCA‑TTCCAC‑GATGACCTCCTCCAACAGGACTCACCAGAGGAGGACtt                                                                                                        >  2:406889/1‑238 (MQ=255)
                                                                                                                 cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGGCGAGGGCTTG‑CGCAACGTCGGCCCGCGGTGAGCGATAGCG‑ATGGA‑CCA‑TTTCACTGATAACCCCCAGCAAAGGAACTTCCCAGAGGGAGACTTCACGc                                                                                                   >  1:934156/1‑238 (MQ=255)
                                                                                                                   gTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGCCTAATTTTCTTTTCAGC‑GTGTCGGGGGGTTG‑CGCAAAGTCGGCCTGGGGGGGGCGATAGCG‑GTGTA‑CCA‑TTTCACTGATAACCGCCACCAAAGGAACTTCCCAAGGGATGACTTCCCGcccc                                                                                                >  1:90289/1‑237 (MQ=255)
                                                                                                                      gCAGACGCTAAACGATGGTTTCTGGTACTGC‑GCGAGAAAACCGGGGACGGGGGCAAAAATCC‑ATCAGAGGGAGGCGCGACATATCAGGACGGGGGGTAAAC‑GCGGAGAGGTCGTCTAATTTGCTTTCCGGC‑GTGTCGAGGGCTTG‑CCAACCGCCGGCCGGAGGGGGGAGATAGCG‑AGGTA‑CAT‑TTTCACGGAGAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACggg                                                                                             <  2:438034/239‑1 (MQ=255)
                                                                                                                        aGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCTGTGagcgata                                                                                                                                                                   <  1:207612/170‑1 (MQ=255)
                                                                                                                        aGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAA‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCTGTGagcgata                                                                                                                                                                   >  2:207612/1‑170 (MQ=255)
                                                                                                                              ttAAGGTTGGTTTCTGGTACTAG‑GCGAAAAAAAGGGGGA‑GGGTGGCAAAAACCAATCAGAGGGAGGGGCTAAATAACAAGACGGGGGGTGACC‑CCGTAGAAGTCGACTATGTTGATTTCCGGG‑GGGCGGGGGGTTGG‑CGCAACGTCGGCCTGCGGTGGAGGTTAGCG‑TGGTA‑CCA‑TTTCACTAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTa                                                                                      <  2:341389/238‑1 (MQ=255)
                                                                                                                              taaaggATGGTTTCTTG‑ACATCGGCGAGAAAAACGGGGCGGGGTGCAAAAATCC‑AACAGGTGGAAGAGCGAAATATCAAGACCGGGGGTTAAC‑GCGGAAAAGGCGACTAATTTGTTTCCCAGC‑GTGTGGAGGGTTTG‑CGAAACGCCGGCCGGCGGTGGGCGATAGCG‑AGGTA‑CCA‑TTTCCCGAAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTa                                                                                      <  2:383272/233‑1 (MQ=255)
                                                                                                                                      ggTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGGCGACTCAATTTCTTTCCAGA‑GGGGCGGGGGGCTG‑CGCAAAGCCGGCCTGCGGGGAGCGGAAGCG‑AGGTA‑ACA‑ATTTTCTCATAAACCCCAGCAAAGGGACTTCCCACAGGAGGGATTTACTCCACGGGCCCCTTTAAAGaaaa                                                                              >  2:981972/1‑236 (MQ=255)
                                                                                                                                        tttCTTGTCCTCC‑GGAAAAAAACCGGCGCGGGGGGCAAAAACCC‑ATCAGATGGAAGCGGGAAATACAACGACCGGCGTTTACC‑GCCGAGAAGCCGACTAATTTGCTTCCCAGG‑GTGTCGAGGGTTTG‑CGAACCGCCGGCCTGCGGTGGGCGTTAGCG‑ATTTA‑CCA‑TTTCCCGGAAACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATc                                                                            <  2:1059599/238‑1 (MQ=255)
                                                                                                                                          tgttgcacggc‑gcgAAAAAACCGGCGAGGGGTGCAAAAACCC‑ATAAGATGGAAGCGCGACATACCAAGACGGGGGTTTAAC‑GCCGAGAAGGCGACTAATTTGCTTCCCGGC‑GTGTCGGGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAg                                                                         <  2:1030287/229‑1 (MQ=255)
                                                                                                                                                aCTTC‑GCGAGAAAAACGGCGACGGGTCCCAAAACCC‑ATCAGATGGAAGCGCGACATATCAAGCCCGGCGTTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTGG‑CGAAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGc                                                                    <  1:838957/238‑1 (MQ=255)
                                                                                                                                                    c‑gcgAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGGCGACTAATTTGCTTTTCAGC‑GGGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGAGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGATTGCCATTa                                                               >  1:858265/1‑239 (MQ=255)
                                                                                                                                                    c‑gcgAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCCAGAAAGGGGATAAATTGCTTTCCAGA‑GGGTCGGGGGGTTGCCCCAACG‑CGGCCCGCCGGGGGGCATAGCG‑GTGGA‑CCATTTTCCCTATCAACCCCCCGAAAAGAACCTTCCC‑GGGGATTGATTCCCGCCCCCGGCCCCCTAAACGAATATCCaccctggccttt                                                               >  2:129666/1‑228 (MQ=255)
                                                                                                                                                    c‑gcgAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGCTAAC‑GCCCAGAAGTCGACTAATTTGCTTTCCAGC‑GCGTCGGGGGCTTG‑CGCAAAGTCCGCCTGCCGTGAGCGATAGCG‑GTTGA‑CCA‑TTTCACTGATAACCGCCCGCAAAGGAACTTCCCCGAGGATGACTTCACACCCCCGGCCCCCTTAACGAATATCCCACCTTCCATTa                                                               >  2:176576/1‑239 (MQ=255)
                                                                                                                                                    c‑gcgAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGATAAC‑GCCGGCAAGGCGACGAATTTGATTTCCACC‑GCGGCGGGGGGTCG‑CGGAAAGCCGGCCTGCGGTGTGGGGTAGAG‑GGGTA‑GCC‑TTTCACCGATAACCCCCCCGCAGGGAAAGTACCAAAGGGTGACTTTACGCCACCGACCGGCTAAACGAaaataactacaggcaatg                                                               >  2:589659/1‑225 (MQ=255)
                                                                                                                                                       cgAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAA‑CCCGCGAAGACGACTACTTTGGTTTTCAGC‑GGGTCGGGGGCTTG‑CGCAACATCGGCCCTCGGGGAGCGATAGCG‑AAGGA‑CCA‑TTTTACTGCTAACCACCAGCCAAAGAACATCCCCCAGGAGGAATTCACGCCCACGGCAGCATAAAcaacgaatcatgtctcccacacg                                                             >  2:1191393/1‑217 (MQ=255)
                                                                                                                                                        gataaagagcgcggcgGGTCCACAAATCC‑ATGCTATGGAAGGTCGACATAGCAAGCCCGGTGGTTAAC‑GCAGGAAAGCCGTCTTCTTCTCTTTCCAGG‑GTGTGGGGGTTGTG‑GCCACCGTCGGGTCGGGGGGGCCGATAACT‑ATGTC‑CCA‑TTACACGGATATCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTGGCCATTATCg                                                            <  2:742607/229‑1 (MQ=255)
                                                                                                                                                                  ggggcGGGTGAAAAAACCC‑ATAATAGGGAAGCGCGACATACAAAGCCTGGCGGTTAGC‑GCCGAGAGTCCGATTAGTGTGCTTTCCGGC‑GTGCCGAGGGGTGG‑CGCAACTCCGGCCGGGGGGGAGCGAT‑GCGAATGTC‑CCT‑TTTCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACg                                                   <  2:779208/234‑1 (MQ=255)
                                                                                                                                                                                             ggAAGCGCGACATATCAAGACCGGCGGATAAA‑GCCGCGAAATCGAATAATTTTCTTTCCAGC‑GTGTCGGGGGCCTG‑CGCAACGGCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTTACTTATAACCGGCAGCAAAGGGACTTTCCAGAGGAAGACTTCAAGCCACGGCCCGACTAAACGAATATCCAACTTGCCCATATAGTTGGACCCTGTGGCTTTTTCCGGGTTAAAGCggc                          >  2:897926/1‑236 (MQ=255)
                                                                                                                                                                                                aGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGGCGACTAATTTGCTTTCCAGC‑GGGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAAACGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAAAATCCAGCTTGCCCTTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTc                      >  1:1072813/1‑238 (MQ=255)
                                                                                                                                                                                                aGCGCGACATATCAAGACCGGCGGTTAAA‑GCCCGGAAGTCGGCTAATTTTCTTTCCCGC‑GTGTCGAGGGCTTT‑CGCCACGTCGGCCTGCGGGGTGCGCTAGCG‑GTGTA‑CCA‑TTTCCCTGATAACCGCCCGCAAAGGAACTTCCCCGAGGAGGACTTTACCCCACCGGCCGCCTTAACGAAAACTCCGCCTGCCCTTATCTTTGTGCACGGGGACTTGTTCCAGGGTAAAGCgggaacc                      >  2:943205/1‑233 (MQ=255)
                                                                                                                                                                                                      ccATATAAGGAGCGGTGGTAAAC‑GCGGAAAATGCGCTTAATGTGCTTCCCGGG‑GTTTGGGGGCCTGG‑CGCAGCGTCGGCCGGCGGTGGGCAATGGCG‑ATGTA‑CCA‑TTTCCCGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                <  2:1087425/237‑1 (MQ=255)
                                                                                                                                                                                                      caaagaacaagccGGGGGGTAACGGCCGAGATGCCGAATAATTTGCTTTCGAGG‑GGTTCGAGGTCTTC‑GGCAACGTCGGCCGGGGGTGAGCATTAGGT‑GTGTA‑CCT‑TTCCCCGGATAGCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                <  2:576046/228‑1 (MQ=255)
                                                                                                                                                                                                      aCATATAAATCCCGGCGGTTACG‑GCCGAGAAGTCGACTAATTTGCTTTCCGGG‑GTGTCGGGGGCTTG‑CGCAACGGCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                <  1:34652/238‑1 (MQ=255)
                                                                                                                                                                                                         tatCAAGCCGGGCGGTTACC‑GCCGAGAAGTCGCTTAATTTGTTTTCCGGC‑GTGTCGAGGGTTTG‑CAAAACGTCGGCCTGCGGTGAGCAATAGCA‑AGGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAg            <  2:1072813/239‑1 (MQ=255)
                                                                                                                                                                                                                    gcggTAAAC‑GCCGAGAAGTCGACTAATTTGCGTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGTTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt  <  1:347951/238‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                
TCCGGGCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTC‑GCGAGAAAAACGGCGACGGGTGCCAAAATCC‑ATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAC‑GCCGAGAAGTCGACTAATTTGCTTTCCAGC‑GTGTCGAGGGCTTG‑CGCAACGTCGGCCTGCGGTGAGCGATAGCG‑ATGTA‑CCA‑TTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGT  >  NC_000913/989662‑990109

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