Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	766,609	A→T	8.7%	Q34L (CAG→CTG)	metB →	cystathionine gamma‑synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	766,609	0	A	T	8.7%	131.2 / 2.9	46	Q34L (CAG→CTG)	metB	cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base A (24/18);  new base T (2/2);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGA  >  NC_002947/766464‑766756                                                                                                                                                   |                                                                                                                                                     cGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGAACCGTATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGAc                                                                                                                                                   <  3:168725/149‑1 (MQ=255) cGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGAc                                                                                                                                                   >  6:107231/1‑149 (MQ=255)         cATCGGTGTACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATc                                                                                                                                            <  1:299014/148‑1 (MQ=255)         cATCGGTGTACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATc                                                                                                                                            <  1:303279/148‑1 (MQ=255)          aTCGGTGTACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGt                                                                                                                                          >  5:229716/1‑149 (MQ=255)              gtgtACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGt                                                                                                                                       >  4:107161/1‑148 (MQ=255)                gtACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCTGACCACGATCGGGGTctg                                                                                                                                    >  6:180195/1‑147 (MQ=255)                  acacAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAgcg                                                                                                                                  >  6:275872/1‑149 (MQ=255)                   cacaACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCTGACCACGATCGGGGTCTgcgt                                                                                                                                 <  5:180195/149‑2 (MQ=255)                   cacaACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAgcgc                                                                                                                                 <  5:107231/149‑1 (MQ=255)                    acaACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGtgcc                                                                                                                                <  1:38839/149‑1 (MQ=255)                      aaCTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACACGTACAACGCCACGCAGACGCCGATCGTGGTCAGcgccgc                                                                                                                              >  4:65603/1‑149 (MQ=255)                           aGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTAt                                                                                                                          <  6:276737/148‑1 (MQ=255)                            gCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATgg                                                                                                                        <  6:176393/149‑1 (MQ=255)                              tGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGc                                                                                                                       >  2:262676/1‑148 (MQ=255)                                ggCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCt                                                                                                                      <  8:217840/147‑1 (MQ=255)                                   gggAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGTGAGCAGGTACAGCAAACCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACaa                                                                                                                  <  1:85188/149‑1 (MQ=255)                                              atGAACAACAACGAAAAATTAGAGGGAAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGc                                                                                                     >  2:6887/1‑149 (MQ=255)                                                gAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCt                                                                                                   >  5:13391/1‑149 (MQ=255)                                                           aaaaaTTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTa                                                                                               <  7:253728/142‑1 (MQ=255)                                                           aaaaaTTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTa                                                                                               >  8:253728/1‑142 (MQ=255)                                                                ttCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGgcg                                                                                     >  2:315359/1‑147 (MQ=255)                                                                             ctgcGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGTGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTCGCGCAGGGCAAGCAGcc                                                                       <  1:6887/148‑1 (MQ=255)                                                                                  gCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAggg                                                                                >  1:241774/1‑134 (MQ=255)                                                                                  gCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAggg                                                                                <  2:241774/134‑1 (MQ=255)                                                                                   ccGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTc                                                                <  1:262676/149‑1 (MQ=255)                                                                                           gggACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAg                                                        >  2:262832/1‑149 (MQ=255)                                                                                            ggACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGGCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTGCATCTACAGc                                                       >  4:284819/1‑149 (MQ=255)                                                                                             gACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGcc                                                      <  6:229716/149‑1 (MQ=255)                                                                                                gcgGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑a                                                                                                            <  5:298573/92‑1 (MQ=255)                                                                                                gcgGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑a                                                                                                            >  6:298573/1‑92 (MQ=255)                                                                                                  gggTCGTGTGGAGCGGGGAGCAGGTACAGC‑TCCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCAt                                                  >  1:28121/1‑148 (MQ=255)                                                                                                  gggTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCAt                                                  >  4:113307/1‑148 (MQ=255)                                                                                                        tgtgGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAAc                                           <  3:107161/149‑1 (MQ=255)                                                                                                        tgtgGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAAc                                           >  1:62834/1‑149 (MQ=255)                                                                                                         gtgGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTc                                                                >  1:131255/1‑127 (MQ=255)                                                                                                            gAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCa                                       >  6:2557/1‑149 (MQ=255)                                                                                                            gAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCa                                       >  7:157016/1‑149 (MQ=255)                                                                                                              gCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCGCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACg                                     <  5:275872/149‑1 (MQ=255)                                                                                                               cGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACgg                                    >  6:124081/1‑149 (MQ=255)                                                                                                                               c‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTTCGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCCCGGTCGCCACGCTGGAAg                     >  6:67865/1‑148 (MQ=255)                                                                                                                                     gTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCCCGGTCGCCACGCTGGAAGCAAAAc               >  2:200274/1‑149 (MQ=255)                                                                                                                                             cggcGCTGACCACGATCGGGGTCTGCGCCGCTTATGGCTACAACGACGAACGATGCATGGTACGACGTGGCGCAGGCCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCACCCCCACGGACGCCa                              >  1:71499/4‑127 (MQ=255)                                                                                                                                             cggcGCTGACCACGATCGGGGTCTGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCa                              <  2:71499/123‑1 (MQ=255)                                                                                                                                               aCGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCt     >  2:108034/1‑149 (MQ=255)                                                                                                                                                  cAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGa  <  1:315359/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CGATGCTGCATCGGTGTACACAACTCAGCTGGGCGGGAACCGGATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGC‑ACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGAC‑ATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGA  >  NC_002947/766464‑766756

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