Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,139,187	T→G	22.6%	V429G (GTG→GGG)	yegH →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,139,187	0	T	G	22.6%	150.1 / 8.8	69	V429G (GTG→GGG)	yegH	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (40/11);  new base G (3/14);  total (44/25)
Fisher's exact test for biased strand distribution p-value = 1.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCC  >  NC_000913/2138958‑2139420                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           gTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGGGG‑AAgaga                                                                                                                                                                                                                                   >  2:894554/1‑238 (MQ=255)   tATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGAtag                                                                                                                                                                                                                                >  1:375860/1‑237 (MQ=255)      cGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGGGACATTAAGT‑GACGTCACTAAACCCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGAc                                                                                                                                                                                                                              >  1:368084/1‑238 (MQ=255)             cTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCa                                                                                                                                                                                                                      >  1:1123700/1‑239 (MQ=255)                       ctctGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCag                                                                                                                                                                                                            >  1:250798/1‑239 (MQ=255)                       ctctGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCag                                                                                                                                                                                                            >  2:1385951/1‑239 (MQ=255)                       ctctGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACGGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCag                                                                                                                                                                                                            >  1:412958/1‑239 (MQ=255)                        tctGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGATGGTGTTCCCGGAGACCTTGCCGTTGTTACCTGCCCTGGAGCAGGTCC‑GTAATGCCCGCACGCTCTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGGGGGAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTACGTTCACGACAGAACGGG‑AAGAAAGCGAAGCCCGCCATGATATTAaga                                                                                                                                                                                                           >  2:1060686/1‑239 (MQ=255)                                gcgAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AACAGATCGACGCCCGCCATGATATTCAGAAGAATGCg                                                                                                                                                                                                   >  2:869408/1‑239 (MQ=255)                                                   gtggttggttcgcggccggtgggTGTTCCCGGAACCTTTGCGTTTGTTCCCTGCCCGGGGGC‑GTTCCGGTAGGGCCCGCCGGCTTTTTGTTT‑TTGTGTTGGTGGAGTTTGGCT‑GGTGGGAGGGATTTTTGAAATTAGGT‑GAGGTCACGGAAACCATTGCTGGTACCTTACCGAAGGAGGGGA‑AAGGGATGGCGGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCg                                                                                                                                                                                 <  2:625427/217‑1 (MQ=255)                                                   ggTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGa                                                                                                                                                                                >  1:1328025/1‑239 (MQ=255)                                                            tcgctcgCCGCTGGGTTTCCCGAAACCCTTGCGTTTGTTACCTGCCCGGGGACAGTCCC‑GTAATGCCCGCCCGCTTTTTGTTT‑TTGTGGTGGATGAGTTGGGCTCGGAGGAAGGGT‑TTGTGCCA‑TAAGTAGGGTTCTCGGAACCAATTCCCGATAACTCACCAAAGAATGTGA‑AAGAGATCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCata                                                                                                                                                                       <  2:1123700/233‑1 (MQ=255)                                                                    cGCTGGTGTTCCGGAAACCCTTGCCTTTGTTTCCTGCCTGGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACCTCACTAAACCCATTGCCGGTAACTTACCGAACGAAGGGA‑AAGAGATCGCcgcccg                                                                                                                                                                                                                         <  2:1092206/181‑1 (MQ=255)                                                                       tggttctCCCGGAAGCCTTGCCGTGTTTACTTGCCCTGGCGCAGTCCC‑GTATGGCCCTCACGCATTTTTTTT‑GTGGGGTGGATGATTTTGGCTTGGGGAAAGGGA‑TTGTACCATTAAGG‑GACGCCACTGAAACCATGGCCGTTATCTTCCCGAAGGAAGGGG‑AAGACCTGGCGGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTgga                                                                                                                                                             <  2:375860/232‑1 (MQ=255)                                                                       tGGTGTTCCCGGAACCCTTGCCGTTGTTACCTCCCCTGGAGCATTTCG‑GTAATGCCCGCCCGCACTTTGTTT‑TTGTGGTGGATGAGTTGGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GCCGTCACTGAACCCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTgga                                                                                                                                                             <  1:31799/238‑1 (MQ=255)                                                                          tgctCCCGGACCCCTCGCCGTTGTTACCCGCCCGGCACAATTTTG‑TTCATCCCGGCTCGCACTTTTCTT‑TTGTGGTCGTTGAGTTTGGCTCGGTGGAGGGGT‑TTGTGAAATAAAGT‑GACGTTACTGAAACCATCGCCGGTAGCTTCCCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGAt                                                                                                                                                         <  2:1006422/236‑1 (MQ=255)                                                                              cccGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATAtt                                                                                                                                                                                                               >  1:1092206/1‑181 (MQ=255)                                                                                                   aCCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGtgccgctgcc                                                                                                                                  >  2:922200/1‑237 (MQ=255)                                                                                                           tGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGACGGGG‑AAGAGATCGACGCCGGCCATGACATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAGTGGTCATATGCCGCTGCAGGGTCGGGTGCAATATGTGCCGCTTCCGCTGGCTGaa                                                                                                                        >  2:1250250/1‑239 (MQ=255)                                                                                                                cattccc‑GTAATGCCCTCTCTTACTTGGTTT‑TTGGGGTGGATGATTTGTGCGGGGGAGAGGGGT‑TGGTGATATAAAGT‑ACCGCCACTGAACCCTGTGCCGGTAACTTACCGAAGGAAGTGG‑AAGAGACCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAAc                                                                                                                    <  2:588033/233‑1 (MQ=255)                                                                                                                 tgTTGC‑GTACGGCCCGCA‑GCATTTTGGTTGTTGGGGGGGATGGGTGGGGCTGGGGGGAAGGGA‑TTGTGACATAAAGA‑GACGTCACAGAACCCATTGCGGTAAACTTCCGGAAGGAGGGGG‑AAGAGACCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGt                                                                                                                  <  2:385412/238‑1 (MQ=255)                                                                                                                  gTTCC‑TTAATGCCCGCTCGCATTTTGTTT‑TGGTGGTGGAGTATTTTGCCTCGGGGGAAGGAA‑TTGCAAAAGTAAGG‑GCCGCCACAGAAACCTCCGCCGTTAACTCACCAAACGAGGGGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGt                                                                                                                  <  2:15352/238‑1 (MQ=255)                                                                                                                  gTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTg                                                                                                                 >  1:888842/1‑239 (MQ=255)                                                                                                                  gTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTg                                                                                                                 >  2:287423/1‑239 (MQ=255)                                                                                                                   ttCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGa                                                                                                                >  2:134938/1‑239 (MQ=255)                                                                                                                        gTAATGCCAGCGCGCTGTTCTTTT‑TGGGGGGGGATGGGTCTGTCTGGGTGGAATAGA‑TGGTGCAATTAGGT‑GACATGACAAAAACATCTGCCGGTAACTAACCGAACGAGGGAG‑AAGCGACCAACCCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAatat                                                                                                            <  2:791203/239‑1 (MQ=255)                                                                                                                              cccGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACAcc                                                                                                      >  2:1308903/1‑239 (MQ=255)                                                                                                                              cccGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACAcc                                                                                                      >  1:849694/1‑239 (MQ=255)                                                                                                                                cgcacgcaCTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACAAt                                                                                                    >  2:964841/1‑239 (MQ=255)                                                                                                                                   acgcacTTTGCTT‑TTGTGGTGGATGCGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGATCCTGGACGACGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTg                                                                                                  >  1:1257346/1‑238 (MQ=255)                                                                                                                                   acgcacTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAata                                                                                                                                              >  1:885481/1‑194 (MQ=255)                                                                                                                                    cgctttttGTTT‑TTGTGGTGGATGGGTTGGGCGGGGTGGAAGGAT‑TTGGAACTTTAGGT‑ACCGTCAATGAAACCTTGCCCGGTAACTCACCGAAGGAAGGGG‑AAGAGATCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGcc                                                                                                <  2:888842/234‑1 (MQ=255)                                                                                                                                            tgttgtggtggggGAGGGGTTTGGCGGGGGGAAAGGGA‑TTGTGAAATAAAGG‑GGCGCCACAAAACCCAGTCCGTGAAACTCACGAAAGGAGGGGA‑AATAGAGCGCGCCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTAATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGt                                                                                         <  2:1238199/232‑1 (MQ=255)                                                                                                                                             cTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTg                                                                                       >  1:1390357/1‑239 (MQ=255)                                                                                                                                               tttgggggtggATGAGTTGGGCTGGGGGGAAGGGA‑TTGTGAAATTAAGG‑GACGTCACGGAAACCTTGGCGGTTAACTTCCGAAAGGAAGTGG‑AAGAGACCGCCGCCCGCCATGATATTCAGAAGAATGCGGa                                                                                                                                                                                                 <  2:312769/126‑1 (MQ=255)                                                                                                                                               t‑ttGGGGTAGTTTAGTTTGGCGCGGAGGAAGGAA‑TTGTGAAATTAAGT‑GACGTTATTGAAACCTTTGCCGGTAACTTACCGAACGAAGGGG‑AAGATATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGAt                                                                                     <  1:968024/239‑1 (MQ=255)                                                                                                                                                  tggggtggATGAGTTTGGCGCGGTGGAGGGTT‑TTGTGACATTAAGT‑GACGTCGCTGAAACCATTGCCGGTAACTTCCCGAACGAAGGGG‑AAGAGATCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATgg                                                                                   <  1:221038/236‑1 (MQ=255)                                                                                                                                                      gtggATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAata                                                                                                             >  1:1193266/1‑209 (MQ=255)                                                                                                                                                         gATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTgg                                                                                                                                                                                       >  1:312769/1‑132 (MQ=255)                                                                                                                                                                ttGGCGCGGTGGAGGGTA‑TTGGGACATTAAGG‑GACGTCACTGAACCAATTGCGGGTAACTTACCAAACGAAGTGG‑AAGAGTTGGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCg                                                                      <  2:1328025/238‑1 (MQ=255)                                                                                                                                                                 tGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTa                                                                    >  1:905771/1‑239 (MQ=255)                                                                                                                                                                 tGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTa                                                                    >  1:386208/1‑239 (MQ=255)                                                                                                                                                                        gTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGAGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGGCGGCGAATGGTCATATGCCGCTGGAGGATCTGGGGCCATATGTGCCGCTGCCGCTGGATGACAAACGTGAATATCCCACCATTGCGGGGCTGTTGATGGAGTATTTGCAGCGTATTCCaga                                                             >  2:866261/1‑237 (MQ=255)                                                                                                                                                                                   ttGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCgaa                                                   >  2:1187327/1‑239 (MQ=255)                                                                                                                                                                                    tgtgACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCTTGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCgaag                                                  >  1:919840/1‑239 (MQ=255)                                                                                                                                                                                    tgtgACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCgaag                                                  >  2:695997/1‑239 (MQ=255)                                                                                                                                                                                    tgtgACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCgaag                                                  >  1:89382/1‑239 (MQ=255)                                                                                                                                                                                    tgtgACATTAAGC‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCgaag                                                  >  1:39321/1‑239 (MQ=255)                                                                                                                                                                                    tgtcACATTAAGG‑GACGTCACTAAACCCATGGCCGGTAACTTACGAAACGAAGTGG‑AAGAGACCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCgaag                                                  <  1:894554/239‑1 (MQ=255)                                                                                                                                                                                          aataaagggccgccacgaaaaCCTTTCCCGTTAACTTCGCGAAGGAAGTGG‑AAGAGAGCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGtt                                             <  1:869408/221‑1 (MQ=255)                                                                                                                                                                                          aTTAAGT‑CACGTCGCTGAAACCATTGCCGGTAACTTACCGAGCGAAGTGA‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTc                                            <  2:394960/239‑1 (MQ=255)                                                                                                                                                                                           ttAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCa                                           >  2:932422/1‑239 (MQ=255)                                                                                                                                                                                           ttAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCa                                           >  2:451120/1‑239 (MQ=255)                                                                                                                                                                                            tAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAg                                          >  2:796527/1‑239 (MQ=255)                                                                                                                                                                                                   aCGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTgggg                                    >  1:563121/1‑239 (MQ=255)                                                                                                                                                                                                     gTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAatat                                                                                                            >  2:1312344/1‑165 (MQ=255)                                                                                                                                                                                                     gTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAatat                                                                                                            <  1:1312344/165‑1 (MQ=255)                                                                                                                                                                                                     gTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTgggggt                                  >  2:810547/1‑238 (MQ=255)                                                                                                                                                                                                        aCTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGatta                               >  2:1252602/1‑239 (MQ=255)                                                                                                                                                                                                         cagAAACTTTTGCGGGTACTTTCCCGAACGAAGGGG‑AAGAGAGCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGattat                              <  1:1269079/237‑1 (MQ=255)                                                                                                                                                                                                         ataaacccATGGCCTGTACTTTACC‑AACCAAGGGGAAAGAGTGCGCCCCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGattat                              <  1:246400/233‑1 (MQ=255)                                                                                                                                                                                                             acccATTGCGGGTAACTTACCGAAAGAGGGGG‑AAGAGATCGACGGCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGattat                              <  1:1362672/233‑1 (MQ=255)                                                                                                                                                                                                                  ttgcagatacctgcccGAACGAAGGGA‑AAGGACTCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTaaa                     <  1:1208154/225‑1 (MQ=255)                                                                                                                                                                                                                  ttGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTggggg                                   >  1:101455/1‑225 (MQ=255)                                                                                                                                                                                                                  ttGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTaaa                     >  2:202670/1‑239 (MQ=255)                                                                                                                                                                                                                   tcccGATAACTTGCTGAACGAAGTGG‑AAAGGGTCCACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGTAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTaaa                     <  1:1385951/236‑1 (MQ=255)                                                                                                                                                                                                                   tGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTaaa                     >  1:912377/1‑238 (MQ=255)                                                                                                                                                                                                                                   aacgaaGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGaaa     >  1:1134725/1‑238 (MQ=255)                                                                                                                                                                                                                                     cgaaGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGcc  >  2:140472/1‑239 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           GTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCC‑GTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGA‑TTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCC  >  NC_000913/2138958‑2139420

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