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	96.5%	93.7 / ‑6.4	29	S58T (TCG→ACG)	plsX	putative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (10/18);  total (10/19)
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.40e-01

GTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAG  >  NC_000913/1147553‑1148031                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 gTGATGAAGCTTAGTGAGGATTTTCCCCAGTCAACGGGGGAAAGACCAAACCGGGCGGCTACGATACCTTGACACGTCTAGCCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAAGTCTAATTCGCAACTCACTCTTCTTTTATTCGGCAATTCCGACGCCACCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAg                                                                                                                                                                                                                                       <  1:213457/251‑1 (MQ=255) gTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGc                                                                                                                                                                                                                                         >  1:149019/1‑249 (MQ=255)                    aTTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGcagtcag                                                                                                                                                                                                                    <  1:398088/251‑1 (MQ=255)                                     ggggAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGc                                                                                                                                                                                                    <  2:458840/250‑1 (MQ=255)                                      gggAAAGACCAAACCGGGCGGCGCCGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGccc                                                                                                                                                                                                  <  2:140734/251‑1 (MQ=255)                                      gggAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTGGCTAAAGCTGACTTTGAACAACGGACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGccc                                                                                                                                                                                                  <  1:506885/251‑1 (MQ=255)                                           aGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCt                                                                                                                                                                                             >  1:396669/1‑251 (MQ=255)                                                 aaCCGGGCGGCGACGCTACCTTGACACGTCTAACCCTGGCGTTAGATGTCAAGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCAATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCa                                                                                                                                                                                        >  1:351983/1‑250 (MQ=255)                                                            gacgaTACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCgcgc                                                                                                                                                                            >  2:367326/1‑251 (MQ=255)                                                             acgaTACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACCGAATTCTAATTCGCAACTCACTCTACTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGACAGTGATGCCCGGCCTTCGCAAGCTATCCGCGcc                                                                                                                                                                           >  1:234818/1‑251 (MQ=255)                                                                                        cGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCGCTCTTCTTTTAGTCGTCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGGCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTGCAATGCGCGTGG‑cc                                                                                                                                               <  2:194918/251‑1 (MQ=255)                                                                                           tAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCt                                                                                                                                             <  1:304170/250‑1 (MQ=255)                                                                                                  cATGGGAGGGGATTTTGGCCCTTCCGGTACAGTGCCGTCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCGGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTGCAATGCGCGTGG‑CCCTGGAGCTg                                                                                                                                      <  2:322752/250‑1 (MQ=255)                                                                                                                  ggCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAgc                                                                                                                     <  2:216132/251‑1 (MQ=255)                                                                                                                             gACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATGCCGACGCCATCACGCCATTACGTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCGATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGGCGAGCGCAAGCCtgtg                                                                                                          <  2:542798/251‑1 (MQ=255)                                                                                                                                  tGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGTTGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGt                                                                                                     >  1:392279/1‑251 (MQ=255)                                                                                                                                   gCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGt                                                                                                     <  1:165900/250‑1 (MQ=255)                                                                                                                                              tGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATAc                                                                                          <  1:72102/250‑1 (MQ=255)                                                                                                                                               gCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCg                                                                                        <  1:75308/251‑1 (MQ=255)                                                                                                                                                cAGGCACTGAATTCTAATTCGCCACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTCAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCCGTGCCGGTAATACCgg                                                                                       >  2:487062/1‑251 (MQ=255)                                                                                                                                                  ggCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCgggg                                                                                     >  1:314749/1‑251 (MQ=255)                                                                                                                                                        gAATGCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGGCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTg                                                                                <  1:487062/250‑1 (MQ=255)                                                                                                                                                               aaTTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGc                                                                         <  2:91915/250‑1 (MQ=255)                                                                                                                                                                                tcttTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTc                                                        >  1:480121/1‑250 (MQ=255)                                                                                                                                                                                    gtggTCGGCAAGTCCGAGGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGGTTGCAGAGTATGCCTGCGGAGTCAGTTATCGCCAGTGATGCCCGGCCGTCGCAAGCTATCCGCGCCAGGCGTGGTAGTTCAAGGCG‑GTGGGCCCTGGAGCTGGGGAAAGAGGGTCGGGCGCAAGCCTGTGTCAGTGCCGGTAATTCCGGGGCGCTGGTGGGGCTGGCAAAATTATTACTCAAGc                                                    <  2:192647/247‑1 (MQ=255)                                                                                                                                                                                         cGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCTAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTgg                                              >  1:451827/1‑251 (MQ=255)                                                                                                                                                                                              aTTCCGACGCCATCCCGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAgggg                                         <  2:396669/251‑1 (MQ=255)                                                                                                                                                                                                  cgcctccATCTCGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTg                                     <  2:245207/246‑1 (MQ=255)                                                                                                                                                                                                                                      tGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAg  <  1:221178/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 GTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGG‑CCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAG  >  NC_000913/1147553‑1148031

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