Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,061,202	A→G	N118S (AAC→AGC)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,061,202	0	A	G	97.4%	134.5 / NA	39	N118S (AAC→AGC)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (12/26);  minor base T (0/1);  total (12/27)
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

TCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGT  >  NC_000913/3061073‑3061337                                                                                                                                  |                                                                                                                                        tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                                >  2:229387/1‑139 (MQ=255) tCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGt                                                                                                                                >  1:667927/1‑139 (MQ=255)         tATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTggcgggcg                                                                                                                        >  2:261575/1‑139 (MQ=255)          attaagatttttCAACAGCAAAAAAGTGATACGGTTGTGATCGCCGTCGCCTCGCCGCCGTGCAAAAAGGTCTTACCGTTACGTTTGACCTCGCCACTCGCCGTGGCTACGACTCCGATATGCCGCGCGTGGCGGGCGa                                                                                                                       <  2:202702/132‑1 (MQ=255)                   tttCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCaa                                                                                                              >  2:102568/1‑139 (MQ=255)                         cAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCggg                                                                                                        <  1:261575/139‑1 (MQ=255)                           gCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                      <  2:353182/139‑1 (MQ=255)                           gCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                      <  1:459028/139‑1 (MQ=255)                           gCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCg                                                                                                      <  2:627716/139‑1 (MQ=255)                                       aaCGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGaca                                                                                          <  1:442140/139‑1 (MQ=255)                                       aaCGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGaca                                                                                          <  2:751127/139‑1 (MQ=255)                                       aaCGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGaca                                                                                          <  2:212315/139‑1 (MQ=255)                                         cGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGcc                                                                                                                                       <  2:382366/92‑1 (MQ=255)                                            tttttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                     <  2:417369/139‑1 (MQ=255)                                            tttttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTg                                                                                     <  1:542394/139‑1 (MQ=255)                                              tttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGa                                                                                   <  1:196499/139‑1 (MQ=255)                                               ttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGaa                                                                                  <  2:416583/139‑1 (MQ=255)                                                        tAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGaaa                                                                         >  2:538237/1‑139 (MQ=255)                                                          aCCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGt                                                                       <  2:63468/139‑1 (MQ=255)                                                           ccTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTc                                                                      <  2:199791/139‑1 (MQ=255)                                                               gccgccGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGt                                                                  >  1:607344/1‑139 (MQ=255)                                                                ccgccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGtt                                                                 <  1:593755/139‑1 (MQ=255)                                                                 cgccgGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTc                                                                <  1:322244/139‑1 (MQ=255)                                                                        cAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGa                                                         <  2:607344/139‑1 (MQ=255)                                                                             aGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCg                                                    <  1:102568/139‑1 (MQ=255)                                                                             aGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCg                                                    <  1:229387/139‑1 (MQ=255)                                                                                       gTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATaaaa                                          <  2:771912/139‑1 (MQ=255)                                                                                         tGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATg                                        >  2:741802/1‑139 (MQ=255)                                                                                                gACCTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGttt                                 >  2:211417/1‑139 (MQ=255)                                                                                                   cTTGCCACCCACCGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGa                              <  2:428434/139‑1 (MQ=255)                                                                                                               cGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGtt                                                                 <  2:383891/92‑1 (MQ=255)                                                                                                               cGTGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGtt                                                                 >  1:383891/1‑92 (MQ=255)                                                                                                                 tGGCTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgc                <  1:741802/139‑1 (MQ=255)                                                                                                                    cTACGACTCCGATAGCCCGCGCGTGGCGGGTGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGTCCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCa             <  2:237447/139‑1 (MQ=255)                                                                                                                    cTACGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCa             <  1:547861/139‑1 (MQ=255)                                                                                                                      aCGACTCCGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGt           <  1:538237/139‑1 (MQ=255)                                                                                                                             cGATAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCa    >  2:410248/1‑139 (MQ=255)                                                                                                                               aTAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGt  >  1:610196/1‑139 (MQ=255)                                                                                                                               aTAGCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGt  >  2:660140/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                        TCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGT  >  NC_000913/3061073‑3061337

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