Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,139,187	T→G	20.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	20.6%	152.3 / 6.9	66	V429G (GTG→GGG)	yegH	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (40/10);  new base G (3/12);  total (44/22)
Fisher's exact test for biased strand distribution p-value = 4.07e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCC  >  NC_000913/2138958‑2139420                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           gTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGGGG‑AAgaga                                                                                                                                                                                                                                   >  2:894554/1‑238 (MQ=255)   tATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGAtag                                                                                                                                                                                                                                >  1:375860/1‑237 (MQ=255)      cGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGGGACATTAAGT‑GACGTCACTAAACCCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGAc                                                                                                                                                                                                                              >  1:368084/1‑238 (MQ=255)             cTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCa                                                                                                                                                                                                                      >  1:1123700/1‑239 (MQ=255)                       ctctGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCag                                                                                                                                                                                                            >  2:1385951/1‑239 (MQ=255)                       ctctGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCag                                                                                                                                                                                                            >  1:250798/1‑239 (MQ=255)                       ctctGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACGGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCag                                                                                                                                                                                                            >  1:412958/1‑239 (MQ=255)                        tctGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGATGGTGTTCCCGGAGACCTTGCCGTTGTTACCTGCCCTGGAGCAGGTCCGTAATGCCCGCACGCTCTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGGGGGAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTACGTTCACGACAGAACGGG‑AAGAAAGCGAAGCCCGCCATGATATTAaga                                                                                                                                                                                                           >  2:1060686/1‑239 (MQ=255)                                gcgAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AACAGATCGACGCCCGCCATGATATTCAGAAGAATGCg                                                                                                                                                                                                   >  2:869408/1‑239 (MQ=255)                                                   ggTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGa                                                                                                                                                                                >  1:1328025/1‑239 (MQ=255)                                                            tcgctcgCCGCTGGGTTTCCCGAAACCCTTGCGTTTGTTACCTGCCCGGGGACAGTCCCGTAATGCCCGCCCGCTTTTTGTTT‑TTGTGGTGGATGAGTTGGGCTCGGAGGAAGGGTTTGTGCCA‑TAAGTAGGGTTCTCGGAACCAATTCCCGATAACTCACCAAAGAATGTGA‑AAGAGATCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCata                                                                                                                                                                       <  2:1123700/233‑1 (MQ=255)                                                                    cGCTGGTGTTCCGGAAACCCTTGCCTTTGTTTCCTGCCTGGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACCTCACTAAACCCATTGCCGGTAACTTACCGAACGAAGGGA‑AAGAGATCGCcgcccg                                                                                                                                                                                                                         <  2:1092206/181‑1 (MQ=255)                                                                       tggttctCCCGGAAGCCTTGCCGTGTTTACTTGCCCTGGCGCAGTCCCGTATGGCCCTCACGCATTTTTTTT‑GTGGGGTGGATGATTTTGGCTTGGGGAAAGGGATTGTACCATTAAGG‑GACGCCACTGAAACCATGGCCGTTATCTTCCCGAAGGAAGGGG‑AAGACCTGGCGGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTgga                                                                                                                                                             <  2:375860/232‑1 (MQ=255)                                                                       tGGTGTTCCCGGAACCCTTGCCGTTGTTACCTCCCCTGGAGCATTTCGGTAATGCCCGCCCGCACTTTGTTT‑TTGTGGTGGATGAGTTGGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GCCGTCACTGAACCCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTgga                                                                                                                                                             <  1:31799/238‑1 (MQ=255)                                                                          tgctCCCGGACCCCTCGCCGTTGTTACCCGCCCGGCACAATTTTGTTCATCCCGGCTCGCACTTTTCTT‑TTGTGGTCGTTGAGTTTGGCTCGGTGGAGGGGTTTGTGAAATAAAGT‑GACGTTACTGAAACCATCGCCGGTAGCTTCCCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGAt                                                                                                                                                         <  2:1006422/236‑1 (MQ=255)                                                                              cccGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATAtt                                                                                                                                                                                                               >  1:1092206/1‑181 (MQ=255)                                                                                                   aCCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGtgccgctgcc                                                                                                                                  >  2:922200/1‑237 (MQ=255)                                                                                                           tGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGACGGGG‑AAGAGATCGACGCCGGCCATGACATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAGTGGTCATATGCCGCTGCAGGGTCGGGTGCAATATGTGCCGCTTCCGCTGGCTGaa                                                                                                                        >  2:1250250/1‑239 (MQ=255)                                                                                                                cattcccGTAATGCCCTCTCTTACTTGGTTT‑TTGGGGTGGATGATTTGTGCGGGGGAGAGGGGTTGGTGATATAAAGT‑ACCGCCACTGAACCCTGTGCCGGTAACTTACCGAAGGAAGTGG‑AAGAGACCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAAc                                                                                                                    <  2:588033/233‑1 (MQ=255)                                                                                                                 tgTTGCGTACGGCCCGCA‑GCATTTTGGTTGTTGGGGGGGATGGGTGGGGCTGGGGGGAAGGGATTGTGACATAAAGA‑GACGTCACAGAACCCATTGCGGTAAACTTCCGGAAGGAGGGGG‑AAGAGACCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGt                                                                                                                  <  2:385412/238‑1 (MQ=255)                                                                                                                  gTTCCTTAATGCCCGCTCGCATTTTGTTT‑TGGTGGTGGAGTATTTTGCCTCGGGGGAAGGAATTGCAAAAGTAAGG‑GCCGCCACAGAAACCTCCGCCGTTAACTCACCAAACGAGGGGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGt                                                                                                                  <  2:15352/238‑1 (MQ=255)                                                                                                                  gTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTg                                                                                                                 >  2:287423/1‑239 (MQ=255)                                                                                                                  gTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTg                                                                                                                 >  1:888842/1‑239 (MQ=255)                                                                                                                   ttCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGa                                                                                                                >  2:134938/1‑239 (MQ=255)                                                                                                                       gTAATGCCAGCGCGCTGTTCTTTT‑TGGGGGGGGATGGGTCTGTCTGGGTGGAATAGATGGTGCAATTAGGT‑GACATGACAAAAACATCTGCCGGTAACTAACCGAACGAGGGAG‑AAGCGACCAACCCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAatat                                                                                                            <  2:791203/239‑1 (MQ=255)                                                                                                                             cccGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACAcc                                                                                                      >  1:849694/1‑239 (MQ=255)                                                                                                                             cccGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACAcc                                                                                                      >  2:1308903/1‑239 (MQ=255)                                                                                                                               cgcacgcaCTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACAAt                                                                                                    >  2:964841/1‑239 (MQ=255)                                                                                                                                  acgcacTTTGCTT‑TTGTGGTGGATGCGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGATCCTGGACGACGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTg                                                                                                  >  1:1257346/1‑238 (MQ=255)                                                                                                                                  acgcacTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAata                                                                                                                                              >  1:885481/1‑194 (MQ=255)                                                                                                                                   cgctttttGTTT‑TTGTGGTGGATGGGTTGGGCGGGGTGGAAGGATTTGGAACTTTAGGT‑ACCGTCAATGAAACCTTGCCCGGTAACTCACCGAAGGAAGGGG‑AAGAGATCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGcc                                                                                                <  2:888842/234‑1 (MQ=255)                                                                                                                                           tgttgtggtggggGAGGGGTTTGGCGGGGGGAAAGGGATTGTGAAATAAAGG‑GGCGCCACAAAACCCAGTCCGTGAAACTCACGAAAGGAGGGGA‑AATAGAGCGCGCCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTAATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGt                                                                                         <  2:1238199/232‑1 (MQ=255)                                                                                                                                            cTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTg                                                                                       >  1:1390357/1‑239 (MQ=255)                                                                                                                                              tttgggggtggATGAGTTGGGCTGGGGGGAAGGGATTGTGAAATTAAGG‑GACGTCACGGAAACCTTGGCGGTTAACTTCCGAAAGGAAGTGG‑AAGAGACCGCCGCCCGCCATGATATTCAGAAGAATGCGGa                                                                                                                                                                                                 <  2:312769/126‑1 (MQ=255)                                                                                                                                              t‑ttGGGGTAGTTTAGTTTGGCGCGGAGGAAGGAATTGTGAAATTAAGT‑GACGTTATTGAAACCTTTGCCGGTAACTTACCGAACGAAGGGG‑AAGATATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGAt                                                                                     <  1:968024/239‑1 (MQ=255)                                                                                                                                                 tggggtggATGAGTTTGGCGCGGTGGAGGGTTTTGTGACATTAAGT‑GACGTCGCTGAAACCATTGCCGGTAACTTCCCGAACGAAGGGG‑AAGAGATCGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATgg                                                                                   <  1:221038/236‑1 (MQ=255)                                                                                                                                                     gtggATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAata                                                                                                             >  1:1193266/1‑209 (MQ=255)                                                                                                                                                        gATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTgg                                                                                                                                                                                       >  1:312769/1‑132 (MQ=255)                                                                                                                                                               ttGGCGCGGTGGAGGGTATTGGGACATTAAGG‑GACGTCACTGAACCAATTGCGGGTAACTTACCAAACGAAGTGG‑AAGAGTTGGCCGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCg                                                                      <  2:1328025/238‑1 (MQ=255)                                                                                                                                                                tGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTa                                                                    >  1:905771/1‑239 (MQ=255)                                                                                                                                                                tGGCTCGGTGGAAGGGATTGTGACATTAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTa                                                                    >  1:386208/1‑239 (MQ=255)                                                                                                                                                                       gTGGAAGGGATTGTGACATTAAGT‑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)                                                                                                                                                                                        aTTAAGT‑CACGTCGCTGAAACCATTGCCGGTAACTTACCGAGCGAAGTGA‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTc                                            <  2:394960/239‑1 (MQ=255)                                                                                                                                                                                         ttAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCa                                           >  2:451120/1‑239 (MQ=255)                                                                                                                                                                                         ttAAGT‑GACGTCACTGAAACCATTGCCGGTAACTTACCGAACGAAGTGG‑AAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCa                                           >  2:932422/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)                                                                                                                                                                                                                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)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           GTTATCGACCTGCTACAACAGTCTCTGCGCGGCGAACCGCTCAACCTGCGGGTGTTGATTCGCCAGCCGCTGGTGTTCCCGGAAACCTTGCCGTTGTTACCTGCCCTGGAGCAGTTCCGTAATGCCCGCACGCACTTTGCTT‑TTGTGGTGGATGAGTTTGGCTCGGTGGAAGGGATTGTGACATTAAGT‑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: ‑