Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,9590AC27.6% 66.2 / 11.9 65V133G (GTA→GGA) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/15);  new base C (0/16);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 3.92e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGCATCGGCAT  >  NC_000913/989710‑990208
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                
cccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑c                                                                                                                                                                                                                                                                <  2:244409/251‑1 (MQ=255)
 ccATCGGCGTGAGGGAACGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑cc                                                                                                                                                                                                                                                               <  1:178808/251‑1 (MQ=255)
   atcggctggAGGGAAAGCTGGCGCGCCAGATCGTAGTTGCGGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTGCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGTCTAATTTGCTTACCAGCGTGTCGAGGGCTTGCGCAACGGCGGCCTGCGGTG‑AGCT‑ATATCGATGTG‑CCAt                                                                                                                                                                                                                                                             <  2:511896/243‑1 (MQ=255)
     cGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCAtt                                                                                                                                                                                                                                                            <  2:56198/250‑1 (MQ=255)
      ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTTGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑c                                                                                                                                                                                                                                                         <  1:64983/251‑1 (MQ=255)
         gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCCGCGTGTCGGGGGCTTGCGCAACGGCGGGCCGCGGGG‑AGCG‑ATAGCGATGTA‑ACCTTT‑CACt                                                                                                                                                                                                                                                      >  1:352193/1‑251 (MQ=255)
           gAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTg                                                                                                                                                                                                                                                     <  2:466263/250‑1 (MQ=255)
           gAGGGAAAGCCGACGCCCCGGCTCGTAGTGTCTGTTGCCCACAAACCAGGAGTCCTGTTGCAGACGCTTAACGATGGGTTCTTGGACTTCGCGAGAAAAACGGCGACGGGGGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑TCCG‑AGAAGTCGTCTAATTGGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGGT‑CCATTT‑CACTg                                                                                                                                                                                                                                                     <  2:321203/250‑1 (MQ=255)
                  aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGc                                                                                                                                                                                                                                             <  2:234079/251‑1 (MQ=255)
                                        tGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAACCACCG‑‑GAAATCGGCTAATTTGCTTTTCAACGTGGCGGGGGCTTGCGCCACGTCGGGCCGCGGGG‑AGGG‑GTAGCGGTGTA‑CCCTTT‑CCATGGTAACCGCCAGC‑AA‑AGG‑AAACTCCCAGAgg                                                                                                                                                                                                                    >  1:482125/1‑251 (MQ=255)
                                                cccACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTTCCGCGGGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGG‑GGCG‑GTAGCGATGTA‑CCATTT‑CACTGATAACCGCCCGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTc                                                                                                                                                                                                            >  1:233883/1‑251 (MQ=255)
                                                     aaaCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAATCGACTTATTTTCTTTCCCGCGTGTCGAGGGCTTGCGCAACCTCGGCCTGCGGGG‑GGCG‑ATAGCGATGTA‑TCATTT‑TCCTGATAACAGCCAAC‑AA‑AGG‑AAATTCCCCGAGGATGTCTTCACGcc                                                                                                                                                                                                       >  1:570064/1‑251 (MQ=255)
                                                         cAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAA‑GCCG‑CGAAGTCGACTAATTTGCTTTCCCGCGTGTCGAGGGCTTGCGCAACGGCGGCCTGCGGGG‑AGCG‑ATAGCGATGTA‑CCCTTT‑CACTGATAACCGCCCGC‑AA‑AAG‑AACTTCCCAGAGGATGACTTCACGCCACgg                                                                                                                                                                                                   >  1:459989/1‑251 (MQ=255)
                                                          aGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑GGAAAACGCCAAATTTGATTTCCAGCGTGTCCAGGGGTTTCGCAACAACGGCCTCCGGTG‑TGCG‑ATAAGGGTGTT‑CCCTTT‑TCCTGATAAACCCCAGA‑AA‑AGG‑AACCTcccctagggtgccttcacccccacgg                                                                                                                                                                                                  >  1:86019/1‑228 (MQ=255)
                                                               tCCTGTTGCAGCCGTTTAACTATGGTTTCTTGTATTCCGCGAAAAAACCGGCGACGGGTCCCAAAATCCAACAG‑ATGGAAGGGCGACATATCAAGCCGGGCGG‑TTAAC‑CCCG‑AGAAGTCGACTATTTTCCTTCCCAGCTTGTCGAGGGCTTGCCCAACGGCGCCCTGCGGTG‑GGCG‑TTAGCGATTTC‑CAATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGc                                                                                                                                                                                              <  1:391868/250‑1 (MQ=255)
                                                                 cggtgggaGACGCTTACCGTGGGTTTCTGGTTCTCCGCGAGAAAAAGGGCGACGGGTGCAAAAATCCATCAG‑GTGAAAGCGCGACATCTAAAGCCCGGCGG‑TTAAC‑GCGG‑GGAATACGCATATTTGGCTTTCCAGCGTGTCGGGGGTTGGCTCAACGCCGGCCGGCGGTG‑AGCG‑GTGGCTATGTA‑CCATTT‑CCCGGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGc                                                                                                                                                                                              <  2:357063/241‑1 (MQ=255)
                                                                             cTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑CGAAGGCGACTAATTTGCTTTCCCCCGGGGTGAGGGCTTGCGCCACGTCGGCCTTCGGTG‑GGCG‑GTAGCGATGTG‑CCATTT‑TACTGATAACCGCCCGC‑AACAGG‑A‑CGTCCCTGAGGGAGAATTCACGCCCCCGGCCGCTTaaaaaacaaaccc                                                                                                                                                                               >  2:100001/1‑241 (MQ=255)
                                                                             cTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAA‑GCCG‑AGAAGGCGACTAATTTGGTTTCCAGCGTGTCGAGGGGTTTCGCACCGTCGGCCTGCGGGG‑GGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑CA‑AAG‑AACTTCCCCAAGGATGACTTTAAGCCACGGGGCCCTTAAACGAAAATcc                                                                                                                                                                                >  1:474844/1‑250 (MQ=255)
                                                                                   gttGGTTTTTGGTACTTGGGGAGAAAACCGGCGAGGGGGGCAAAAACCCA‑CAGGATGAAGGCGCGATATACCAGGCCCGGCGG‑TAAAC‑GCGG‑AGAAGTCGACTATTTGGCTTCCCAGGGGGTCGGGGGTTTGCGCAACGGCGGCCTGCGGGG‑GGCA‑ATAGCGATGTC‑CCATTT‑CCCTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGc                                                                                                                                                                         <  2:221526/249‑1 (MQ=255)
                                                                                   gATGTTTTTTTGTTCTTCGAGAAAAAACCGGCGACGGGTGCAAAAACACATAAG‑AGGGAAGCGAGAAATATAAAGACCGGGGG‑TTAAC‑GCGG‑AGAAGTGGACAAATTTGCTTTCCAGGGGGTGGGGGGTGTGCGCAACGTCGGCCGGGGGTG‑GGCG‑ATAGAGATGTC‑CCATTT‑AACGGAAAACCGCCAGC‑AA‑AGG‑AACGTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGc                                                                                                                                                                         <  1:573128/251‑1 (MQ=255)
                                                                                   gATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCGCGGAAG‑CGCCTAATTTGCTTTTCAGCGGGTCGGGGGGTTGCGCAACCTCGGCCCGCTGGG‑AGCG‑GTGTCGATGAAACAATTT‑CAC‑GATAACAACCAGC‑AA‑AGA‑AACTTCCCAGAAGATGACTTTACGCCACGCGCCCCTTTAAcacaaatacagtcttc                                                                                                                                                                         >  2:163789/1‑236 (MQ=255)
                                                                                      ggTTTCTGGTACTTCGCGAGAAAACCGGCGCCGGGTGCCAAAACCCATCAG‑ATGAAAGCGCGAAATATCAAGCCCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCGGGGTGTCGAGGGTTTGCGAAACGTCGGCCTGGGGTG‑AGCG‑ATAGCGATTTC‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CAt                                                                                                                                                                     <  1:53571/251‑1 (MQ=255)
                                                                                           cTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAA‑GCCC‑AGAAGTCGACCAATTTTCTTTCCAGCGGGTCGCGGGCTTGGGCCACGTCGGCCTGAGGGG‑AGCG‑ATAGCGCTGTT‑CCATTT‑CCCTGATAACCGCCCGC‑CA‑AAG‑AACATCCCAGAGGGTGACTTCACGCCCCGGGCCCCCTAAACGAAAATCCAGCTTGc‑cgatcggg                                                                                                                                                                >  1:382195/1‑244 (MQ=255)
                                                                                                aCTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTg‑g                                                                                                                                                           <  1:51295/250‑1 (MQ=255)
                                                                                                  ttCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTTAC‑GCCG‑AGAAGTCGACTAATTTTCTTTCCAGCGGGTCGAGGGATTGCGCAACCTCGGCCTGCGGTG‑GGGG‑ATAGGGATGTA‑CCCTTT‑CCCTGATAACCCCCCGC‑AC‑AGG‑AACTTCCCAGAGGGTCACTTCACGCCACGGGCCGCTTACAAGAATATCCCACTTTC‑AATTATCGTTG‑Gac                                                                                                                                                         >  1:320798/1‑250 (MQ=255)
                                                                                                            aaaCGGCGCGGGGTGCAAAAACCCATAAG‑TGGGAAGCGCGACATATCAAGACCGGCGG‑TTAGC‑CCCG‑AGAAGTCACCTATTTGGCTTTCCAGGGGGGCGGGGGCTGGCGAAGCGTCGGCC‑GGGGTGGAGCG‑ATAGCGATGTA‑CAATTT‑CACTAATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTg                                                                                                                                              <  1:201399/251‑1 (MQ=255)
                                                                                                               cGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑GTAAC‑CCCG‑AGAAGGAGACTAATTTGTTTTTCAGCGCGTCGCGGGCGTGCGCCACGACGGCCTGCCGGG‑AGCG‑CCAGCGGGGTA‑CCCATT‑CACCGAAAACCACCAGC‑AA‑AGG‑AAATTCCCCGAGGATGACTTCACCCCCCCGGGCGCGTAAACGAAAATCCACCGTGC‑CATTTTCGCTT‑GACACGGTGAgttgttt                                                                                                                                           >  2:488699/1‑250 (MQ=255)
                                                                                                                ggCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGGGTTAA‑‑GCCG‑AGAAGTAGGCGAATATGCTTTCCAGCGGGTCGGGGGGTTGGGCCACGTCGGCCCGAGGGG‑GGGG‑GTATAGGTGTA‑CACTTT‑TACTCATAAACCCCAGC‑AA‑AGGGAATTTCCAAG‑GGGGGACTTCCCGCCCCGGGCCGCCTAAAAGAGTAACCAGCCNGCGCATT‑TCGTTG‑GACCCGGGGACTTGTTcc                                                                                                                                          >  2:277208/1‑251 (MQ=255)
                                                                                                                    aCGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCCGCGGTG‑AGCG‑ATTGCGGTGTA‑CCCTTT‑CACTGATAACCCCCCGC‑AA‑AGG‑AACTTCCCAGAGGGTGACTTCACGCCACGGGCCCCCTAAACGAAAATCCCGCTTGC‑CCATATCCTTG‑GAAACGGGGACTTGTTTCgggg                                                                                                                                      >  1:339387/1‑250 (MQ=255)
                                                                                                                      gggTGCCAAAATCCATCAG‑AGGGAAGCGCGACATACCAAGCCCGGCGG‑TTAAC‑GCCG‑AAAAGTCGACTAATTGGCTTTCCGGCGTGTCGGGGGCTTGCGAAACGTCGGCCGGGGGTG‑AGCG‑ATAGCGAGGTC‑CCATTT‑CCCTGAAAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTa                                                                                                                                    <  1:463046/251‑1 (MQ=255)
                                                                                                                       ggTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑ATAAA‑GCCC‑AGAAGGCGGCTAATTTTCTTTCCAACGTGTCGGGGGCTTTCGCCAAGGCGGCCTGTGGGG‑AGCG‑ATAGCGATGGA‑CCATTT‑CACTGCTAACCCCCCGCCAA‑AGA‑AACTCCCCAAAGG‑TGACTTTCCACCACCGGCCGCCTAAAACAAAATTCACCTTTC‑CATTATCGTTG‑GAAACGGTGGCTTGTTCCggggta                                                                                                                                    >  2:310967/1‑247 (MQ=255)
                                                                                                                           ccAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑CCCC‑AGAAGACGGCTAATTTGGTTTCCAGCGCGTCGAGGGGTTGCGCAACGCCGGCCTGCGGGG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACCGATAACCGCCCGC‑AA‑AGG‑AACTTCCCCGCGGATGACTTAACCCCACGGGCCGCTTAAACGAATATCCACCTTGC‑CATTATCGTTGAGACACGG‑GCATTGTTCCGGGGTaaac                                                                                                                                 >  2:90045/1‑248 (MQ=255)
                                                                                                                            caaaacccacaaG‑ATGGAAGCGCGACATATCAGGCCCGGGGG‑TTAAC‑GCCG‑AGACGTCGACTATTTTGCTTCCCGGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTC‑CCATTT‑CACGAATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCgg                                                                                                                              <  1:100635/240‑1 (MQ=255)
                                                                                                                              aaaTCCACCAG‑GTGGAGCCGCGACATATCAAGACCGGCGT‑TTAAC‑GCCG‑AGAAGTCGTCTAATTGGCTTTCCAGCGTGTCGGGGGCTTGCGCACCGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTC‑CCATTT‑CCCTGATATCCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGGTCCGGGTTAAAGCGGaa                                                                                                                            <  1:444579/251‑1 (MQ=255)
                                                                                                                                    aTCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGGAGTCGACTAATTTTCTTTTCCGCGTGTCGGGGGCTTGCGCAACGTCGGCCTGCGGGG‑AGCG‑ATAGCGCTGTA‑CCATTT‑CCCTGATAACCGCCCGC‑AA‑AAG‑AACTTCCCAAAGGATGACTTCACGCCCCCGGCCGCTTAAACAAATATCCAACTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGGTAAAGCGGAACTcgcg                                                                                                                      >  1:377425/1‑251 (MQ=255)
                                                                                                                                     gagg‑ATGGAACCTCGACATTGCAAGACCGGCGG‑GAACC‑GCAG‑AGAAGTCCACTAATTTCCTTTCAGGGGTGTCGGGGTCTTGCGAAACTTCGGCCGGCGGTG‑GGCA‑ATAGCGATGTC‑CCTTTT‑CACTGAAAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGt                                                                                                                     <  2:380428/248‑1 (MQ=255)
                                                                                                                                      ccagaggaaaGCGCGCCATA‑CAAGGCCGGCGGGTTACC‑GCCG‑AGAAGCCGACTATTTTGCTTTCCAGCGTGTGGGGGGCTTGCGAAACGGCGGCCTGGGGGG‑AGCG‑ATAGCGATGTC‑CCATTT‑CCCTGATACCCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGt                                                                                                                     <  1:198928/243‑1 (MQ=255)
                                                                                                                                          aTGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGGAGGAGAACAATTTGCTTTCCACCGGGTTGAGGGCTTGGGCCACGTCGGCCCTCGGTG‑GGCGGTTAG‑GATGGA‑CCATTT‑TTCTGATAACCGCCCGC‑AA‑AGGGAACT‑CCCCGAGGGTGACATCACGCCACGGGCCGCTTAAAAGCAAATTCCGACTGC‑CATTTTCGTTT‑GGCACGGTGGCTTGTTTCGGGTTTAAggcgaacatgcctaaac                                                                                                                 >  2:100984/1‑235 (MQ=255)
                                                                                                                                                            caagccgggggttgaaC‑GCGG‑AGAAGTCGTTTAATTGTCTTCCCAGGGTGTGGGGGTCTGGCGCACCGTCGGCCTGCGGGG‑AGCG‑AGAGCGATTTT‑CCATTT‑CCCGGATACCCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTGAAGATAGTCGGcc                                                                                                <  2:572042/237‑1 (MQ=255)
                                                                                                                                                               ccgggcgg‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGTTTTCACGCGTGTCGAGGGCTTGCTAAACGTCGGCCGGCGGTG‑AGCG‑ATAGCGAGGTC‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaa                                                                                            <  1:45081/248‑1 (MQ=255)
                                                                                                                                                                    gggcgtttaccgcggaaaattCGACTAATTTGCTTCCCGGCTTGTGGAGGGTTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCAATTTC‑CAATTC‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaag                                                                                          <  1:193635/231‑1 (MQ=255)
                                                                                                                                                                         tAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCGGCTTGTCGGGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTC‑CCATTC‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggc                                                                                    <  1:494475/250‑1 (MQ=255)
                                                                                                                                                                                cG‑AGAAGTCGACTAATTTCCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGG‑GGCG‑ATAGCGATGTC‑CCATTT‑CCCTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcagggcag                                                                             <  1:454578/251‑1 (MQ=255)
                                                                                                                                                                                       gtcgcatatttcgttttCCGGCGTCTGGAGGGCTTGCGCAACGGCGCCCTGGGGTG‑AGGG‑ATAGCAATGT‑‑CCATTTTCACTGAAACCCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGcc                                                                       <  1:124558/237‑1 (MQ=255)
                                                                                                                                                                                           aCTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCACCGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATACCCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTg                                                                   <  1:310967/251‑1 (MQ=255)
                                                                                                                                                                                                  tgctctccggggggccGAGGGTTTGCCAAACGTCGGCCTGGGGGG‑GGGG‑ATGGCGATGTC‑CCATTC‑CACTGAAAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtc                                                             <  1:100984/235‑1 (MQ=255)
                                                                                                                                                                                                     tttCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑ACCTTCCCAGAGGATGACTTCACGCCCCGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTACAGAAAGGCAGGGCAGACAGCCACTGATATtcatc                                                          >  2:168581/1‑250 (MQ=255)
                                                                                                                                                                                                       tCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCatcatc                                                       >  1:363880/1‑251 (MQ=255)
                                                                                                                                                                                                            cGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACCCGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCAACTGATATTCATCATCCTGcg                                                  >  1:428411/1‑250 (MQ=255)
                                                                                                                                                                                                             gtgccGAGGCCTTGCGCCACGTCGGCTTGGGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CGCTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                                                 <  1:345139/247‑1 (MQ=255)
                                                                                                                                                                                                               gtCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGgc                                               >  1:100493/1‑251 (MQ=255)
                                                                                                                                                                                                                  taGGGCTTGCGCAACGTCTGACTGCGGTG‑AGCG‑AAAGAGATGTA‑CAATTT‑CACTGATACCCGCCAGC‑AA‑ATG‑AACTTCACAGAGGATGCCTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGA‑CATTATCGTTG‑GACACGGAGACTTGTTCAGGGTTAACGCGGAACTCGCGTAAACAGTTCAGATAGTCGGCCTTAAGGAAAGGCAGGGCAGAAAGCCCCGAATATTCATCATCCTTCAGgcgc                                             >  2:404990/2‑250 (MQ=255)
                                                                                                                                                                                                                  gAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGAAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGgcgc                                             >  2:176665/1‑250 (MQ=255)
                                                                                                                                                                                                                  gAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCACAGAGGATGACTTCACGCCACGGGCAGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCAACTGATATTCATCATCCTGCAGgcgc                                             >  1:6664/1‑250 (MQ=255)
                                                                                                                                                                                                                         tGCGCAACGTAGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgctg                                     >  1:128710/1‑251 (MQ=255)
                                                                                                                                                                                                                            gcAACGTCGGCCTGCGGTG‑CGGC‑AGAGCGATGTA‑CCATTT‑CAGTAATAACCGCCAGC‑AA‑AGG‑AACGTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCAGCATCCTGCAGGCGCAGGTGCTGCAt                                  <  1:90302/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGCCTGAGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTg                            >  2:442271/1‑250 (MQ=255)
                                                                                                                                                                                                                                     gCCTGCGGTG‑AGCT‑ATAGCAATGTC‑CCTTTT‑CCCTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAAc                         <  1:277208/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggggtgagcgaaggcattGTA‑CCATTC‑CACTAAAAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGt                     <  1:100001/234‑1 (MQ=255)
                                                                                                                                                                                                                                           ggTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑CACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTc                   >  1:206399/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gg‑AGCG‑ATAGCGATGTA‑CCATTT‑ACCTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCAc                 <  1:341660/250‑1 (MQ=255)
                                                                                                                                                                                                                                              g‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACg                >  2:464494/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      aagcgttGTC‑CCATTG‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGc          <  1:488697/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        gcgatGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGCAt        >  1:380125/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            tGTA‑CCATTT‑CGCTGAAAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGCATCgg     <  2:459989/250‑1 (MQ=255)
                                                                                                                                                                                                                                                               a‑CCATTG‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGCTGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTCGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCCGTTAAGATAGACGGCCATAAAGAAAGGCAGGGCAGAAAGCCACTGATATACATAATCATGCAGGCGCAGGGGCTGCATCGCCTGACCCTGTTCCCGACTAGCATGGGCAt  >  2:23515/1‑250 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                
CCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAG‑ATGGAAGCGCGACATATCAAGACCGGCGG‑TTAAC‑GCCG‑AGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTG‑AGCG‑ATAGCGATGTA‑CCATTT‑CACTGATAACCGCCAGC‑AA‑AGG‑AACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGC‑CATTATCGTTG‑GACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACCTGTTCACGAATAGCATCGGCAT  >  NC_000913/989710‑990208

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: ‑