Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	760,485	G→A	G594S (GGC→AGC)	sucA →	subunit of E1(0) component of 2‑oxoglutarate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	760,485	0	G	A	92.9%	97.8 / ‑6.5	44	G594S (GGC→AGC)	sucA	subunit of E1(0) component of 2‑oxoglutarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (23/16);  total (26/18)
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

AAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTACA  >  NC_000913/760348‑760623                                                                                                                                          |                                                                                                                                           aaGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGt                                                                                                                                       <  1:53267/143‑1 (MQ=255)           gAAGCGCCTGCAGGAGCTGGCTAAACGCATCAGGACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAAcc                                                                                                                             <  2:129651/143‑1 (MQ=255)                  cTGCAGGAGCTGGCGACACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCACGGTTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAGACTGTTCGCCTGGAGCGGTGCGGAACACCTCGCGTa                                                                                                                      >  1:598867/1‑143 (MQ=255)                        gAGCTGGCGAAACGCATCAGCACGGAGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATAGCGATCACCAGGCGATGCCTGCCGGTGAGAAACTGTTCGACAGGAGCGGTGCGGCAAACCCCGCTTACGCCAc                                                                                                                >  2:275950/1‑143 (MQ=255)                         aGCTGGCGAAACGCAACAGCACGGTGGCGGAAGCAGCTAAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCCATGGCGACCGGTGAGACACTGTTCGACTGGAGCGGTGCGGAAAAACGCGAGTACGCCACg                                                                                                               >  2:687217/1‑143 (MQ=255)                            tGGCTAAACGCACCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTg                                                                                                            <  2:511448/143‑1 (MQ=255)                              gCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAAATGTTTGGCTGGGGCGGTGCGGAAAAACTCGCGTAGGCCACgacggt                                                                                                          >  2:126232/1‑138 (MQ=255)                                  aaCGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGTGAAACTGTTCGACTGGCGCGGTGCGGAAGACCTCGCTTACGCCACGCTAGTtgat                                                                                                      >  1:688822/1‑143 (MQ=255)                                    cGCATCAGCACGGTGCCGGAAGCAGTTGACATGCAGTCTCGCGTTCCCAAGATTTATGGCGATCGCCAGGCGATCGCTGCCGGGGAGAAACTGTTCGACTGGGGCGGGGCGGAGCACCCCGCTTCCGCCGCGCGGGtgcgtgc                                                                                                    >  1:335226/1‑137 (MQ=255)                                     gCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGTAGCGGTGCGGAAAAAATCGCTTACGCCACGCAGGTtgatgac                                                                                                   >  1:593803/1‑142 (MQ=255)                                      cATCAGCACGGTGCAGGAAGCAGTTGAAATGCAGTCTGGTGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTTAGAAACTGTTCGACTGGAGAGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAg                                                                                                  <  2:392794/143‑1 (MQ=255)                                         cAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCa                                                                                               >  1:490745/1‑143 (MQ=255)                                          aGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAACCCTCGCTTACGCCACGCTGGTTGATGAAggtct                                                                                              >  1:183700/1‑140 (MQ=255)                                           gCACGGTGCCGGAAGCAGTTGACATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCCACAGGCGGTGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCCTACGCCACCCCGCCTCATGAAAGCAtt                                                                                             >  2:556918/1‑143 (MQ=255)                                           gCACGGTGCCGGAAGCAGTTGAAATGCAATCTCGCGTAACCACGATGTATGGCGATCGCCAAGCGATGGCTGCCGGTGAGAAACTGTTCGACTGAAGCGGTACGGAAAACCTCGCTTACGCCACGCTGGTTGACGCAggactg                                                                                             >  2:512449/1‑139 (MQ=255)                                              gggTGCTGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGGGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCg                                                                                          <  1:766436/142‑1 (MQ=255)                                                     ggAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGcc                                                                                   >  1:745280/1‑143 (MQ=255)                                                      gAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATAGCGATCGCCAGGCGCTGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGACAACCTCGCTTACGCCACGCTGGTTGATGACGGCAGTCCGGGTCGccg                                                                                  >  2:847284/1‑142 (MQ=255)                                                          cAGTTGAAACGCAGGCCCGCCTTGCCAAGAGTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCg                                                                              <  1:150381/143‑1 (MQ=255)                                                           aGTTGAAATGCAGTCTCGCGTTGCCAAGACTTAGGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACATCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCgg                                                                             >  1:626028/1‑143 (MQ=255)                                                            gTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAACCACCGCTTACGCCACGCGGGTTGATGAAGCGCCACCGGTTCGACTGGCggg                                                                            >  2:617162/1‑143 (MQ=255)                                                               gAAATGCAGTCTCGCGTTGCCAAGACTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGa                                                                         <  2:724144/143‑1 (MQ=255)                                                                aaaTGCAGTCTCGCGTTTCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGaa                                                                        <  1:751423/143‑1 (MQ=255)                                                                        tctcGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGCTGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCgg                                                                >  2:452601/1‑143 (MQ=255)                                                                        tctcGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCgg                                                                >  1:112412/1‑143 (MQ=255)                                                                            gggTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTTGCCTGTCGGGTGAAGACTCCGGTcgc                                                            <  2:432619/141‑1 (MQ=255)                                                                            gcgTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTcgc                                                            >  2:811520/1‑143 (MQ=255)                                                                             ggTTGCCAAGATTTATGGCGGTCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTcgcg                                                           <  2:299741/142‑1 (MQ=255)                                                                             cgTTGCCAAGTTTTATGGCGTTCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTcgcg                                                           <  1:217312/143‑1 (MQ=255)                                                                                           gtgtggttCGCCAGGCGATGGCTGCCGGTGAGGAATTGTTCGACTGGGGCGGTGCTGAAAACCTCGCTTACGTCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCAc                                             <  2:342301/136‑1 (MQ=255)                                                                                             ggCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGGCTGTCGGGTTAACACTCCGGTCGCGGTACCTTCTTCCACCg                                           <  2:697007/143‑1 (MQ=255)                                                                                                 aTCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGGCACGCCGGCTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCAc                                       >  1:313097/1‑143 (MQ=255)                                                                                                    gCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGTTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTTTCAGGTGAAAACTCACGTCTCGGTACTTTTTTTCAACCCCAcgca                                    >  1:521514/1‑142 (MQ=255)                                                                                                    gCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGCCGGCATTCCGGTTCGCCTGTCCGGTGAATACTCCGGTCGCGATACCTTCTTCAACCGCCAcgct                                    >  1:639201/1‑142 (MQ=255)                                                                                                    gCCAGGCGATGGCTGCCGGCAAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGACGACTCCGGTCGCGGTGCCTTCTTCCACCGCCAcgcc                                    >  2:181792/1‑142 (MQ=255)                                                                                                     ccAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCgg                                   >  1:683097/1‑143 (MQ=255)                                                                                                     ccAGGCGATGGCTGCCGGAGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATAACGGTTCGCCTGGCGGGTGAAGACTCCGGTCGCGAGACCATCATCCACCGCCAcggcg                                   >  2:42302/1‑140 (MQ=255)                                                                                                       acgCGATGGCTGCCGGTGAGAAACTGTTAGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCCTTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGGCCACTCTCGCGCCGCCACGCGGTg                                 >  2:344182/3‑143 (MQ=255)                                                                                                         gCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGAt                               <  2:112412/143‑1 (MQ=255)                                                                                                                cTGTCGTTGAGAATCTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTCGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAAc                        <  2:278615/143‑1 (MQ=255)                                                                                                                  gCCGGTGACATACTGTTCGACTGGAGCGGTGCGGAACACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCACCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCa                      <  2:216082/143‑1 (MQ=255)                                                                                                                              cTGTTCTACTGGTGCGTTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATGCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTc          <  1:716970/143‑1 (MQ=255)                                                                                                                                  tCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATGCTCAACCCGTCTGACTGTTTCatt      >  1:125451/1‑141 (MQ=255)                                                                                                                                      cTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGTATTCCGGTTCGCCTGTCGGGCGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTaca  <  1:121213/143‑1 (MQ=255)                                                                                                                                          |                                                                                                                                           AAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTACA  >  NC_000913/760348‑760623

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