Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,485,644	C→T	G59R (GGG→AGG)	yhfA ←	OsmC family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,485,644	0	C	T	97.4%	130.4 / ‑4.3	38	G59R (GGG→AGG)	yhfA	OsmC family protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (15/22);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GTGCGTAAACAGGCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCCTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACACCATTTCCATCGGACTTGG  >  NC_000913/3485567‑3485719                                                                              |                                                                            gtgCGTAAACAGGCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCtttt                                                                          >  2:665281/1‑81 (MQ=255)      tAAACAGGCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTg                                                                       >  2:665293/1‑79 (MQ=255)        aaCAGGCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAg                                                                    >  2:665311/1‑80 (MQ=255)        aaCAGGCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAg                                                                    <  2:665246/80‑1 (MQ=255)          cAGGCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGAt                                                                 <  1:665247/81‑1 (MQ=255)           aGGCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATc                                                                <  2:665248/81‑1 (MQ=255)            ggCGTGGT‑CCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCg                                                               >  1:665303/1‑80 (MQ=255)             gCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCg                                                               >  2:665315/1‑80 (MQ=255)             gCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGa                                                              <  2:665250/81‑1 (MQ=255)             gCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGa                                                              <  2:665249/81‑1 (MQ=255)               gtggtgCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGaaa                                                            >  2:665299/1‑81 (MQ=255)                  gtgCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAAcc                                                          <  1:665252/80‑1 (MQ=255)                  gtgCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACca                                                         <  2:665251/81‑1 (MQ=255)                    gCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACcaca                                                       >  2:665294/1‑81 (MQ=255)                    gCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCATGACGCCTTTTTTGCAGGATCGAAACcaca                                                       >  1:665307/1‑81 (MQ=255)                     ccTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACAt                                                      <  2:665254/81‑1 (MQ=255)                     ccTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACAt                                                      <  2:665253/81‑1 (MQ=255)                          tCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATg                                                 >  2:665320/1‑81 (MQ=255)                           cgcgGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGa                                                   <  1:665255/78‑1 (MQ=255)                                  ttCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCa                                         <  1:665256/81‑1 (MQ=255)                                  ttCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCa                                         <  2:665257/81‑1 (MQ=255)                                       agGTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAaccac                                    <  1:665258/81‑1 (MQ=255)                                         gTCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCAccc                                  <  2:665259/81‑1 (MQ=255)                                          tCAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACccg                                 >  1:665295/1‑81 (MQ=255)                                           cAATTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACccgc                                >  1:665305/1‑81 (MQ=255)                                            aaTTTTACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACccgcc                               <  1:665260/81‑1 (MQ=255)                                              ttttACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCcgccgc                             <  2:665261/81‑1 (MQ=255)                                                ttACTTCACAATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCgccgcc                            <  2:665262/80‑1 (MQ=255)                                                      cacaATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCaa                       <  2:665263/79‑1 (MQ=255)                                                      cacaATCGACCACATCCTGACG‑CCTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAacac                    >  2:665301/1‑81 (MQ=255)                                                        caATCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACAcc                   <  2:665264/81‑1 (MQ=255)                                                           tCGACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACACCAtt                <  1:665265/81‑1 (MQ=255)                                                             gACCACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACACCATTTc              >  1:665309/1‑81 (MQ=255)                                                               ccACATCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACACCATTTCCa            >  1:665316/1‑81 (MQ=255)                                                                cacaTCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACACCATTTCCAt           >  1:665306/1‑81 (MQ=255)                                                                 acaTCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACACCATTTCCATc          >  1:665327/1‑81 (MQ=255)                                                                 acaTCCTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACACCATTTCCATc          <  1:665266/81‑1 (MQ=255)                                                                      cTGACGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACACCATTTCCATCGGACt     >  2:665304/1‑81 (MQ=255)                                                                         aCGCCTTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACACCATTTCCATCGGACTTgg  <  1:665267/81‑1 (MQ=255)                                                                              |                                                                            GTGCGTAAACAGGCGTGGTGCCTCTTCGCGGCGTTCAGAGGTCAATTTTACTTCACAATCGACCACATCCTGACGCCCTTTTTGCAGGATCGAAACCACATCGATGGCACTGCAACCACCCGCCGCCATCAACACCATTTCCATCGGACTTGG  >  NC_000913/3485567‑3485719

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