Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	989,953	C→G	9.3%	R135P (CGC→CCC)	pncB ←	nicotinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	989,953	0	C	G	9.3%	122.9 / 2.3	61	R135P (CGC→CCC)	pncB	nicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/27);  new base G (3/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACC  >  NC_000913/989716‑990186                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                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)                                    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)                                                                                                                                                                                                                                                  gatagcga‑tGTCCCTTTTCCCTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAAcc  <  1:1307475/239‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                GGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAG‑AAAAAC‑GGCGACGGGTGCCAAAATCCATC‑AGATGGAAGCGCGACATATCAAGACCGGCGGTT‑AA‑CGCCG‑AG‑AAGTCGACTAATTTGCTTT‑CCAGCGTGTCG‑AGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGA‑TGTACCATTTCACTGATAACCGCCAGCAAAGGAACTT‑CCCAGAGGAT‑GA‑CTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTT‑CCGGGTTAAAGCGGAACTCGCGTAACC‑AGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTGAACC  >  NC_000913/989716‑990186

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