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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	989,945	0	A	G	30.2%	69.7 / 8.7	61	S138P (TCA→CCA) ‡	pncB	nicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (22/17);  new base G (6/12);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATC  >  NC_000913/989711‑990177                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               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:935805/239‑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: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)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               CCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATC  >  NC_000913/989711‑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: ‑