Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	766,626	A→T	9.8%	S40C (AGC→TGC)	metB →	cystathionine gamma‑synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	766,626	0	A	T	9.8%	140.5 / 5.6	51	S40C (AGC→TGC)	metB	cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base A (23/23);  new base T (3/2);  total (26/25)
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.83e-01

ACAACTCAGCTGGGCGGGAACCGG‑ATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGAG  >  NC_002947/766483‑766757                                                                                                                                                 |                                                                                                                                    acaACTCAGCTGGGCGGGAACCGGAATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAgcgc                                                                                                                                 >  7:200487/1‑149 (MQ=255)   aaCTCAGCTGGGCGGGAACCGG‑ATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGTTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGcgccgc                                                                                                                              <  1:220162/149‑1 (MQ=255)        aGCTTGGCGTGAACCGG‑ATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGTGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTAt                                                                                                                          <  3:133797/148‑1 (MQ=255)             ggCGGGAACCGG‑ATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTa                                                                                                                     >  6:169782/1‑148 (MQ=255)                  gAACCGG‑ATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGa                                                                                                               <  6:140585/149‑1 (MQ=255)                     ccGG‑ATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACAt                                                                                                            >  4:235212/1‑149 (MQ=255)                              gAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATg                                                                                                       >  1:213645/1‑146 (MQ=255)                              gAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATg                                                                                                       <  2:213645/146‑1 (MQ=255)                               aacaacaacGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTg                                                                                                   >  2:279535/1‑149 (MQ=255)                                 caacaacGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGt                                                                                                 <  1:279535/149‑1 (MQ=255)                                        gAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGt                                                                                          >  2:266248/1‑149 (MQ=255)                                                 gAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAgg                                                                                  <  8:200487/148‑1 (MQ=255)                                                  aGGGCAATGCTGCTGCAGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGTCTACAAAGACATCGATGCCTGGTACGACGTGGCGCAggg                                                                                 >  6:276437/1‑148 (MQ=255)                                                    ggCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTg                                                                                         <  7:220710/138‑1 (MQ=255)                                                    ggCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTg                                                                                         >  8:220710/1‑138 (MQ=255)                                                    ggCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGc                                                                                >  7:297739/1‑147 (MQ=255)                                                                 ccGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTc                                                                 <  3:235212/149‑1 (MQ=255)                                                                     tgtgGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCt                                                             <  1:266248/149‑1 (MQ=255)                                                                        ggggACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTa                                                            >  6:72559/1‑147 (MQ=255)                                                                               cgGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGATCTCGATCGTGGTCAGCGCGGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCAt                                                   <  1:103478/149‑1 (MQ=255)                                                                                gggTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCAt                                                   >  4:286862/1‑148 (MQ=255)                                                                                 ggTCGTGTGGAGCGTGGACCAGGGACATCACCCGTAAAACGCCACGCAGACCCCGTTCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTTGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGa                                                 <  5:276437/149‑1 (MQ=255)                                                                                 ggTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTc                                                                 <  7:106353/133‑1 (MQ=255)                                                                                 ggTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTc                                                                 >  8:106353/1‑133 (MQ=255)                                                                                   tCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGc                                               <  5:169782/149‑1 (MQ=255)                                                                                         ggAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGa                                                 <  2:76721/141‑1 (MQ=255)                                                                                         ggAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGa                                                 >  1:76721/1‑141 (MQ=255)                                                                                              ggggAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGt                                    <  1:320946/149‑1 (MQ=255)                                                                                                   gCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGc                                                     <  1:118984/127‑1 (MQ=255)                                                                                                   gCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGc                                                     >  2:118984/1‑127 (MQ=255)                                                                                                      ggTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGa                                                 <  1:85164/128‑1 (MQ=255)                                                                                                      ggTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGa                                                 >  2:85164/1‑128 (MQ=255)                                                                                                           aGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTgg                         >  6:122037/1‑147 (MQ=255)                                                                                                                ccGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAg                                                                                   >  4:129515/1‑84 (MQ=255)                                                                                                                ccGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAg                                                                                   <  3:129515/84‑1 (MQ=255)                                                                                                                     cAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAg                                                         <  5:164933/105‑1 (MQ=255)                                                                                                                     cAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAg                                                         >  6:164933/1‑105 (MQ=255)                                                                                                                      aaCGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGc            <  5:122037/149‑1 (MQ=255)                                                                                                                         gCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCg                                                                       >  1:45691/1‑87 (MQ=255)                                                                                                                         gCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCg                                                                       <  2:45691/87‑1 (MQ=255)                                                                                                                          ccACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAg        >  5:250054/1‑149 (MQ=255)                                                                                                                           cACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGc       >  3:89083/1‑149 (MQ=255)                                                                                                                            gcGCTGACCACGATCGGGGTCTGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGaa                       >  6:213373/2‑132 (MQ=255)                                                                                                                            gcGCTGACCACGATCGGGGTCTGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGaa                       <  5:213373/131‑1 (MQ=255)                                                                                                                             cGCTGACCACGATCGGGGTCTGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCt                           <  7:206372/127‑1 (MQ=255)                                                                                                                             cGCTGACCACGATCGGGGTCTGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCt                           >  8:206372/1‑127 (MQ=255)                                                                                                                             cGCTGACCACGATCGGGGTCTGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTg     >  7:148877/1‑149 (MQ=255)                                                                                                                                aGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGAg  <  3:286862/149‑1 (MQ=255)                                                                                                                                aGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGAg  >  8:130295/1‑149 (MQ=255)                                                                                                                                        aTCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCAt                                                   <  6:176931/92‑1 (MQ=255)                                                                                                                                        aTCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCAt                                                   >  5:176931/1‑92 (MQ=255)                                                                                                                                                 |                                                                                                                                    ACAACTCAGCTGGGCGGGAACCGG‑ATATGAACAACAACGAAAAATTCGAGGGCAATGCTGCGGCCGGTGTGGGGACGCGGGTCGTGTGGAGCGGGGAGCAGGTACAGCACCCGTACAACGCCACGCAGACCCCGATCGTGGTCAGCGCCGCTTATGGCTACAACGACATCGATGCCTGGTACGACGTGGCGCAGGGCAAGCAGCCGGGCTTCATCTACAGCCGCATGAGCAACCCCACGGTCGCCACGCTGGAAGCAAAACTGGCCGAGCTGGAG  >  NC_002947/766483‑766757

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