Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,865,065	C→A	A403A (GCG→GCT)	yidE ←	putative transport protein YidE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,865,065	0	C	A	89.5%	52.5 / 0.7	19	A403A (GCG→GCT)	yidE	putative transport protein YidE
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (8/9);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGACCGAGAGGAACAGCACGATCCCCAGCTCCCGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCCGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATGCCAATAAACACTGGCAGCATCTGAACCTGTTGCAGTTTTTGCTGCGCATTCCCCAGC  >  NC_000913/3864936‑3865200                                                                                                                                  |                                                                                                                                        aaCGACCGAGAGGAACAGCACGATCCCCAGCTCCCGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCa                                                                                                                               >  1:231557/1‑140 (MQ=255) aaCGACCGAGAGGAACAGCACGATCCCCAGCTCCCGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCa                                                                                                                               >  2:41567/1‑140 (MQ=255)   cGACCGAGAGGAACAGCACGATCCCCAGCTCCCGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGt                                                                                                                             >  2:196147/1‑140 (MQ=255)          gagGAACAGCACGATCCCCAGCTCCCGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAAcg                                                                                                                      <  1:273950/140‑1 (MQ=255)                        tCCCCAGCTCCCGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCt                                                                                                        >  2:305691/1‑140 (MQ=255)                           ccAGCTCCCGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCCGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGc                                                                                                     >  1:79147/1‑140 (MQ=255)                           ccAGCTCCCGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGc                                                                                                     <  1:305691/140‑1 (MQ=255)                                tCCCGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACaaa                                                                                                <  1:41567/140‑1 (MQ=255)                                  ccGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGa                                                                                              >  2:431683/1‑140 (MQ=255)                                        cgcgAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGaa                                                                                        <  2:132394/140‑1 (MQ=255)                                                         gcggcATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCCGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACg                                                                       >  2:209302/1‑140 (MQ=255)                                                                   ccAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGa                                                             <  1:337957/140‑1 (MQ=255)                                                                     aGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATg                                                           <  1:196147/140‑1 (MQ=255)                                                                     aGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATg                                                           <  2:330824/140‑1 (MQ=255)                                                                       tACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATGcc                                                         >  2:410767/1‑140 (MQ=255)                                                                            cTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATGCCAATaa                                                    >  2:45362/1‑140 (MQ=255)                                                                               gccgatactgccgatacGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATGCCAATAAaca                                                 >  1:110157/1‑140 (MQ=255)                                                                                        gccgatacGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATGCCAATAAACACTGgcagca                                        <  1:410767/140‑1 (MQ=255)                                                                                                                   aaTCAGCGGTCCGCCCGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATGCCAATAAACACTGGCAGCATCTGAACCTGTTGCAGTTTTTGCTgcg             <  2:79147/140‑1 (MQ=255)                                                                                                                     tCAGCGGTCCGCCAGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATGCCAATAAACACTGGCAGCATCTGAACCTGTTGCAGTTTTTGCTGCGCa           <  1:431683/140‑1 (MQ=255)                                                                                                                              cgccagccAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATGCCAATAAACACTGGCAGCATCTGAACCTGTTGCAGTTTTTGCTGCGCATTCCCCAGc  >  1:125038/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                        AACGACCGAGAGGAACAGCACGATCCCCAGCTCCCGCAGCGCGAGGTTGGCGCTTGGCGGCATAAACCAGTACAGCTTGCCGATACTGCCGATACGCCCGAGGATCAACGCCATAATCAGCGGTCCGCCCGCCAGCCCCAGTTTCAACGCGGCCGGGAATCCTGGCACAAAGACGGGAATAGAACCTAACAATACGCCTAGCCCGATGCCAATAAACACTGGCAGCATCTGAACCTGTTGCAGTTTTTGCTGCGCATTCCCCAGC  >  NC_000913/3864936‑3865200

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