Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	762,739	G→A	*406* (TAG→TAA)	sucB →	dihydrolipoyltranssuccinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	762,739	0	G	A	98.2%	170.3 / ‑4.1	55	*406* (TAG→TAA)	sucB	dihydrolipoyltranssuccinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (29/25);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

CCTGCCGATGATGTACCTGGCGCTGTCCTACGATCACCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAGTAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCA  >  NC_000913/762613‑762845                                                                                                                               |                                                                                                           ccTGCCGATGATGTACCTGGCGCTGTCCTACGATCACCGTCTGATCGATGGTCGCGAATCCGTGGGTTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGttt                                                                                                <  1:445525/139‑1 (MQ=255)       gatgatGTACCTGGCGCTGTCCTACGATCACCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTg                                                                                          <  1:1223242/139‑1 (MQ=255)         tgatgTACCTGGCGCTGTCCTACGATCACCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCa                                                                                        <  2:1232132/139‑1 (MQ=255)         tgatgTACCTGGCGCTGTCCTACGATCACCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCa                                                                                        <  2:248073/139‑1 (MQ=255)                     cgcTGTCCTACGATCACCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAGTAGTTTAAGTTTCACCTGCACTGTAGACCGGa                                                                            <  2:1132595/139‑1 (MQ=255)                       cTGTCCTACGATCACCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATa                                                                          <  2:829789/139‑1 (MQ=255)                                    aCCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCa                                                                                        <  1:894567/112‑1 (MQ=255)                                    aCCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGcc                                                             >  1:1137760/1‑139 (MQ=255)                                    aCCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGcc                                                             >  1:1148066/1‑139 (MQ=255)                                    aCCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATACGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCa                                                                                        >  2:894567/1‑112 (MQ=255)                                        tCTGATTGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTtct                                                         >  1:234329/1‑139 (MQ=255)                                           gATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCg                                                      >  2:88512/1‑139 (MQ=255)                                            aTCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCgg                                                     >  2:477756/1‑139 (MQ=255)                                            aTCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGACAGCTGCTGGACGTGCAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCgg                                                     >  2:1201870/1‑139 (MQ=255)                                                 tGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAAtt                                                >  1:1219666/1‑139 (MQ=255)                                                    tCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGa                                                                                 >  2:562943/1‑103 (MQ=255)                                                    tCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGa                                                                                 <  1:562943/103‑1 (MQ=255)                                                      gcgAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTc                                                                                               <  1:465779/87‑1 (MQ=255)                                                      gcgAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTc                                                                                               >  2:465779/1‑87 (MQ=255)                                                      gcgAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGc                                           >  2:543450/1‑139 (MQ=255)                                                      gcgAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGc                                           >  2:233101/1‑139 (MQ=255)                                                       cgAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTg                                                                                          >  2:650887/1‑91 (MQ=255)                                                       cgAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTg                                                                                          <  1:650887/91‑1 (MQ=255)                                                           tCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTTTAATAGTTTAAGTTTCAc                                                                                             <  2:730532/84‑1 (MQ=255)                                                           tCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCAc                                                                                             >  1:730532/1‑84 (MQ=255)                                                           tCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTg                                               <  1:723335/130‑1 (MQ=255)                                                           tCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTg                                               >  2:723335/1‑130 (MQ=255)                                                                      ccTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGAc                           >  2:474260/1‑139 (MQ=255)                                                                      ccTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGAc                           >  1:379314/1‑139 (MQ=255)                                                                            aaCGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGt                                                                                                  <  1:408957/62‑1 (MQ=255)                                                                            aaCGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGt                                                                                                  >  2:408957/1‑62 (MQ=255)                                                                              cGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTg                                                                                          <  1:781807/68‑1 (MQ=255)                                                                              cGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTg                                                                                          >  2:781807/1‑68 (MQ=255)                                                                                aTCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATc                 >  2:1260654/1‑139 (MQ=255)                                                                                 tCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAgg                                                                       >  2:171192/1‑84 (MQ=255)                                                                                 tCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAgg                                                                       <  1:171192/84‑1 (MQ=255)                                                                                  cAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAg               <  2:1219666/139‑1 (MQ=255)                                                                                          tGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACgg       <  1:88512/139‑1 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGa                                                                            >  2:152287/1‑69 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGa                                                                            <  1:152287/69‑1 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAGGCCTGATGCGACGCTGa                            >  2:169515/1‑117 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGa                            <  1:169515/117‑1 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACggg      >  2:587305/1‑139 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACggg      <  2:1148066/139‑1 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACggg      >  2:248817/1‑139 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACggg      <  2:234329/139‑1 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACggg      <  1:587305/139‑1 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACggg      <  1:248817/139‑1 (MQ=255)                                                                                           gCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGACTACggg      >  2:1274259/1‑139 (MQ=255)                                                                                                   aTCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGacca  >  2:260581/1‑135 (MQ=255)                                                                                                   aTCCGACGCGTCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGacca  <  1:260581/135‑1 (MQ=255)                                                                                                             tCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACggg      >  2:1075481/1‑121 (MQ=255)                                                                                                             tCTGCTGCTGGACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACggg      <  1:1075481/121‑1 (MQ=255)                                                                                                                       gACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCtt                    <  1:1072581/97‑1 (MQ=255)                                                                                                                       gACGTGTAATAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCtt                    >  2:1072581/1‑97 (MQ=255)                                                                                                                               |                                                                                                           CCTGCCGATGATGTACCTGGCGCTGTCCTACGATCACCGTCTGATCGATGGTCGCGAATCCGTGGGCTTCCTGGTAACGATCAAAGAGTTGCTGGAAGATCCGACGCGTCTGCTGCTGGACGTGTAGTAGTTTAAGTTTCACCTGCACTGTAGACCGGATAAGGCATTATCGCCTTCTCCGGCAATTGAAGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGGACCA  >  NC_000913/762613‑762845

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑