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	94.7%	80.5 / ‑4.8	38	L8P (CTC→CCC)	gltA	citrate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/16);  minor base T (0/2);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

CGCAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAA  >  NC_000913/754309‑754584                                                                                                                                          |                                                                                                                                           cGCAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGTACATCCAGTTAAACAGCTGTATACCAGTTGAGGGTGGGttat                                                                                                                                       >  2:39984/1‑141 (MQ=255)       aTGCGGTTGAAGTGAAGCCTGTGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCtttt                                                                                                                                 <  1:13070/143‑1 (MQ=255)             ttGAAGTGAAGCCTTGGTCGAAGGTGAACACACCTTTTGAACCGAGAGTCCGGATATCTATAACAGCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCa                                                                                                                           <  2:88615/143‑1 (MQ=255)             ttGAAGTGAAGCCTGTGTCAAAGGTGATCACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCa                                                                                                                           <  1:383503/143‑1 (MQ=255)              ttaagtgaagCCTGGGGCAGAGGTGAACACAGCTGTTGGACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAg                                                                                                                          <  2:30631/141‑1 (MQ=255)              tgaagtgaagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAg                                                                                                                          >  1:409200/1‑143 (MQ=255)                   tgaagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTATCAGCACATCCAGCTCAAGAACTGTATCGCCGCTACTGCTGGGTCTTGCTCTTGTATCatcgctc                                                                                                                     >  2:719319/1‑137 (MQ=255)                       gCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACCTCCAGTTCAGCAGCTGTATTCCCAATGAAAGTGGATATAAAATTTGTATCAGACATTTAAg                                                                                                                 >  1:15104/1‑143 (MQ=255)                                          acacCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTAtt                                                                                              <  1:743764/143‑1 (MQ=255)                                          acacCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCATTAACGACTGAtt                                                                                              >  1:320729/1‑143 (MQ=255)                                          acacCTTTTGAACCGAGAGGACGGATATCAATAACATCTTGACCCAGTGTGCTTTTTAGCACATCTAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTAtt                                                                                              <  1:340048/143‑1 (MQ=255)                                             ccTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGTGGGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCg                                                                                           <  2:552540/143‑1 (MQ=255)                                                      ccGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTGGTATCAGCCATTTCAGGTCTCCTTCGCGCCTTATTGCGTAAAACTga                                                                                  >  1:31716/1‑142 (MQ=255)                                                       tgagaGTACGGATATCAATATCATGTTGACCCAGCGTGCATTTCAGCACATCGAGTTCAACAGCTGTAGCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGcc                                                                                 <  1:765351/142‑1 (MQ=255)                                                           agTACGGATATTAATATCATCTTGACCCAGCGTGCCTGTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGaa                                                                             <  1:4267/143‑1 (MQ=255)                                                                ggATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGATCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTaa                                                                        >  1:60556/1‑143 (MQ=255)                                                                ggATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCATATATATAAGACATTTAAGGTCTCATTAGAGCCTTATTGCGTAAGAGTACCGGAACTTaa                                                                        >  1:600597/1‑143 (MQ=255)                                                                   tatCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGATCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAAtt                                                                     >  2:535725/1‑143 (MQ=255)                                                                    atCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGGCTCCTTAGCGCCTTATTGACTAATACTGCCGGAACTCAAAtca                                                                    >  2:370546/1‑141 (MQ=255)                                                                       aaTAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGcc                                                                 >  2:89382/1‑143 (MQ=255)                                                                           aCATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCg                                                             <  1:71389/143‑1 (MQ=255)                                                                            cATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGCACTTAAATTTGCCTTCGc                                                            >  1:528365/1‑143 (MQ=255)                                                                                    cccAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCCTTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCTGAACTTAAATTTGCCTTCGCACaacaac                                                    >  2:364659/1‑143 (MQ=255)                                                                                        gCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTgg                                                <  2:527536/143‑1 (MQ=255)                                                                                             gcTTGCAGCAGATCCAGTTCAGCAGCTGTATCCTCGTTGATGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTa                                           <  2:771424/142‑1 (MQ=255)                                                                                              ctctgtgCACATCCAGTTCATAAGCTGTCTCCACGTTCACGTTGTGTCTTGCTTTTGTATCAGCCATTTAAAGTATTCTTCCCGCCTTAGTGCGTAAGACTGCCGGAACTTATATTTGCCTTCGCACATCAACCTGGCTTTAc                                          <  2:254427/137‑1 (MQ=255)                                                                                              cTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCAATTTACGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCAGGAACTTAAATTTGCCTTCCCACATCACCCTGGCatttc                                          >  2:265566/1‑141 (MQ=255)                                                                                                       caTCCAGTTCAATAGCTGTATCCCCGTTTAGGCTGTTTTTTGCTTTTGTATCATCTATTTAAGTTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGtttttt                                 <  1:754057/143‑1 (MQ=255)                                                                                                       caTCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGtttttt                                 <  1:305207/143‑1 (MQ=255)                                                                                                        aTCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTGGCTGTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTa                                >  2:785950/1‑143 (MQ=255)                                                                                                         tCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCGCATCAACCTGGCTTTACCCGTTTTTTAt                               >  1:79296/1‑143 (MQ=255)                                                                                                                  aCAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTcgc                      <  1:762100/143‑1 (MQ=255)                                                                                                                  aCAGCTGTATCCCCGTTGAACGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGACTTATTGCATAAGACTGCCGGAACTTAAATTTGCCTTCGCACACCAACCCAGCTCTACCCCTTTGTGATTTTGcgcgt                      >  2:122463/1‑142 (MQ=255)                                                                                                                    aGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTATCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTGTACCCGAATTTTATTTGTCTctccg                    >  2:183681/1‑143 (MQ=255)                                                                                                                      cTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGGCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGGACTTAAATTTGCCTTCGCACATCACCCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGct                  <  2:87088/143‑1 (MQ=255)                                                                                                                       tGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGctc                 <  1:107427/143‑1 (MQ=255)                                                                                                                                     gggtgggtTTTGCTTTTTTATCAGTCATTTAAGGTCTCCTTAGCGCCTTAAAGCGGAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGa   <  1:98100/143‑1 (MQ=255)                                                                                                                                      ggTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCCCATCCACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGaa  >  1:787344/1‑143 (MQ=255)                                                                                                                                      ggTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGAACATAGCCCGTGCGTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGaa  >  1:579204/1‑143 (MQ=255)                                                                                                                                          |                                                                                                                                           CGCAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAA  >  NC_000913/754309‑754584

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