Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	754,446	A→G	L8P (CTC→CCC)	gltA ←	citrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	754,446	0	A	G	97.5%	111.6 / NA	40	L8P (CTC→CCC)	gltA	citrate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (27/12);  minor base C (1/0);  total (28/12)
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

GGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAA  >  NC_000913/754313‑754584                                                                                                                                      |                                                                                                                                           gcttgaggTTGAAGTGAATCCTGGGTCAAAGGTTAACTCTCCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACAGCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCtt                                                                                                                                   <  1:242107/137‑1 (MQ=255)            aagtgaagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCCGcc                                                                                                                        >  2:22430/1‑143 (MQ=255)               tgaagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACGTCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTCTCAGCCAtt                                                                                                                     >  1:509397/1‑143 (MQ=255)               tgaagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCAtt                                                                                                                     >  1:596625/1‑143 (MQ=255)                 aagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCGTGCATCACACCCAGTTCAACAGCTGTCTCTCCGTAGAGGGTGGCTTTTGTATTTGTGTCACACATTTa                                                                                                                   >  2:456445/1‑143 (MQ=255)                   gCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCCACAGCTGTATCCCCGTGGAGGGAGGGTTTTGATCTTGTACCCGCCATTCAAg                                                                                                                 >  2:496840/1‑143 (MQ=255)                   gCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGCTGCGGGTGGGTCTTACTTTTGTATCCGCCATCTAAg                                                                                                                 >  2:496295/1‑143 (MQ=255)                     cTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGGTTTGCTTTTGTATCAGCCATTTAAGGt                                                                                                               >  1:336710/1‑143 (MQ=255)                        ggTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCCGCGTGCCATTCAGCACATCCAGTTCAACAGCTGTGTCCCCGCAGAGGGTGCGTTTTGCTTACGTATCAGCCATATAAggacac                                                                                                            >  2:369752/1‑139 (MQ=255)                                gTGATCACAGCTTTTTAGCCTGGAGAACAGATAGCAATACCTTCATGACCCAGCGTGCCTGTCAGCACATCCAGGTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAgcgc                                                                                                    <  1:151335/143‑1 (MQ=255)                                           ttttGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGGGGCGGGTGGGTTTTGCTTTTGTGTCAGCCATTTAAGGTCTCCTTAGCGCCCTCTTggtta                                                                                         >  1:453519/1‑139 (MQ=255)                                                  ccGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGc                                                                                  >  1:908913/1‑143 (MQ=255)                                                   cGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTCTTGCGTAAGACTGcc                                                                                 >  2:1067705/1‑143 (MQ=255)                                                       agTCCGGATATCAATAACATCCTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGaa                                                                             <  1:223615/143‑1 (MQ=255)                                                        gTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCGCCTTAGCGCCTTATTGCGTAAGACTGCCGGAAc                                                                            >  1:872658/1‑143 (MQ=255)                                                        gTACGGAGATGAATAACATCTCGACCCAGCGTGCCTTTCAGGACATCTAGTTCATCAGCTGTATCCCCGTTGAGGGTGGGTTTTGCGTCTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAAc                                                                            <  2:450609/143‑1 (MQ=255)                                                            ggATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAGCATaa                                                                        >  2:146560/1‑143 (MQ=255)                                                            ggATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACCGCCGCACCTTaa                                                                        >  2:642835/1‑143 (MQ=255)                                                            ggATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTGTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTaa                                                                        <  1:402750/143‑1 (MQ=255)                                                             gagatCAATAACAGCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTaaa                                                                       <  2:812141/140‑1 (MQ=255)                                                             gATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGTGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTaaa                                                                       <  2:510074/143‑1 (MQ=255)                                                              atatCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAAt                                                                      <  1:630221/143‑1 (MQ=255)                                                                 tCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTCTCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTCTTTCGTTACACTGCCGGATCTTCAAttct                                                                   >  2:972476/1‑141 (MQ=255)                                                                   aaTAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTATCTTTGcc                                                                 >  1:763459/1‑143 (MQ=255)                                                                   aaTAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTccc                                                                 >  2:738104/1‑143 (MQ=255)                                                                    aTAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCt                                                                >  1:997907/1‑143 (MQ=255)                                                                      agtaTCTTGTGCCAACGTGCCTTTAAGCAAATCTAGTTCAACAGCTATATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTc                                                              <  1:490910/140‑1 (MQ=255)                                                                        cATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGCCTGCCGGCACTTAAATTTGCCTTCGc                                                            >  2:647343/1‑143 (MQ=255)                                                                          tCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGcac                                                          >  2:831825/1‑143 (MQ=255)                                                                            ttGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTCCCTTCGcacct                                                        >  2:446704/1‑141 (MQ=255)                                                                                          gtttGAGCACATCCAGTTCAACTGCTGTGTCACCGTTGAGTGTGGGTTTTGTTTTTGTATTAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCTGGAACTTAAATTTGCCTTCGTACATCAACCTGGCTTTAc                                          <  2:224840/142‑1 (MQ=255)                                                                                               aGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTCAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTCTGCCCGtt                                     >  1:619520/1‑143 (MQ=255)                                                                                                gCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGttt                                    >  1:151410/1‑143 (MQ=255)                                                                                                gCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGttt                                    >  1:151514/1‑143 (MQ=255)                                                                                                  acaTCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTCAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCCttttt                                  >  2:423918/1‑143 (MQ=255)                                                                                                     tCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCCACCTGGCTGTACCCGTGTTTTag                               >  1:179883/1‑142 (MQ=255)                                                                                                          ttCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCCCATCAACCTGGCTTTACCCGTTTTTTATTTGGc                          >  2:717504/1‑143 (MQ=255)                                                                                                                       tCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTATACCCGTTTTTTATTTGGCTCGCCGCTCtgtg             <  2:594237/143‑1 (MQ=255)                                                                                                                             ttGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTCACGTCCCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAagag       <  1:369732/143‑1 (MQ=255)                                                                                                                                  ggTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTACGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGATAGAGGGGaa  >  1:1119119/1‑143 (MQ=255)                                                                                                                                      |                                                                                                                                           GGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAA  >  NC_000913/754313‑754584

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: ‑