Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,147,792	T→A	S58T (TCG→ACG)	plsX →	putative phosphate acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,147,792	0	T	A	88.9%	59.8 / 2.5	27	S58T (TCG→ACG)	plsX	putative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (12/12);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

GCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGC‑GTCCGGCGCTGGTGACGGTATTACCACATCAGCAAA  >  NC_000913/1147551‑1148035                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     gCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGt                                                                                                                                                                                                                                                 <  1:102763/247‑1 (MQ=255)   gTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGc                                                                                                                                                                                                                                             >  1:322589/1‑249 (MQ=255)        gaAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGattat                                                                                                                                                                                                                                      <  1:84701/251‑1 (MQ=255)                             ccAGGCAACTGGGGAAAGACCAAAACGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGc                                                                                                                                                                                                                 <  1:83846/251‑1 (MQ=255)                                   aaCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGa                                                                                                                                                                                                           <  1:206207/251‑1 (MQ=255)                                                  aaaCCGGGAGGCGACGCTACCTTGACGCGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGAAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGACGGCAATTCCGACGCCATCACGACATTACTTGTTAAAGCTTACTTTGAACATCGTAAGCGTCTGCAGATTATTCCTGCGCAGACAGTTATCGCACGTGATGCAAGGCCttagc                                                                                                                                                                                             >  2:366471/1‑247 (MQ=255)                                                   aaCCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCCGTTATCGCCAGTGATGCCCGGCCTTCGCa                                                                                                                                                                                            >  1:237711/1‑250 (MQ=255)                                                       gggCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTa                                                                                                                                                                                       <  1:363949/251‑1 (MQ=255)                                                                     ccTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGt                                                                                                                                                                         >  1:305064/1‑251 (MQ=255)                                                                                      cTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAATCAGtt                                                                                                                                                                                                                      >  2:18396/1‑189 (MQ=255)                                                                                      cTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAATCAGtt                                                                                                                                                                                                                      <  1:18396/189‑1 (MQ=255)                                                                                                aTGTCCTGGGAGGGGATTTTAGCCCTCCCGTGACAGTGCCTGCAGCATTGCAGGCACTGTATTCTAATTCGCAACTCACTCTTCTTTTAGTAGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTAAGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCATTCGCAAGCTATCCGCGCCCGTCGTGGGAGTTAAATGCGCGTGGCCCTgg                                                                                                                                                >  2:114413/1‑249 (MQ=255)                                                                                                aTGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTgg                                                                                                                                                >  1:284499/1‑249 (MQ=255)                                                                                                   tCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGtt                                                                                                                                                                                                                      <  1:324100/176‑1 (MQ=255)                                                                                                   tCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGtt                                                                                                                                                                                                                      >  2:324100/1‑176 (MQ=255)                                                                                                                                         gcagcaTTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTCCTTGCTAAAGCTAACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCCAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGCGCTGGTGAAAGAGGGTCGAGCGCAAGCCTGTGTCacgatcgg                                                                                                     >  2:384654/1‑244 (MQ=255)                                                                                                                                         gcagcaTTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGAGCAAGCCTGTGTCa                                                                                                            <  1:384654/244‑1 (MQ=255)                                                                                                                                              aTTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATa                                                                                                <  1:327618/251‑1 (MQ=255)                                                                                                                                                           aatcttAATCCGCAACACTCTCTTTTTATAGCCGGCAAGACCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCTTCTGCAGATTATTCCGGCGCAGTCAGTTATCGCCAGTGATGCCCTGCCTTCGCAAGCTACCCGCGCCAGTCGTGTGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTTTCAGTGCCGGTAATTCCGGGGCGCTGAt                                                                                   <  2:260169/246‑1 (MQ=255)                                                                                                                                                                                     tttAGTCGGCAATTCCGCCGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTACTGCGCAGTCAGTTATCGCCAATGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTAGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGTTgg                                                                                 >  2:186772/1‑227 (MQ=255)                                                                                                                                                                                     tttAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAATGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATgg                                                                                 <  1:186772/227‑1 (MQ=255)                                                                                                                                                                                                   ccGACGCCATCCCGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGAtt                                           <  1:77573/251‑1 (MQ=255)                                                                                                                                                                                                     gACGCCATCACGCCATTATTGGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                         <  1:176655/251‑1 (MQ=255)                                                                                                                                                                                                                                     cTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGATGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAAC‑TTCCGGCGCTGGAGACGGTATTACCACATc        >  2:152846/1‑251 (MQ=255)                                                                                                                                                                                                                                       ttGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCACGCTATACGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGGGAAAGAAGGTCGAGCGCAAGCCTGTTTCAGTGCCGGAAATACCGGGGCACTGATGGAGCTGGCAAGATTATTACTCAAGCCCCTGGAGGGGATTGGGC‑GTCCGGCGCTGGTGACGGTATTACcactac        >  2:316674/1‑246 (MQ=255)                                                                                                                                                                                                                                           acaacGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGC‑GTCCGGCGCTGGTGACGGTATTACCACATCAGCaaa  <  1:161701/251‑1 (MQ=255)                                                                                                                                                                                                                                            caacGTACGCGACTGCAGATAATTCCTGAGCAGTCAGTTATAGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGGGGGAGTTCAATGCGCATGGCCCTGGAGCTGATGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGGGCCGGTAATCCAGGGGAGCTGAAGGGGGTGGGAAAATAATTAAACAAGACCATGGAAGGGATT‑AGAAGTCCGTCGCTGGTGACGGTATTAcccactcaacaaa  >  2:229595/1‑239 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     GCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGC‑GTCCGGCGCTGGTGACGGTATTACCACATCAGCAAA  >  NC_000913/1147551‑1148035

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