Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	124,241	C→T	Q409* (CAG→TAG)	aceE →	pyruvate dehydrogenase E1 component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	124,241	0	C	T	93.2%	113.3 / ‑4.6	45	Q409* (CAG→TAG)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (16/25);  total (17/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGGTCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCT  >  NC_000913/124102‑124378                                                                                                                                            |                                                                                                                                          tggtggTCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAgg                                                                                                                                        <  1:494275/143‑1 (MQ=255) cggtggtCACTGGCCGAAGAAAATGTAGGCTGCATTGAAGACAGCGCAGGAAACCCAAGGCAAAGCGACAGGAATGCTGGCTCATACCATTAACGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAgg                                                                                                                                        <  2:626001/142‑1 (MQ=255)      gtCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAg                                                                                                                                   >  1:543436/1‑143 (MQ=255)              ccGACGAAAATCTACCGTGCATTAAACGAGGCGCAGGAAACCAAAGGCAAAGCGATAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGaa                                                                                                                           <  1:1004364/143‑1 (MQ=255)                   gaaAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAACTGAACATgg                                                                                                                      >  1:576496/1‑143 (MQ=255)                    aaaaTCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGa                                                                                                                     >  1:683865/1‑143 (MQ=255)                      aaTCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACg                                                                                                                   <  1:178280/143‑1 (MQ=255)                               tggaTTCAAGAAAGCGCAGGAAACCAAAGGCAACGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTc                                                                                                          <  1:776449/140‑1 (MQ=255)                                        gAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCcgcg                                                                                                 >  1:605391/1‑143 (MQ=255)                                              gccggACACCACCGGCAAAGCTACAGGACTCCTTGCTCATACCATTAACGGTGACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGtt                                                                                           <  1:197603/140‑1 (MQ=255)                                                  gAAACCAAGGGCAAAGCGACAGTAATCCGTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCaa                                                                                       <  2:570810/143‑1 (MQ=255)                                                  gAAACCAAAGGCCAAGCGACAGTAATCCTTGCTCATACCAAGAATGGCCGCTGCAGAGGCGACGCGGCTGCAGCTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCaa                                                                                       <  1:40005/143‑1 (MQ=255)                                                  gAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCaa                                                                                       >  1:1049610/1‑143 (MQ=255)                                                      ccAAAGGCAAAGCGACAGTAATGCTTGCTCGTACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTCAAAACATCTCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAAtgtg                                                                                   <  2:14515/143‑1 (MQ=255)                                                        aaaGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGcc                                                                                 >  2:827737/1‑143 (MQ=255)                                                         aaGGCAAAGCGACGGTAGTCTTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCg                                                                                <  2:696160/143‑1 (MQ=255)                                                                 gCGACAGTAATCCTTGCTCATCCCATTACAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGa                                                                        <  1:473984/143‑1 (MQ=255)                                                                    aCAGCAATCCTTGCTCATACCATTAAAGGTTCCGGCACGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGc                                                                     <  1:482975/143‑1 (MQ=255)                                                                     cAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCa                                                                    >  2:369626/1‑143 (MQ=255)                                                                        tAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGCCGCGGCTGAAGGTAACAACCTCGCGCACTAGGTTCAGAAACTGAACATGGACGGTGTGCGTCATATCCGCGACCGGTTCAATGTGCCGGTGTCTGATGcgggt                                                                 >  1:524530/1‑139 (MQ=255)                                                                        tAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGGCTGATGCAggt                                                                 >  1:644711/1‑141 (MQ=255)                                                                        tAATCCTTGCGCATGCCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGat                                                                 >  2:815358/1‑143 (MQ=255)                                                                         aaTCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGata                                                                >  1:486136/1‑143 (MQ=255)                                                                           tCCTGGCTCATACCATTGAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATc                                                              <  1:941596/143‑1 (MQ=255)                                                                                gCTCATTCAACTAAAGGTTACGGGATTCGCGTCGCGGCTGAAGGTAAAAAACTCGGGCACGAGGGTAAGAAAATGAACATGTACGGTGTGCGTCATATCCGCGAGCGTTTAAATGTGCCGGTGTCTGATGCAGATATCGaaaa                                                         <  2:639632/143‑1 (MQ=255)                                                                                  tCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAAc                                                       <  1:369626/143‑1 (MQ=255)                                                                                    aTACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCCATGTGCCGGTGCCTGATGCAGATATCGAAAAACTg                                                     >  1:65647/1‑143 (MQ=255)                                                                                    aTACCATTAAAGGTTACGGCATGGGCGACGCGGCGGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACCTGGTCGGTGTGCGTCTTATCCGTGACCGTTTCAATGTGGCAGTGTCTGATGGAAATATCGAAACACTg                                                     >  1:143740/1‑143 (MQ=255)                                                                                     tACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGc                                                    <  2:486136/143‑1 (MQ=255)                                                                                          ttAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTAc                                               <  1:894472/143‑1 (MQ=255)                                                                                            aaaGGTTACGGCATGGGCGACGCGGCGGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTCCGCCATATCCGCGCCCGTTTCAATGCGCCGGTGTCTCATCCAGACATCGAAAAACTGCCGTACAt                                             >  2:1019572/1‑143 (MQ=255)                                                                                             aaGGTTACGCCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATc                                            <  2:459582/143‑1 (MQ=255)                                                                                              aGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATCTCGAAAAACTGCCGTACATCa                                           >  2:182733/1‑143 (MQ=255)                                                                                               ggggaCGGCAGGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCAc                                          <  1:506012/139‑1 (MQ=255)                                                                                               ggTTGCGGCATGGGCGACGCGGGTGAAGGTAAAGACAGCGCGCACTAGGTGAGGAAAATGAACATGGACGGGGTGCGTCATATCCGCGACCGTGTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCAc                                          <  2:594896/143‑1 (MQ=255)                                                                                               ggTTACGGCATGGGCTACGCGGCTGAATGTAAAAACGTCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCAc                                          <  2:985171/143‑1 (MQ=255)                                                                                                              gactaggCTGAATGGATAATCATCTCTCACTAGGGGAAGAACATGAACTTCGTCTGGCTGCGTCATATCCGCGACCGTTTCAACGTTACGGTGTCTCGTGCGGATATCGAAAAACTGCCGTACATCACCTTCCCGGTAGGTTc                           <  2:698475/138‑1 (MQ=255)                                                                                                               aCGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAAGGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCt                          <  2:53749/143‑1 (MQ=255)                                                                                                                cgcgGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTcc                                     >  1:859009/1‑131 (MQ=255)                                                                                                                cgcgGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTcc                                     <  2:859009/131‑1 (MQ=255)                                                                                                                  cgGGTGAAGGTAAAAACATCGCGCACTAGGGTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTTTCGGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGaa                       <  1:707713/143‑1 (MQ=255)                                                                                                                      cgAAGGGAACAGCATCGCGCACCAGGTGAAGAAGATGAACATGGACGGTGGGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATACAGATATCGAAAAGCTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGc                   <  2:182352/142‑1 (MQ=255)                                                                                                                          ggTAAAAACATTGCCCACCACTTTAGTAAGATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCTCCTTCCCGGAAGGTTCTGAAGAGCATAc               <  2:4854/143‑1 (MQ=255)                                                                                                                                    gcgcgCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCAc     <  1:966473/142‑1 (MQ=255)                                                                                                                                       cgcACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATGTGGACGcg  >  1:1002235/1‑142 (MQ=255)                                                                                                                                            |                                                                                                                                          TGGTGGTCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACGCT  >  NC_000913/124102‑124378

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