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%	114.1 / ‑7.0	41	L8P (CTC→CCC)	gltA	citrate synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base G (15/24);  minor base T (1/0);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCGCAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAA  >  NC_000913/754306‑754585                                                                                                                                               |                                                                                                                                            aTTCGCAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACAGCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCGGTATCCGCGTAG‑Agggtgggt                                                                                                                                           >  1:458413/1‑143 (MQ=255)    cGCAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATAGATACCATCTTGACCAAGTCTGCCCGACAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGTGtttt                                                                                                                                        >  1:383994/1‑143 (MQ=255)               gTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATc                                                                                                                             >  2:536997/1‑143 (MQ=255)                ttGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCa                                                                                                                            >  1:538577/1‑143 (MQ=255)                     gtgaagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCAt                                                                                                                       <  1:734706/143‑1 (MQ=255)                       gaagCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAACAACATCTTGACCCAGCGTGCCTTTCA‑CAAATCATTTTGTATCGGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCAttt                                                                                                                     >  1:666318/1‑143 (MQ=255)                          gCCTGGGTCAAACGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAg                                                                                                                  <  1:884073/143‑1 (MQ=255)                                   aaaGGTGAACACACCTTTTGCACCGAGAGTACGGATAGCAAGAACATCTTGACCCAGCGTGCCTTTCAGCAGAGCTAGTT‑CGGCAGCTGTAGCCCCGTTG‑AGGGGGGGTTTTGCTGTGGTATCAGCCATTTAAGGTCTCCTTa                                                                                                         <  1:652012/143‑1 (MQ=255)                                   aaaCGTCAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTa                                                                                                         <  1:197467/143‑1 (MQ=255)                                    aaGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAg                                                                                                        <  2:198018/143‑1 (MQ=255)                                         gAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCCGCCATTTAAGGTCTCCTTACCGCCt                                                                                                   >  1:32072/1‑143 (MQ=255)                                            cacacCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCGTTAt                                                                                                >  1:801477/1‑143 (MQ=255)                                               acCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCCCCGTCCCTTTCAGCACATCCAGTT‑CAACAGCTGTATACCCGTTG‑ACGGCCGACCTTGCTCTTGTAGCAGCCATGTCAGGTCCCCCTCGCGGCTTATTGc                                                                                             >  2:190500/1‑143 (MQ=255)                                               acCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCGTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGATTTGGCTTTTGTATCAGCCATTTCAGGTCTCCTTAGCACCTTAttac                                                                                             >  2:19767/1‑141 (MQ=255)                                                 cTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGt                                                                                           <  1:517906/143‑1 (MQ=255)                                                 cTTTTGAACCGAGAGAACGGATATCAACAACATCTTGACCCAGCGTGCCTTTCAGCAAATCCAGTT‑CAATAGCTGTATCCCCGTTG‑ATGGTGGTTTTTGCTTTTGTATCAGCTATTTTAGGTCTCCTTAGCGCCTTATTGCGt                                                                                           <  1:830067/143‑1 (MQ=255)                                                         ccGCGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGc                                                                                   <  1:274269/143‑1 (MQ=255)                                                         ccGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGc                                                                                   <  1:536997/143‑1 (MQ=255)                                                         ccGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCACCGTTA‑CGGCTGGGTTTTGCTTTTGGATCCGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGc                                                                                   >  1:570546/1‑143 (MQ=255)                                                                   ggATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑CGGGTGGGTTTTGCTTCTCTACCCGCCCTTTCCGGTCTCCTTAGCGCCTTCTTGCGTAGGAATGCCGTAACTTaa                                                                         >  1:61000/1‑143 (MQ=255)                                                                      cagcAATACCATCTTGCCCCGGCGCGCCTCCCAGAACATCGTGTT‑CAGCAGCTGTTTCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCGCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAAtt                                                                      <  1:31843/140‑1 (MQ=255)                                                                       atCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGGTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAAttt                                                                     >  1:598069/1‑143 (MQ=255)                                                                           aTAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGAATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTGAAATTTGCCt                                                                 >  2:277622/1‑143 (MQ=255)                                                                              aCATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCg                                                              <  2:258512/143‑1 (MQ=255)                                                                                 tCTTGCCCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTCTACCAGGTG‑ACGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGcac                                                           <  1:207520/143‑1 (MQ=255)                                                                                 tCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGcac                                                           <  2:538577/143‑1 (MQ=255)                                                                                              tGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCtt                                              <  1:325112/143‑1 (MQ=255)                                                                                                 cTTTGAGCGCATCCAGTG‑CCACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGTTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTAc                                           <  2:24335/143‑1 (MQ=255)                                                                                                        cacaTCCAGTT‑CAACAGCTGTATCCCCGTTGTCAGGATATGAGTGC‑ATTGTATCAGCCATTTAAGGTCTCCTTATCTCCTGATGGCGCAAGACTGCCGGGACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGtgtt                                    >  2:155108/1‑143 (MQ=255)                                                                                                         acaTCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGAGCCTTATTGCGTAGGACTGCGGGAACTTAAATTTGCCTTCGCACATCAACCAGGGCTTACCCGtgtgt                                   >  2:331817/1‑139 (MQ=255)                                                                                                             ccAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTAtt                               <  2:328526/143‑1 (MQ=255)                                                                                                              cAGTT‑CAACAGCTGTATCCCCGTTT‑CGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTAttt                              <  2:262893/143‑1 (MQ=255)                                                                                                                   ttaaCAGCTGTATCCCCGTTG‑AGTGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCt                          <  1:796939/141‑1 (MQ=255)                                                                                                                           tGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTa                                                                                                         >  1:823764/1‑56 (MQ=255)                                                                                                                           tGTATCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTa                                                                                                         <  2:823764/56‑1 (MQ=255)                                                                                                                              agctccGTTG‑AGGGTGGGTTTGGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTGATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCACCCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCtgt               <  2:895676/139‑1 (MQ=255)                                                                                                                              aTCCCCGTTG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCtgt               <  2:796041/143‑1 (MQ=255)                                                                                                                                      tG‑AGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAgg       <  1:659409/143‑1 (MQ=255)                                                                                                                                      tG‑AGGGGGGGTTTTGCTTTTGTATCAGCCAGTTGAGGTCGCCTGAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAgg       <  1:369917/143‑1 (MQ=255)                                                                                                                                         aGGGTGGGGTTTGCTTTTGTATCAGCCGTTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAgggg     <  2:911209/143‑1 (MQ=255)                                                                                                                                          gggtgggtTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGa    >  1:471013/1‑143 (MQ=255)                                                                                                                                            gTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTTTGTGAAAGAGGGGaaa  >  1:729865/1‑143 (MQ=255)                                                                                                                                            gTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCGCTGTGAAAGAGGGGgaa  >  1:742880/1‑143 (MQ=255)                                                                                                                                               |                                                                                                                                            ATTCGCAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTT‑CAACAGCTGTATCCCCGTTG‑AGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAA  >  NC_000913/754306‑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: ‑