Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,442,298	G→A	intergenic (‑122/‑84)	msrA ← / → tamA	methionine sulfoxide reductase A/translocation and assembly module subunit TamA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,442,298	0	G	A	96.7%	99.4 / ‑4.0	30	intergenic (‑122/‑84)	msrA/tamA	methionine sulfoxide reductase A/translocation and assembly module subunit TamA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (19/10);  total (20/10)
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.91e-01

CAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAGGATATTCAG  >  NC_000913/4442221‑4442373                                                                              |                                                                            cAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACtt                                                                          >  2:918026/1‑81 (MQ=255)  aaGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACttt                                                                         >  1:918010/1‑81 (MQ=255)           aTTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCa                                                                 >  1:918028/1‑80 (MQ=255)            tttGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCata                                                               >  1:918041/1‑81 (MQ=255)                cacacaATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTAc                                                           >  2:918015/1‑81 (MQ=255)                      aTAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGAtt                                                     <  2:917983/81‑1 (MQ=255)                              cGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGtt                                               >  2:918032/1‑79 (MQ=255)                               gTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAg                                            >  1:918024/1‑81 (MQ=255)                                 gAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGcc                                          >  1:918030/1‑81 (MQ=255)                                  aaTAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGcc                                          <  1:917984/80‑1 (MQ=255)                                  aaTAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCg                                         >  1:918021/1‑81 (MQ=255)                                    tAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGtt                                       >  2:918017/1‑81 (MQ=255)                                      ggTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGttt                                      <  1:917985/80‑1 (MQ=255)                                            aaGTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCg                               >  1:918034/1‑81 (MQ=255)                                              gTAAAACATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGtt                             <  2:917986/81‑1 (MQ=255)                                                aaaaCATAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTc                            >  1:918035/1‑80 (MQ=255)                                                      tAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAgggg                      >  1:918038/1‑80 (MQ=255)                                                      tAAGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGc                     <  2:917987/81‑1 (MQ=255)                                                        aGAAATATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCa                    <  1:917988/80‑1 (MQ=255)                                                          aaaTATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAgg                  >  1:918027/1‑80 (MQ=255)                                                           aaTATTTGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGa                 >  2:918036/1‑80 (MQ=255)                                                                 tGGGGTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAg           >  1:918046/1‑80 (MQ=255)                                                                 tGGGGTTTAGCCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAg           >  2:918044/1‑80 (MQ=255)                                                                  ggggTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAgg          <  2:917989/80‑1 (MQ=255)                                                                    ggTTTAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAGGata       <  1:917990/81‑1 (MQ=255)                                                                      tttAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAGGATAtt     <  1:917991/81‑1 (MQ=255)                                                                       ttAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAGGATAtt     >  2:918029/1‑80 (MQ=255)                                                                       ttAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAGGATAtt     >  2:918033/1‑80 (MQ=255)                                                                        tAGTCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAGGATATTc    <  2:917992/80‑1 (MQ=255)                                                                           tCTACTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAGGATATTCAg  >  1:918040/1‑79 (MQ=255)                                                                              |                                                                            CAAGAGGCGTATTTGCACACAATAACCACCGTGAATAGGTCTAAAGTAAAACATAAGAAATATTTGGGGTTTAGTCTGCTTTTTAATCCATATTACTGGATTTTTGTTAAGCCGTTTAACGGCGTTCCAGGGGCAGGAAAAAAGGATATTCAG  >  NC_000913/4442221‑4442373

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