Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,346,840	(T)8→7	coding (45/105 nt)	ymiB →	putative protein YmiB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,346,833	0	T	.	96.8%	118.6 / ‑1.7	31	coding (38/105 nt)	ymiB	putative protein YmiB
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (13/17);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATG  >  NC_000913/1346695‑1346955                                                                                                                                           |                                                                                                                           cGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtt                                                                                                                           <  2:556738/140‑1 (MQ=255) cGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCACTGtt                                                                                                                           >  1:476541/1‑140 (MQ=255)    tCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttt                                                                                                                        <  2:491232/140‑1 (MQ=255)      cGAGAAAGGTGTACAATTTTGTCGACTCCCCCAGGGTTTCCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttttt                                                                                                                      <  2:326983/140‑1 (MQ=255)        agaAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTATCAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttat                                                                                                                    <  2:387167/140‑3 (MQ=255)                        ttGTGGACTCACTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCGAAAAGTATTGCTTTTTTT‑TATTTTTATTTGGaagc                                                                                                   >  2:159418/1‑138 (MQ=255)                                  ccTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGa                                                                                         >  1:568469/1‑140 (MQ=255)                                            ttCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTT‑TATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAg                                                                               >  2:462511/1‑140 (MQ=255)                                            ttCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAg                                                                               >  2:460633/1‑140 (MQ=255)                                            ttCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAg                                                                               >  1:517757/1‑140 (MQ=255)                                                  atgCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCa                                                                         >  1:464315/1‑140 (MQ=255)                                                       tatCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTGTATATTTAATTGGAATACAGTTTTGGAAAATCGCCAGTTTTCAAACGa                                                                    >  2:197056/1‑140 (MQ=255)                                                           ctccGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGt                                                                <  2:568469/140‑1 (MQ=255)                                                           ctccGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGt                                                                <  1:352888/140‑1 (MQ=255)                                                             ccGGTTAACAGCGGAGTCCGCGTTTAAATAAAAGGCACACCATGACTTGCATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCa                                                              <  2:223917/140‑1 (MQ=255)                                                             ccGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCa                                                              >  1:381193/1‑140 (MQ=255)                                                               ggTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCCGTTTTTGAAAAACGCCAGTTTTCAAACGAAAGTCAGt                                                            >  2:464324/1‑140 (MQ=255)                                                               ggTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGt                                                            <  2:476541/140‑1 (MQ=255)                                                                       gCGGAGTCCGCGGTTAAATAAACGGAACAACGTTAGTTAAATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAAt                                                    <  2:289166/140‑1 (MQ=255)                                                                         ggAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCt                                                  >  1:48556/1‑140 (MQ=255)                                                                           aGTCCGCGGTTAAATAAAAGGATCATCATGACTTACAGGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGc                                                <  1:261467/140‑1 (MQ=255)                                                                           aGTCCGCGGTTAAATAAAAGGAACACCATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGc                                                <  1:268176/140‑1 (MQ=255)                                                                           aGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGc                                                <  2:464315/140‑1 (MQ=255)                                                                             tCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTCAATACCTTG‑TTTTTTTCTATATTTTTTTTTAAACCAGCTTTGGGAAAAAGCCCGTTTTCAAACGGAAGTCAGTTAAAAAATCTGGCt                                              >  2:128191/1‑140 (MQ=255)                                                                             tCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCATTTTTTGAAAAACGCCAGGTTGCAACAGTAAGTCAGTTAAAAAATCTGCCt                                              >  1:289166/1‑140 (MQ=255)                                                                             tCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCt                                              <  2:485452/140‑1 (MQ=255)                                                                                   gTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGatat                                        <  1:128191/140‑1 (MQ=255)                                                                                       aataaaAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCGAACGAAAATCAGTTAAAAAATGTTCCTGGATATAACg                                    >  1:171584/1‑140 (MQ=255)                                                                                            aaGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGt                               >  2:446496/1‑140 (MQ=255)                                                                                                  caacaTGACTTACATGATATTATCGGCGGCTCAACCCATCG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCg                         <  2:206794/140‑1 (MQ=255)                                                                                                        gACTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAATAAATCTGGCTGGATGTAACGAGGGTAGAGCGGGGaaa                   >  1:198354/1‑140 (MQ=255)                                                                                                          cTTACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATa                 <  2:381193/140‑1 (MQ=255)                                                                                                            tACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATaaa               <  1:642514/140‑1 (MQ=255)                                                                                                            tACATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATaaa               <  1:159418/140‑1 (MQ=255)                                                                                                              cATGAAATTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACg             <  1:460633/140‑1 (MQ=255)                                                                                                                   aaTTAACGGCGGCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCcca        <  1:446496/140‑1 (MQ=255)                                                                                                                        aCGGCGGCTAAACGCATTG‑TTTCTTTATACATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCccatccat   <  1:197056/140‑1 (MQ=255)                                                                                                                         cggcggCTAAACGCATTG‑TTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATg  >  1:77511/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                           CGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATG  >  NC_000913/1346695‑1346955

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