Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,148,972	C→T	Q72* (CAG→TAG)	fabH →	3‑oxoacyl‑[acyl carrier protein] synthase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,148,972	0	C	T	97.1%	113.5 / NA	35	Q72* (CAG→TAG)	fabH	3‑oxoacyl‑[acyl carrier protein] synthase 3
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/19);  minor base A (1/0);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACCAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGGCGCGCACCTGCGATCCAACCGATCGTG  >  NC_000913/1148734‑1149212                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 ataACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTg                                                                                                                                                                                                                                       <  2:202656/250‑1 (MQ=255) ataACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGa                                                                                                                                                                                                                                      <  2:43214/251‑1 (MQ=255)   aaCCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGAt                                                                                                                                                                                                                                     <  1:11528/250‑1 (MQ=255)    aCCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCg                                                                                                                                                                                                                                   >  1:316948/1‑251 (MQ=255)     ccGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGAGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACACGTACCGGGATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGTTTGGCCTGATCGt                                                                                                                                                                                                                                  >  1:314989/1‑251 (MQ=255)           aGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTGGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGc                                                                                                                                                                                                                            <  2:415489/251‑1 (MQ=255)                      gTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTg                                                                                                                                                                                                                 >  2:408304/1‑251 (MQ=255)                                                     aGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCa                                                                                                                                                                                  <  1:277153/251‑1 (MQ=255)                                                                      aGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACGTTGACGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCCATTGAAATGGCGGGCATTGAGAAAGCCTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTCCGCACGCTTTCCCGAGCTCAGCTTGTCAGATTCCAAGCCTGTTgg                                                                                                                                                                 >  2:234112/1‑251 (MQ=255)                                                                            gACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCAtt                                                                                                                                                            >  1:456583/1‑250 (MQ=255)                                                                                aaCGCCGATTTGGAAAAAATGGTGGACACCTCTGCCGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATc                                                                                                                                                                                                                                    >  2:198475/1‑174 (MQ=255)                                                                                aaCGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAATCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATc                                                                                                                                                                                                                                    <  1:198475/174‑1 (MQ=255)                                                                                    ccGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGc                                                                                                                                                   <  2:11129/251‑1 (MQ=255)                                                                                                    ggtggACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCAGGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTg                                                                                                                                   <  2:456583/251‑1 (MQ=255)                                                                                                        gaccccTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTgcagc                                                                                                                               <  1:331858/247‑1 (MQ=255)                                                                                                                 gACGAGTGGATTGTCACTCGTACCGGTATCCGCGCACGACACATTGCCGCGCCCAACGAAACCGTTTCAACCATGGGCTTTGAAACGGCGACAAGCGCAATTGAGATGGCGGGCATTGCGTAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCCACGCATGTTGGGCATTAAAGGCTGCCCGGCATTTGACGTTGCAGCAGCCTgcgc                                                                                                                      >  1:374460/1‑251 (MQ=255)                                                                                                                        ggATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTc                                                                                                               <  1:149483/251‑1 (MQ=255)                                                                                                                                    gTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCGGCGCAGGTTTCACCTATGCAtt                                                                                                    >  2:131695/1‑250 (MQ=255)                                                                                                                                     tACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGATTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTACAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTa                                                                                                   >  1:291248/1‑250 (MQ=255)                                                                                                                                        cGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCg                                                                                               <  1:408304/251‑1 (MQ=255)                                                                                                                                            aTCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGc                                                                                           >  1:12027/1‑251 (MQ=255)                                                                                                                                            aTCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGCTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGc                                                                                           >  2:11999/1‑251 (MQ=255)                                                                                                                                                 cgAACGCCACATTCCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGAGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACAAGATTGGCCTGATCGTTGAGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGCTTCAAAGCATGGTGGGCATTAACGGTTGCCCGGCATTTGACGTTGCAGCAGCCTACGCAGGTTTCACCTATGCCTTACGCGTAGCCGATc                                                                                      >  2:334143/1‑251 (MQ=255)                                                                                                                                                        cacaTTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCCCCTATGCATTAAGCGTAGCCGATCAATACGt                                                                               <  2:192461/251‑1 (MQ=255)                                                                                                                                                                               cGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAg                                                                                                                                                                                         >  1:362032/1‑122 (MQ=255)                                                                                                                                                                               cGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAg                                                                                                                                                                                         <  2:362032/122‑1 (MQ=255)                                                                                                                                                                                       ccATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGgtc                                                <  1:323790/251‑1 (MQ=255)                                                                                                                                                                                               tttgaacgcgcgACACGCCCACTAGAGATGCCCGGAATTCAGAAAGCCTAGATTGGCCTGCTGGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGCGCGCAGCTTGTCAGACTCAAAGAATGTTGGGCCTTACAGGTTACCCGGCATTTGACGTTGCAGCAGACTGCGCAGGTTTCCCCTATGCATTAAGCGTAGCCGAACATTACGAGAGCTCTGGGGCGGGGAAGTATGACCTGGGCGTCCGTTc                                        >  2:177393/10‑251 (MQ=255)                                                                                                                                                                                                 tGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTcc                                       <  2:12027/250‑1 (MQ=255)                                                                                                                                                                                                       ggcgACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAACGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTAc                                >  1:67132/1‑251 (MQ=255)                                                                                                                                                                                                           acacGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTCGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGgc                            <  2:383550/251‑1 (MQ=255)                                                                                                                                                                                                           acacGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGAACGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGgc                            <  1:234112/251‑1 (MQ=255)                                                                                                                                                                                                                                     tGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGaaa                                                                           >  2:32821/1‑178 (MQ=255)                                                                                                                                                                                                                                     tGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGaaa                                                                           <  1:32821/178‑1 (MQ=255)                                                                                                                                                                                                                                     ggagaAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCTCGCTTTCCCGGGTGCATCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGTCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCTTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCGCTGGTCGTCGGTTCCGATGTACTGGCGCGCACCTGCGATCCAACCGATCGTg  <  2:161817/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 ATAACCGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACCAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGGCGCGCACCTGCGATCCAACCGATCGTG  >  NC_000913/1148734‑1149212

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