Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,032,787	T→C	V99A (GTA→GCA)	yigZ →	UPF0029 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,032,787	0	T	C	95.3%	143.0 / NA	44	V99A (GTA→GCA)	yigZ	UPF0029 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base C (16/25);  minor base G (2/0);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAACCGAATATACTTTGCAATGTGAATATCATCAGT  >  NC_000913/4032644‑4032933                                                                                                                                                |                                                                                                                                                   ccATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCTGGGAAATTACCGCTGTGGCAGTgcg                                                                                                                                              <  2:146283/150‑1 (MQ=255)      gcgtggcgtgGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGctact                                                                                                                                         <  1:519829/150‑1 (MQ=255)             gtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTAcggcgg                                                                                                                                  >  1:62581/1‑150 (MQ=255)             gtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACggcggc                                                                                                                                 <  1:44301/151‑1 (MQ=255)              tggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGGAATTACCGCTGTGGCAGTGCGCTACTAcggcgg                                                                                                                                  >  1:265203/1‑149 (MQ=255)               gggTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCat                                                                                                                               <  1:4807/151‑1 (MQ=255)                   cgcgGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTg                                                                                                                           <  2:375802/151‑1 (MQ=255)                    gcgGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCCGCCAATGGGCAGCGGCGGCGGGGGAATTACCGCTGTGGCAGGGCGCTACTACGGCGGGATAtt                                                                                                                            >  2:393603/1‑149 (MQ=255)                      gggTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAAACCGTGCCCGCCCAGCTAATGGGCAGCGGCGGCGGGGGAAATACCCCTGTGGGAGCGCGCTACTACCGCGGCATATTGCtt                                                                                                                        >  2:334433/1‑151 (MQ=255)                          gcgcCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCa                                                                                                                    <  2:77874/151‑1 (MQ=255)                             ccGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGGAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTTGCACCg                                                                                                                 >  1:640549/1‑151 (MQ=255)                               ggATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGt                                                                                                               <  1:581206/151‑1 (MQ=255)                               ggATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGt                                                                                                               <  1:464510/151‑1 (MQ=255)                                           aCAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCCCGCCCAGCTAATGGGCAGCGGCGGCGGGGAAATTACCGGCGTGGGAGTGGGCTCCCACGGCGGCATATTGCTTGGCCCCGGCGGGTTAGTGaa                                                                                                    >  1:569688/1‑150 (MQ=255)                                                  ggTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATg                                                                                            <  1:367551/151‑1 (MQ=255)                                                           gacgacGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGgcggcggcg                                                                                   <  1:182699/151‑1 (MQ=255)                                                                 ggggAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGGGGGTTAGTGAAAGCGTATGGCGGCGGCGGGaa                                                                               >  1:543337/1‑149 (MQ=255)                                                                   ggAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAg                                                                           <  1:461210/151‑1 (MQ=255)                                                                      gCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGgcg                                                                        <  2:499949/151‑1 (MQ=255)                                                                       ccGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGgcgc                                                                       >  1:242157/1‑151 (MQ=255)                                                                           cggGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTgcg                                                                   <  2:382810/151‑1 (MQ=255)                                                                              gAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCa                                                                <  2:58616/151‑1 (MQ=255)                                                                                 cGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCt                                                             <  1:557915/151‑1 (MQ=255)                                                                                        tAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACAGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGAc                                                       >  1:608458/1‑150 (MQ=255)                                                                                           aCCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCAGATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCaa                                                   <  2:512929/151‑1 (MQ=255)                                                                                            ccGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAAc                                                  <  2:254741/151‑1 (MQ=255)                                                                                             cGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACg                                                 >  1:31540/1‑151 (MQ=255)                                                                                               aTGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGc                                                >  1:224265/1‑150 (MQ=255)                                                                                                    cGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACg                                          >  2:14551/1‑151 (MQ=255)                                                                                                     gCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTTAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACTACCCAACGCAAGAc                                           <  2:628412/149‑1 (MQ=255)                                                                                                     gCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGc                                         <  2:19978/151‑1 (MQ=255)                                                                                                            tAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAAc                                  >  1:28544/1‑151 (MQ=255)                                                                                                             aaTGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCCCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAAcc                                 <  2:310117/151‑1 (MQ=255)                                                                                                             aaTGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAAcc                                 <  1:184069/151‑1 (MQ=255)                                                                                                               tGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCTGTGAGTTAGTAAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTaa                                   <  2:479155/147‑1 (MQ=255)                                                                                                                    aGCGGCGTCGGGGAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAACCGAatat                           >  2:407978/1‑150 (MQ=255)                                                                                                                             ggggAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAACCGAATATACTTTGc                    <  2:455615/148‑1 (MQ=255)                                                                                                                             ggggAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAACCGAATATACTTTGc                    >  1:455627/1‑148 (MQ=255)                                                                                                                                gAAATTACCGCTGTGGCAGTGCGCTACTACGGCGGCATCTTGCTTGGCACCGGGGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTCACGACCCAACGCAAGACGCCATTAACCGAATATACTTTGCAAt                 >  2:45519/1‑148 (MQ=255)                                                                                                                                  aaTTACCGCTGGGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAACCGAATATACTTTGCAATGTGaa            <  2:7878/151‑1 (MQ=255)                                                                                                                                      aCCGCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAACCGAATATACTTTGCAATGTGAATatc        >  2:431400/1‑151 (MQ=255)                                                                                                                                         gCTGTGGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAACCGAATATACTTTGCAATGTGAATatcatc     <  2:508389/151‑1 (MQ=255)                                                                                                                                           tgtgGCAGTGCGCTACTACGGCGGCATATTGCGTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAAACAGGCGCTGCGACAGCTAACGACCCAACGCAAGACGCCATTAACCGAATATACTTTGCAATGTGAATatcatc     >  2:411706/1‑149 (MQ=255)                                                                                                                                            gtgGCAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAACCGAATATACTTTGCAATGTGAATATCATCAGt  <  1:196446/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   CCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTGCGCCAGCTAACGACCCAACGCAAGACGCCATTAACCGAATATACTTTGCAATGTGAATATCATCAGT  >  NC_000913/4032644‑4032933

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