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.7%	121.7 / NA	45	L8P (CTC→CCC)	gltA	citrate synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (21/22);  total (23/22)
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

CAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGA  >  NC_000913/754311‑754574                                                                                                                                        |                                                                                                                                  cAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTCGAGGGTGGGTTTTGc                                                                                                                            >  2:494506/1‑143 (MQ=255) cAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCCGTCTCCCCGTCGAGGCTGCGTCTTGc                                                                                                                            >  2:384591/1‑143 (MQ=255)  aGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACTTCTTGACCAAGCTTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCt                                                                                                                           <  1:430032/143‑1 (MQ=255)     aTGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCtttt                                                                                                                        <  2:189013/143‑1 (MQ=255)     aTGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCtttt                                                                                                                        <  1:363019/143‑1 (MQ=255)          gTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGAGAGGTTTGATTTTGTAtt                                                                                                                   >  2:285527/1‑142 (MQ=255)                 tgaagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCAtt                                                                                                            <  1:816041/143‑1 (MQ=255)                 tgaagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCAtt                                                                                                            <  2:319664/143‑1 (MQ=255)                   aagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTa                                                                                                          <  1:457740/143‑1 (MQ=255)                    agCCTGGGTCAAAGGTGAACACACCTTGTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCGAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTaa                                                                                                         <  2:622487/143‑1 (MQ=255)                     gCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAg                                                                                                        <  1:499073/143‑1 (MQ=255)                            tCAAAGGTGAACACACCTTTTGAAACGAGAGTACTGATATCAATAACATCTTTATCCATCTTGCCTTTCAGCACATCCAGTTCAAAAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCt                                                                                                 <  1:60677/143‑1 (MQ=255)                                  gTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACACCCAGTTCAACAGCTGTATCCCCGTCCAGGGTGGGTTTTGCTTCTGTACCAGCCATTTAAGGCCTCCTTacccc                                                                                           >  2:622779/1‑139 (MQ=255)                                        tcacCTCTTGATCCGAGAGTACGGATATCAATAACATCTTGCCCCAGCGTGCCTTTCAGCACATCCAGTTCAGCAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTTTCTGCCATTTAAGGTCTCCTTAGCGCCTTAtt                                                                                     <  2:512502/142‑1 (MQ=255)                                           ccTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTAGCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCg                                                                                  <  1:554317/143‑1 (MQ=255)                                            cTTTTAAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGt                                                                                 <  1:704888/143‑1 (MQ=255)                                                       agagTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGCGGGTTTTGCTTTTGTATCAGCCAATTAAGGTCTCCTTAGGGCCTTATTGCGTAAGACTACCgg                                                                      >  2:760909/1‑143 (MQ=255)                                                           tACGGATATCAATAACATCTTTACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCTTTTAAGGTATCCTTAGCGACTTATTGCGGAATACTGCCgagact                                                                  >  2:36967/1‑138 (MQ=255)                                                                  atCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCACCAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTCGCGCCTTATTGCGTAATACTGCCGG‑AAATTAAAttt                                                          >  1:257870/1‑143 (MQ=255)                                                                  atCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCGTTATTGCGTAAGACTGCCGGTAACTAAAAta                                                           >  2:39093/1‑142 (MQ=255)                                                                     aaTAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGcc                                                       >  1:624298/1‑143 (MQ=255)                                                                       tAACATCTTTACCCAGCTTGCCCTTCAACACATCCAGTCCTACAGCTGTTCCCCCTTTGAGTGTGGTTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCCCCTTATTGCGTAAGACTGCCGT‑AACTTAAATTTGCCtt                                                     <  2:154361/143‑1 (MQ=255)                                                                       tAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTAATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCtt                                                     >  1:3019/1‑143 (MQ=255)                                                                        aaCATCTTGTCCCTGCTTTCTTTTCATTATTTCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTc                                                    <  1:279994/143‑1 (MQ=255)                                                                             cTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCtttt                                                                                                                        <  1:312981/71‑1 (MQ=255)                                                                             cTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCtttt                                                                                                                        >  2:312981/1‑71 (MQ=255)                                                                             cTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGCGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCCTCTTTCGTAACACTGCCGG‑AGCCTACACTTTCCTTCCcaca                                               >  2:621322/1‑143 (MQ=255)                                                                                    cAGCGTGGCTTTGAGCACATGCAGTTCAACGGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCCTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCt                                        <  2:692594/143‑1 (MQ=255)                                                                                      acGTGCCTTTCAGCACATCTATTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTgg                                      <  2:103297/142‑1 (MQ=255)                                                                                        gTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAAttt                                                          <  1:672308/121‑1 (MQ=255)                                                                                        gTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAAttt                                                          >  2:672308/1‑121 (MQ=255)                                                                                          gCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTCTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTTCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCttt                                  >  1:230096/1‑143 (MQ=255)                                                                                             tttCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑CACTTAACTTTGCCTTCGCACTTCCTCCTGGCTTTAcc                               >  1:32406/1‑143 (MQ=255)                                                                                               tCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCg                             >  1:721002/1‑143 (MQ=255)                                                                                               tCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTAACCg                             >  1:720302/1‑143 (MQ=255)                                                                                                      aTCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTa                      >  2:348421/1‑143 (MQ=255)                                                                                                         cAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCTCCTTATTCCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTGTACCCGTTTTTTAttt                   >  1:693855/1‑143 (MQ=255)                                                                                                           tgtCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTgg                 <  2:444128/141‑1 (MQ=255)                                                                                                           gTTCATCAGGTGTATCCCCGTTGTGGGTGGGTTTTGCTTTTGTAGCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTgg                 <  2:444954/143‑1 (MQ=255)                                                                                                            ttCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTAtttatc                >  2:61394/1‑140 (MQ=255)                                                                                                               aaCAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTCCCCGTTTTTTCTTTGGCttg             >  2:89707/1‑141 (MQ=255)                                                                                                                aCAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCCGGCTTTACCAGATTTTTATTTgccaagc            >  2:271069/1‑137 (MQ=255)                                                                                                                  aGCTGTATCCCTGTTGAGGGTGGGTTTTGCTTTTGTATCAGCTATTTAAGGTCTCTTTAGCGCCTTATTGCGTAATACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTcgccg          <  1:453098/143‑1 (MQ=255)                                                                                                                    cTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACGCGTTTTTTATTTGGCTCGCCGct        >  2:780707/1‑143 (MQ=255)                                                                                                                          ccccGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGa  <  2:624298/143‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  CAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGG‑AACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGA  >  NC_000913/754311‑754574

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