Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	989,943	T→G	26.5%	S138S (TCA→TCC) ‡	pncB ←	nicotinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	989,943	0	T	G	26.5%	96.7 / 3.3	70	S138S (TCA→TCC) ‡	pncB	nicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (19/28);  new base G (12/10);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATC  >  NC_000913/989707‑990177                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 gtTCCCATCGGCGTGAGGGAAAGCTGACGCGCCAGATCGGTGTTGCTGGTGCCCACAAACCACGATTCCTGGTGCAGACGCTTAACGATGGTTTCTTGTACGTCGCGAG‑AAAAAC‑GGCGACGTGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTg                                                                                                                                                                                                                                                 <  2:783044/238‑1 (MQ=255) gtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTg                                                                                                                                                                                                                                                 <  2:932134/238‑1 (MQ=255)     ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                            <  2:1402452/239‑1 (MQ=255)     ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                            <  1:935805/239‑1 (MQ=255)     ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                            <  1:828783/239‑1 (MQ=255)     ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                            <  1:532474/239‑1 (MQ=255)      cATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGagcgat                                                                                                                                                                                                                                           <  1:603600/239‑1 (MQ=255)          ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                       <  2:1226147/239‑1 (MQ=255)             gttaGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑Tg                                                                                                                                                                                                                                   <  2:1268585/236‑1 (MQ=255)              tGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CACCG‑AG‑AAGGCGAATAATTTGATTT‑TCAGCGGGTCG‑GGGGGTTGGGCAACGTCCGGCTGCGGGGGGCGGAAGGGA‑Tgg                                                                                                                                                                                                                                  >  1:1125113/1‑238 (MQ=255)                     aaGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCC‑AG‑TAGTCGACTAAATTTCTTT‑CCAGCCGGTCG‑GAGGCTTTCGCGACGTTCGCCCGCGGTGtgcggcagccgtggaaccttt                                                                                                                                                                                                                            >  2:813415/1‑222 (MQ=255)                                             tgtGGGGCCCAAAAACCAGGATCCCGGTTGCAGACGCTTAACGATGGTTTCTTGTCCTTCGCAAG‑AAAAAC‑GGCGGCGGGTGCAAAAATCCATC‑AGATGGAAGCGCGACATATCAAGCCCGGCGGTA‑AA‑CGCGG‑AG‑AAGTCGATTATTTTGCTTT‑CCAGCGTGTCG‑GGGGTTTGCGCAACGTCGGCCTGCGGGGAGTGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGaa                                                                                                                                                                                                   <  2:390814/237‑1 (MQ=255)                                             gCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACCAAATTGCTTT‑TCAACGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGaa                                                                                                                                                                                                   >  1:625618/1‑239 (MQ=255)                                              cTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑AGCCG‑AG‑AAGTCGAATAAATTTCTTT‑TCCGCGTGTCG‑AGGGCTTTGGCAACGGCGGCCTGCGGGGAGCGGTAGGGA‑GGTACCCTTTTACTGATAACCCCCCGGCAAGGa                                                                                                                                                                                                    >  2:564720/1‑237 (MQ=255)                                                ggTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑AGCCC‑AG‑GAGTCCACTAATTTTCTTT‑CCAGCGTGTCG‑CGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGA‑GGTACCATTTTACTGCTAACCGCCCGCAAAGGAACt                                                                                                                                                                                                 >  1:270813/1‑238 (MQ=255)                                                    cccACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑GG‑AAGTCGACTAATTTGCTTT‑CCAACGTGTCG‑CGGGGTTGCGCCACGGCGGCCCGCGGTGTGGGAGAGGGA‑TGTACCCTTTTACTGATAAACGCCAGCCCAGGAACttaca                                                                                                                                                                                             >  2:59910/1‑235 (MQ=255)                                                     ccACAAACCAGGATTCCGTTGGAAGACGCTTAACGTTGGTTTCTTGGCCTGCGCGAG‑AAAAAC‑GGGGACGGGGGCAAAAATCCATC‑AGATGGAAGCGAGACATACCAAGACGGGGGGTT‑AA‑CGCCG‑AA‑AGGTCGACTAATTTGCTTT‑CCGGGGTGTCG‑GGGGCTTGCGAAACGTCGGCCTGGGGTGAGCGATAGCGA‑TGTACCATTTCACGGATAACCGCCAGCAAAGGAACTT‑CCCa                                                                                                                                                                                           <  2:694009/238‑1 (MQ=255)                                                          aaCCAGAGTTCCGTTTGCAGACGTTAAACGATGTTTTTTTGTACTTGGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑GGGTGGAAGCGGGCCTTATCAAGACGGGCGGTT‑AA‑CGCGG‑AG‑AAGTCGATTATTTTGCTTT‑CCAGCGTGTCG‑GGGGTTTCGCCAACGTCGGCCGGCGGTGAGCGATAGCGA‑TGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAt                                                                                                                                                                                     <  2:534635/239‑1 (MQ=255)                                                           aCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑AGCCG‑GG‑AAGACGACTAATTTGCTTT‑CCCCCGGGTCG‑AGGGCTTGCGCCACGGCGGCCTGCGGTGGGCGGTAGCGA‑AGGAACATTTTACTGATAACCGCCAGCAAAGGAACTTGCCCAggggtg                                                                                                                                                                                     >  2:884123/1‑237 (MQ=255)                                                                 atgccTGGTGCAGTCTCTTAACGTTGGTTTTTTGTCCTCCGCGGG‑AAAACC‑GGCGGTGGGTGCCAAAACCCCAC‑AGAGGAAGGCGCGAAATATCAAGCCCGG‑GGTT‑AA‑ACCCG‑AGTAGGGCGATTAATGGGCTTT‑CTAGCGTGTCG‑CGGGTTTGCGAACCGCCGGCCTGCGTTGGGCGATAGGGA‑TTTACCATTCCACTCATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTc                                                                                                                                                                             <  2:504532/235‑1 (MQ=255)                                                                     cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑GGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGG‑TGGACCCTTTCAATGATAACCGCCAGCAAAGGGACTT‑TCCAGAGGAT‑GA‑CTTCACGc                                                                                                                                                                         >  2:1288371/1‑238 (MQ=255)                                                                              aCGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAATCGACTAATTTGCTTT‑CCAGCGGGTCG‑AGGGCGTGCGCAACGTCGGCCTGCGGGGAGCGAAAGCGA‑TGTACCCTTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑GTTTCCGCCCCGGGCCg                                                                                                                                                                >  2:1267299/1‑238 (MQ=255)                                                                                  ttACCGAGGTTTCCTTGTCCTTGGCGAG‑AAAAAG‑GGCGACGGGTGCAAAAATCCATC‑AGATGGAAGCGCGACATACAAAGACCGGCGGTT‑AC‑CGCCG‑AA‑AAGTCGACTATTTTGCTTC‑CCAGCGTGTGG‑AGGGCTTGCGAACCGTCGGCCGGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTa                                                                                                                                                            <  2:825562/238‑1 (MQ=255)                                                                                  ttAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑CG‑AAGTCGGATAATTTTCTTT‑CCAGCGCGTCG‑AGGGGGTGCGCAACGGCGGCCCGCGGTGTGCGGGAGCGATTGTACC‑TTTTACTGACAACCGCCCGCCAAGGAACTT‑TCCCGAGGGT‑GACCTTCAC‑CCGCGGGCCGcctca                                                                                                                                                           >  2:127559/1‑235 (MQ=255)                                                                                  ttAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑TCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAAACCCCAGCAAAGGAAATT‑CCCCGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTaa                                                                                                                                                           >  1:958056/1‑239 (MQ=255)                                                                                          ggTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAATCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCCACGTCCGCCTGCGGGGAGCGATAGCGA‑TGGACCCTTTTCCTGATAACCCCCCGCAAAAGAACTT‑CCCCGAGGAT‑GA‑CTTCACGCCACGGGCCCCTTACACGaaaa                                                                                                                                                    >  2:1088685/1‑236 (MQ=255)                                                                                                  gTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑TCAGCGTGTCG‑AGGGCTTGCGCAACGGCGGCCTGCGGGGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑ATTCACGCCACGGGCCGCCTAAACGAATATCCAGCtt                                                                                                                                            >  1:768350/1‑238 (MQ=255)                                                                                                   tACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTTCTTT‑TCAGCGGGGCG‑AGGGGTTGCGCAACGTCGGCCTGCGGTGAGCGGTAGCGA‑TGTACCATTTCACTGATAACCGCCCGCAAAGGAACTT‑CCCCGAGGAT‑GA‑CTTCACGCCACCGGCCGCCTAAACCaaaaatccgcctg                                                                                                                                           >  2:643553/1‑227 (MQ=255)                                                                                                       tgcgcaaaaaaaaC‑GGCGACGGGGGCAAAAATCCATC‑AGATGGAAGCGCGACATACAAAGCCCGGCGGTT‑AA‑CGCCG‑AG‑AGGCCGACTAATTGGCTTT‑CAAGCGTGGCG‑GGGGCTGGCGAAACGGCGGCCTGCGGGGAGCAATAGCGA‑GGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCAt                                                                                                                                       <  2:335686/231‑1 (MQ=255)                                                                                                       tCGGGAG‑AAAAAA‑CGCGGCGGGGCCCCAAAACCATCAAGAGGAAAGCGAGATATA‑GAAGGCCCGGGGTT‑TA‑CGCCGAAG‑AAGTCGATTATTTTGCTTT‑CCGGCGTGGCG‑GGGGTTTCCGCAAGGTCGGCCTGCGGGGGGGGATAGCGA‑TGACCCATTTCACGGAAAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCAt                                                                                                                                       <  2:270813/239‑1 (MQ=255)                                                                                                        tgcgAGAAAAAACGGGGGAGGGGTACCAAATCCAATC‑AGAGGGAGGCGGGACTTATCAAG‑CCGGGGGTTGAC‑CGCGG‑AG‑ATGTCGACTAA‑TTGCTTTCCCGGGGTGTCG‑GGGGTTTCCCAACCGTCGGCGTGGGGGGTGGGATGGGGA‑TGTCCTATTCCCCTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCAt                                                                                                                                       <  2:672822/238‑1 (MQ=255)                                                                                                        cgcgAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGCCTAATTTTCTTT‑TCAGCGGGTGG‑AGGGGTTGCGCAACGTCGGCCCGCGGGGGGCGATAGCGG‑TGTACCATTTCACTTATAACCGCCAGCAAAGGGACTT‑TCCAGAGGATGGA‑CTCCCCCCCCCGGGCC‑CGTTAACCAATATCCACACTGCCATTa                                                                                                                                     >  2:711280/1‑239 (MQ=255)                                                                                                        cgcgAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑AGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGGGTCG‑GGGGGGTGGGCCACGGCGGCCTGCGGGGGGCGATAGCGG‑TGTACCCTTTCAATGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGGT‑GA‑ATTCAAGCCACCGGGCGGCTAAAAGAATATCCCGCCTGCCCTTa                                                                                                                                     >  2:431765/1‑239 (MQ=255)                                                                                                        cgcgAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGCT‑AA‑CGCCC‑AG‑AAATCGACTACTTTGCTTT‑TCAGGGCGTCGTAGGGC‑CGAGCAACAGCGGGCTGCGGTGTGAGCTAGAGG‑TTGACAATTTCACTGGTACACGCCCGAAAAAGAACAT‑GCCAGAGGAT‑GA‑CTGCACGCCAAGGGACGCGTCAACGAATATACCGCCTGCCac                                                                                                                                       >  2:1181303/1‑236 (MQ=255)                                                                                                        cgcgAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGCT‑AA‑AACCC‑AG‑AAGTCGAATAATTTGTTTT‑TCAGCGTGGCG‑GGGGCTTGCGCAAAGTCCGCCTGCGGGGAGCGATAGAGG‑AGGAACATTTTCATCATAACCGCGAGCCAAGGGACGT‑TCCCGAGGAT‑GA‑CTTTACGCCACCGGCCGCCTaaaaagaaatccaacatcgcatta                                                                                                                                     >  2:227534/1‑218 (MQ=255)                                                                                                                 aaaC‑GGCGACGGGTGCAAAAATCCATA‑AGAGGGAAGCGCGACATACCAAGGCCGGCGGTT‑AA‑CGCGG‑AG‑AAGTCGACTAATTTGCTTT‑CCGGCGTGTGG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGGGCGATAGCAA‑TGTACCTTTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGa                                                                                                                             <  2:625618/239‑1 (MQ=255)                                                                                                                               tGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGGTTT‑TCAGCGTGTTG‑AGGGGTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTAACATTTTACTGATAAACCCCAACAAAAGAACTT‑CCCCGAGGCT‑GA‑CTTTCCCCCACGCGCCGCCTCAACCAAAAACCAGCTTGCCCTTATCGTTggcaccggtcgcttgt                                                                                                                >  2:1135099/1‑225 (MQ=255)                                                                                                                                 caaaaaCCCAAC‑AGATGGAAGCGCGAAAAAACAAGACCGGCGGTT‑AC‑CGCCG‑AG‑AAGTCGCCTAATTTGCTTC‑CCAGGGTGTCG‑AGGGCTTGCGAAACGTCGGCCTGGGGTGAGCGATAGCGA‑TGTACCATTTCACTAAAAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGtt                                                                                                               <  2:187718/236‑1 (MQ=255)                                                                                                                                      acccATC‑GGATGGAGGCGCGCATTACCAACAGCGG‑GGTT‑AACCGCGG‑AG‑AGGCCGCCTATTGTGCTTC‑CCGGTGTGGCG‑GGGGCTGGCGCACCGTCGGCCGGGGGGGAGCGATAGCGA‑TGTCCCTTTTCACAGAAAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGt                                                                                                        <  2:183282/237‑1 (MQ=255)                                                                                                                                                          gACATATCAAGACCGGGGGTT‑AG‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑GGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTCCCATTTCACTGAAAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTaa                                                                                     <  1:643553/239‑1 (MQ=255)                                                                                                                                                           aaatatCAGGCCGGGGGGTA‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCGGCGTGTCG‑AGGGTTTGCCAAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTaa                                                                                     <  2:506534/236‑1 (MQ=255)                                                                                                                                                                        ggcggTT‑AA‑CGCGG‑AG‑AAGTCGACTAATTTGTTTC‑CCAGGGTGGGG‑GGGGTTTGCCCACCGTCGGCCGGGGGTGGGCGAGAGCGA‑TGTCCCATTCCCCTGAAAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGt                                                                      <  2:155708/239‑1 (MQ=255)                                                                                                                                                                         gcggTT‑AA‑CGCCG‑GG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑GGGGGTTGCGCAACGGCGGCCTGCGGGGAGCGATAGCGA‑TGTAACATTTCACTGATAAACGCCAGCAAAGGAAATT‑CCCAGAGGGT‑GA‑CTTCACGCCACGGGCCGCCTAAAAGAATATCCCGCCTTCCATTTTCGTTTGACACGGTGGCTTTTTCCCGGGTTAAAG‑GGAACCCGCGTAACCAAGTTaaaaagtt                                                                      >  2:305822/1‑233 (MQ=255)                                                                                                                                                                         gcggTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGGTTT‑CCAGCGCGTCG‑AGGGCCTGCGCCACGTCGGGCTGCGGTGGGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCGGCAAAGGAACTT‑CCCAGAGGGT‑GA‑CTTCACCCCCCCGGCCGCCTAAACGAATATCCCCCCTGCCACTATCGTTGGACACGGTGGCTTGTT‑CCGGGGTAAAACGGAACTCGCGGACCC‑CGTTAAGATAGTc                                                                     >  2:840915/1‑239 (MQ=255)                                                                                                                                                                         gcggTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑TCAGCGTGTCG‑AGGGCTTGCGCAACGGCGGCCCGCGGGGAGCGATAGCGG‑TGTACCATTTCACTGCTAACCCCCAGCAAAAGAACTTCCCCGGAGGA‑‑GA‑ATTCCCCCCCCAGGCCGCTTAAAAGAAAATCCCGCTTTCCATTATCGTTTGACACGGCGGCTTGTT‑TCGGGTTCAAACGGAAACCGCGTAACC‑AATTAAGATAGtt                                                                     >  2:416180/1‑238 (MQ=255)                                                                                                                                                                          cggTT‑AA‑AGCCG‑AG‑AAGTCGACTAATTTTCTTT‑CCAGCGTGTCG‑GGGGCTTGCGCCACGTCGGCCTGCGGTGGGCGGTAGCGG‑TGTTCCATTTTCCTGATAAACGCCCGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAAAGAATATCCAGCCTGCCCTTATCGTTGGACCAGGTGGCTTGTT‑CCGGGGTAAAGCGGAACTCGCGTAACC‑ACGTAAGATAGtag                                                                    >  2:312818/1‑237 (MQ=255)                                                                                                                                                                           ggggtaaccgcgG‑AG‑AGGCCGCCTATTTTGCTTT‑CAGGGGTGTCG‑GGGGCTGGCGAACCGGGGGCCGGCGGGGAGCGATAGCTA‑TGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGt                                                                      <  2:958056/230‑1 (MQ=255)                                                                                                                                                                              t‑AA‑CGCCG‑AG‑AGGTCGCCTATTTTGTTTC‑CCGGGGGGTGG‑GGGGTTTGCGCAACGCCGGCCTGCGGGGGGCAATAGCGT‑GTTACCATTTCACTAATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCt                                                                <  2:887787/239‑1 (MQ=255)                                                                                                                                                                                          g‑aAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTaaagaaag                                                       >  1:566619/1‑239 (MQ=255)                                                                                                                                                                                               cctaCTAATTTCCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGGGGTGGGCGATAGCGA‑TGTACCTTTTCACTAAAAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAggca                                                    <  2:566619/235‑1 (MQ=255)                                                                                                                                                                                                      aTTTGCTTT‑CCATGGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGGACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGaa                                            <  1:1214087/239‑1 (MQ=255)                                                                                                                                                                                                         tgtttt‑CCGGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGGGGTGAGCGATAGCGA‑TGTCCCATTCCACTGATACCCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                                         <  1:1288371/236‑1 (MQ=255)                                                                                                                                                                                                          tgtttccgcgggtggCG‑GGGGCTGGCGAAGCGTCGGCCGGCGGTGGGGGATAGCGA‑TGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                                         <  2:673849/228‑1 (MQ=255)                                                                                                                                                                                                           cTTC‑CCAGCGTGACG‑AGGGCTGGCCCACCGTCGCCCTGCGGTGGGGGATAGCGA‑TGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCa                                       <  1:1200854/239‑1 (MQ=255)                                                                                                                                                                                                                ccAGGGTGTCG‑AGGGCTTGCGCACCGTCGGCCGGCGGTGAGGGATGGCGA‑TGTACCATTCCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGa                                   <  1:220414/239‑1 (MQ=255)                                                                                                                                                                                                                 cAGCGTGTGG‑GGGCTTTGCCCACCGTCGGCCTGCGTGGAGCAATGGCGA‑GTTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGat                                  <  1:437113/239‑1 (MQ=255)                                                                                                                                                                                                                 cAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGat                                  >  1:627893/1‑239 (MQ=255)                                                                                                                                                                                                                       gtgg‑AGGGCTTGCGCAACGCCGGCCTGGGGGGAGCGATGGCGA‑TGTACCTTTACACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtcat                            <  1:127559/236‑1 (MQ=255)                                                                                                                                                                                                                        tCG‑AGGGCTTGCGCAACGTCGGCCTACGGCGAGCGATAGCGA‑TGTACCATTTCACTGATAACAGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATCATCGTTGCACACGGTGACTTGTT‑CCGGGTTAAAGCGGaa                                                                                              >  1:1079744/1‑173 (MQ=255)                                                                                                                                                                                                                            aGGTCTTGCGCAGCGTCGGCCGGGGGGGGGGAATAGCGA‑TGTACCATTTCACTAAAAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTcatcat                         <  2:1268890/238‑1 (MQ=255)                                                                                                                                                                                                                               gCTTGCGCAACGTCGACCGGCGGTGGGCGTGCGCGA‑TTTACCAATCCCCAAATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCt                      <  2:254226/238‑1 (MQ=255)                                                                                                                                                                                                                                tttGCCCAACGTGCGCCGGGGGTGAGCGATGGCGA‑TGTCCCATTTCACTGAAAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTg                     <  1:1324242/237‑1 (MQ=255)                                                                                                                                                                                                                                 ttGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCa                   >  1:798299/1‑239 (MQ=255)                                                                                                                                                                                                                                   gcgcAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                  >  1:178707/1‑238 (MQ=255)                                                                                                                                                                                                                                             ggCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGt          >  2:158918/1‑236 (MQ=255)                                                                                                                                                                                                                                              gCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgct       >  1:687965/1‑238 (MQ=255)                                                                                                                                                                                                                                                cTGCGGTGGTCGATAGCGA‑TGTACCATTTCACGAATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa    <  1:59910/239‑1 (MQ=255)                                                                                                                                                                                                                                                cTGCGGTGGGCGATGGGGA‑TGTCCCATTTCACGGAAAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa    <  1:884123/239‑1 (MQ=255)                                                                                                                                                                                                                                                   cGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAATCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc  >  2:618115/1‑238 (MQ=255)                                                                                                                                                                                                                                                   cGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc  >  2:1078435/1‑238 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 GTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATC  >  NC_000913/989707‑990177

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