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.9%	135.4 / NA	47	L8P (CTC→CCC)	gltA	citrate synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (19/27);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAA  >  NC_000913/754315‑754585                                                                                                                                    |                                                                                                                                            aggCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCtttt                                                                                                                                  <  2:1037182/141‑1 (MQ=255)     ggTCGAAGTCAAGCCTGGGTCAAAGGTCAACACACCTTTTGAACCGAGAGTACGGATATCAATCACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTGTTGCTTTTGTAt                                                                                                                              <  1:19693/143‑1 (MQ=255)         taagtgaagCCTGGGTCAAAGGTGAACATACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGc                                                                                                                          <  1:867416/142‑1 (MQ=255)                 gCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAg                                                                                                                  >  2:705255/1‑143 (MQ=255)                 gCCTGGGTCAAAGGTGAACACACCTGTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAg                                                                                                                  <  1:1010940/143‑1 (MQ=255)                  ccTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGGTGAGGGTGAGTGTTGCTTTTGTATCAGCCATTTAAgg                                                                                                                 >  1:1024449/1‑143 (MQ=255)                   cTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCATTTCAGCACATCCAGCTCACCCGCTGTATCCCCGTTGCTGGCAGCTCTGGGTTTTGTATCCGCCATTTTAGGt                                                                                                                >  2:332039/1‑143 (MQ=255)                       gTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTcc                                                                                                            >  1:625495/1‑143 (MQ=255)                        tCACAGGTGAACACACCTTTTCAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCt                                                                                                           <  2:477541/143‑1 (MQ=255)                              gTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTTTATCCCCGTTGAGGGTTGGTTTTGCTTTTGTATCAGCCCTTTAGAGTCTCCTTAgcgc                                                                                                     >  2:607255/1‑143 (MQ=255)                                 aaCACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCtt                                                                                                  <  1:740494/143‑1 (MQ=255)                                      acCTTTTGAACCGAGAGTACGGATATCAATAATATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCCGCCTTTTAAGGTCTCCTTGGAGCTTTATTgg                                                                                             >  2:354482/1‑142 (MQ=255)                                        cTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCATTTAACATTGAGTTGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGt                                                                                           <  1:524296/143‑1 (MQ=255)                                         ttttGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTa                                                                                          <  2:777137/143‑1 (MQ=255)                                                     aggaCGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGaa                                                                              <  1:792644/140‑1 (MQ=255)                                                           gATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGTGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTaaa                                                                        <  2:203384/143‑1 (MQ=255)                                                              atCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAAttt                                                                     <  1:58228/143‑1 (MQ=255)                                                                 aaTAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTACATTTGcc                                                                  >  2:543795/1‑143 (MQ=255)                                                                  aTAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTACGTTTGCCt                                                                 >  1:182014/1‑143 (MQ=255)                                                                   tAACCTCTTGACCCAGCGTGCCTTTCAGCATATCCAGTTCAATAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCtt                                                                <  2:945304/143‑1 (MQ=255)                                                                     gcgccTTGGCCCAGCGTGCCTTTCGGCACCTCCCGTTCAACAGCGGTAGCCCCGTTGAGGGTGGGTTTTGCTTTTGTAGCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCg                                                              <  1:449393/139‑1 (MQ=255)                                                                             ccccGGCGTGCCTTTCAGCACATCTAGTTCAGCAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCaa                                                      <  2:945502/142‑1 (MQ=255)                                                                                   cGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTCTTGCTTTTGTATCAGCCATTTAAGTTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGc                                                <  1:1019704/143‑1 (MQ=255)                                                                                       ccTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTa                                            >  2:430073/1‑143 (MQ=255)                                                                                       ccTTTAAGCACATCCCGTTCCTCAGCTGTCTTCTCGTTGAGTGTGGTTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTACGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTa                                            <  2:89769/143‑1 (MQ=255)                                                                                          ttCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAACTTTGCCCTCGCACCTCCACCTGGCTTCAccc                                         >  2:452969/1‑143 (MQ=255)                                                                                             aGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTCCTTTTGTATCAGCCATTTTAGGTCTCCTCAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACATGATTTTAACAGtt                                      >  1:61228/1‑143 (MQ=255)                                                                                              gCACATCCAGTTCAACAGCTGTAGCCCCGTTGAGGGTGGGTTTTGCTGTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCGTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGttt                                     <  1:752216/143‑1 (MQ=255)                                                                                               cacaTCCAGTTCAACAGCTGTATCCCCGTAGAGGGTGGGGTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGtttc                                    >  2:406265/1‑142 (MQ=255)                                                                                                    ccAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCGCATCACCCGGGCTTTACCCCTTTCGTAtt                               >  2:941762/1‑143 (MQ=255)                                                                                                     cAGTTCAACAGCTGTATCCCCGTTTAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTAttt                              <  1:354482/143‑1 (MQ=255)                                                                                                     cAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTAttt                              <  1:705255/143‑1 (MQ=255)                                                                                                     cAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTAttt                              <  1:430073/143‑1 (MQ=255)                                                                                                     cAGTTCAACAGCCGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTCAGCGGCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTAttt                              <  1:581892/143‑1 (MQ=255)                                                                                                              aGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTcgccg                     >  1:249997/1‑143 (MQ=255)                                                                                                                cggaaaccccGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGct                   >  1:165192/7‑143 (MQ=255)                                                                                                                cTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGATTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTTGTTCGCCGct                   >  2:670443/1‑143 (MQ=255)                                                                                                                     tCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGACTCTTGTGTTGCTCGCCGCTCtggg              >  2:276221/1‑141 (MQ=255)                                                                                                                         cGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGaacg          >  1:809377/1‑141 (MQ=255)                                                                                                                          gTTGAGTGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAaga         <  1:104345/143‑1 (MQ=255)                                                                                                                          gTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAaga         <  1:101592/143‑1 (MQ=255)                                                                                                                           ttGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTCTACCCGTTTTTTATTTTGCTCGCTGCTCTGTGAAagac        >  2:976788/1‑142 (MQ=255)                                                                                                                           gtGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGTTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTTGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAagag        <  1:1095245/142‑1 (MQ=255)                                                                                                                             gAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCGTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAggg      <  1:137618/143‑1 (MQ=255)                                                                                                                             gAGGGTGGGTTTTCCTTTTGTATCAGCCATTTCAGGTCTCCTTAGCGCCTTAGTGCGTAAGACTGCCGGGACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAggg      <  2:940850/143‑1 (MQ=255)                                                                                                                              aGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGTCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAgggg     <  1:712724/143‑1 (MQ=255)                                                                                                                                 gTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGaaa  >  2:916831/1‑143 (MQ=255)                                                                                                                                    |                                                                                                                                            ATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAA  >  NC_000913/754315‑754585

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