Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,875,235	A→G	8.1%	D312D (GAT→GAC)	ilvA‑II ←	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,875,235	0	A	G	8.1%	140.4 / 2.8	49	D312D (GAT→GAC)	ilvA‑II	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base A (27/18);  new base G (2/2);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

ACTCGGTGATCTGACGCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTATCGTCGTAGATATCCTTGATT  >  NC_002947/5875102‑5875382                                                                                                                                      |                                                                                                                                                    aCTCGGTGATCTAACGCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGCTCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTg                                                                                                                                      >  2:61692/1‑149 (MQ=255)   tCGGTGATCTGACGCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTgg                                                                                                                                     >  6:68768/1‑148 (MQ=255)            tGACGCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGATCGGCACGTTCGGCCACATGGCGCAGGCGAGCGAAATTGACGTTGCCACCGCAGTc                                                                                                                           >  8:252320/1‑149 (MQ=255)             gACGCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGAGCGAAATTGACGTTGGCACCGGAGTc                                                                                                                           >  8:105611/1‑148 (MQ=255)             gACGCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCACGGCACGTTCGGCCACATGGCGCAGGCGATCGAAAATGACGTTGGCACCGGAGTc                                                                                                                           >  3:224518/1‑148 (MQ=255)              aCGCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGt                                                                                                                            <  5:37484/146‑1 (MQ=255)              aCGCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGa                                                                                                                         >  4:75583/1‑149 (MQ=255)              aCGCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGa                                                                                                                         >  3:311881/1‑149 (MQ=255)                gCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATg                                                                                                                       <  4:311881/149‑1 (MQ=255)                                gCAGAAGGCCTTGAAGCTACCCGGGCGCTCATGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTgg                                                                                                        <  7:252320/148‑1 (MQ=255)                                         cTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACg                                                                                               <  6:180077/148‑1 (MQ=255)                                         cTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACg                                                                                               >  5:180077/1‑148 (MQ=255)                                              aGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGtat                                                                                         <  1:213095/149‑1 (MQ=255)                                                 tACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGc                                                                                      >  4:101065/1‑149 (MQ=255)                                                   cccGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCCCCCAGCTCGGCCCGTTCGGCCACCTGGCGCAGGCGGTCGaa                                                                                                                                                >  7:26543/1‑89 (MQ=18)                                                   cccGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGGTCGaa                                                                                                                                                >  7:26536/1‑89 (MQ=255)                                                   cccGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGGTCGaa                                                                                                                                                <  8:26536/89‑1 (MQ=255)                                                   cccGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGGTCGaa                                                                                                                                                <  8:26543/89‑1 (MQ=255)                                                       ggCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACg                                                                                <  5:198703/149‑1 (MQ=255)                                                            tCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTActt                                                                           >  8:284025/1‑149 (MQ=255)                                                               gggATGGTCACGGCGATGATCGCCTCGCTCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTActtctt                                                                        <  3:75583/149‑1 (MQ=255)                                                                  aTGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTGGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGAt                                                                     >  1:17794/1‑149 (MQ=255)                                                                        aCGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGc                                                               >  7:106084/1‑149 (MQ=255)                                                                        aCGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGc                                                               >  1:266513/1‑149 (MQ=255)                                                                        aCGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGc                                                               >  3:296727/1‑149 (MQ=255)                                                                          ggCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCca                                                             >  4:189518/1‑149 (MQ=255)                                                                           gCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCca                                                             >  6:333561/1‑148 (MQ=255)                                                                           gCGATGAGCGCCTCGCGCTTTTCACCGAGCTCGGCACGTTCGGCCTCAGGGGGCAGGCGATCGAAATTGACGTTGGCAGCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCcac                                                            <  4:224518/149‑1 (MQ=255)                                                                                gatCGCCTCGCGCTTTTAACCCAGCTCGGCACGTTCGGCCACATTGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTGGACGTACTTCTTGATGCCAGCCACAcc                                                         >  8:141512/1‑147 (MQ=255)                                                                                    gCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAgcgc                                                   <  5:333561/149‑1 (MQ=255)                                                                                       tCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAgcgcgc                                                 >  4:21322/1‑148 (MQ=255)                                                                                       tCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAgcgc                                                   <  1:147553/146‑1 (MQ=255)                                                                                       tCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAgcgc                                                   >  2:147553/1‑146 (MQ=255)                                                                                         gcgcTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGAt                                                                                                                        >  6:292218/1‑75 (MQ=255)                                                                                             ttttCACCCAGCTCGGCAAGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAgg                                          >  5:354785/1‑149 (MQ=255)                                                                                                  aCCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCgg                                     <  2:17794/149‑1 (MQ=255)                                                                                                    ccAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGt                                    >  5:295720/1‑148 (MQ=255)                                                                                                        cTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACAGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTAAACCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATc                                >  5:39830/1‑148 (MQ=255)                                                                                                            gCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCg                           >  4:331140/1‑149 (MQ=255)                                                                                                                    ggCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTAtcg                    <  1:61692/148‑1 (MQ=255)                                                                                                                      ccACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTAtcgtcg                 <  5:68768/149‑1 (MQ=255)                                                                                                                        acaTGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTATCGTCGTa               <  7:284025/149‑1 (MQ=255)                                                                                                                         caTGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTATCGTCGTAg              <  2:266513/149‑1 (MQ=255)                                                                                                                         caTGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTATCGTCGTAg              <  3:101065/149‑1 (MQ=255)                                                                                                                          aTGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTATCGTCGTAGa             <  3:189518/149‑1 (MQ=255)                                                                                                                             gcgcAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGAGCCAGCAGTCTCGGTGAACGAGCGGGTATCGTCGTAGatat          >  8:277974/1‑149 (MQ=255)                                                                                                                              cgcAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACACCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTATCGTCGTAGATATc         >  5:217194/1‑149 (MQ=255)                                                                                                                                   gCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTATCGTCGTAGATATCCTTGa    <  3:21322/149‑1 (MQ=255)                                                                                                                                     gATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTATCGTCGTAGATATCCTTGAtt  >  6:292955/1‑149 (MQ=255)                                                                                                                                      |                                                                                                                                                    ACTCGGTGATCTGACGCTTGCCGATAGCCTCGCAGAAGGCCTTGAAGCTACCCGGGCGCTCAGGGATGGTCACGGCGATGATCGCCTCGCGCTTTTCACCCAGCTCGGCACGTTCGGCCACATGGCGCAGGCGATCGAAATTGACGTTGGCACCGGAGTCGATGGCCACCAAGGTCTGGCCGCTGACGCCGTATAGCTCGACGTACTTCTTGATGCCCGCCACACCCAGCGCGCCAGCAGGCTCGGTGATCGAGCGGGTATCGTCGTAGATATCCTTGATT  >  NC_002947/5875102‑5875382

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