Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,585,149	Δ1 bp	16.7%	coding (1338/2280 nt)	ydeP ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,585,149	0	G	.	16.7%	88.9 / 19.9	36	coding (1338/2280 nt)	ydeP	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (19/11);  new base . (4/2);  total (23/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01

CTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGAT  >  NC_000913/1585011‑1585292                                                                                                                                           |                                                                                                                                                 cTTAGTGGCTACGTGTACCGCCAAATCTACTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAACTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cag                                                                                                                                       <  2:604052/149‑1 (MQ=255) cTTAGGGGCTACGTGTACCGCCAAATCTAATTGCGTTCACGGCACAGCTCGCGCTTCCCCGGCTGTCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑cag                                                                                                                                       <  2:44703/149‑1 (MQ=255)    aGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGCTCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGCTCAATGCCCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CCGCAt                                                                                                                                    >  2:240728/1‑149 (MQ=255)    aGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCAt                                                                                                                                    <  2:571392/149‑1 (MQ=255)    aGTGGCTACGTGTACCGCCAAATCTAATGGCGTTAACGGTACAGCGCTCGCTTCCCGCTCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCAt                                                                                                                                    <  2:509681/149‑1 (MQ=255)      tGGCTACGTGTACCGGCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTACCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCCGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCCATTG‑CCGCgtgt                                                                                                                                  >  2:434709/1‑149 (MQ=255)         cTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCa                                                                                                                               >  1:700326/1‑149 (MQ=255)          tACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAg                                                                                                                              <  2:80097/149‑1 (MQ=255)            cGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCACGt                                                                                                                            >  2:747414/1‑149 (MQ=255)                  ccGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGa                                                                                                                      >  2:683290/1‑149 (MQ=255)                              aTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCGGGTGCATGAGGTGGGGTGAAg                                                                                                          >  1:643302/1‑149 (MQ=255)                              aTTGCGTCAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                                                                          <  1:586025/149‑1 (MQ=255)                                          gTACGGCGCTCGCTCCCCGGTCTGGGGTTGCCAGCGCAAAATGGCCCCCCAGGCAGATCAAGGCTCGGGCCTGCCCCGTACAGATGGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCa                                                                                              <  2:509962/149‑1 (MQ=255)                                          gTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCa                                                                                              <  2:230630/149‑1 (MQ=255)                                                                tGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCCGACGAGCCAGAAACTCTg                                                                        >  2:760325/1‑149 (MQ=255)                                                                                   ttGCCCCCCATGCAGATCCATGCCCGAGCCTGCCCCGTACATATCGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACCATGCAGACGGTTTGTCGGGGAt                                                     >  2:162320/1‑149 (MQ=255)                                                                                        ccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGa                                                <  1:708341/149‑1 (MQ=255)                                                                                        ccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGa                                                >  1:461443/1‑149 (MQ=255)                                                                                         cccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCCGCCGGGTTCTCGGTGATACCGAc                                               >  2:439955/1‑149 (MQ=255)                                                                                                 gATCAATGCTCGAGCCTGCCCCGTACACAGTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGGGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGCTTCTCGGTGATACCGACGGTTCGGt                                       >  2:559748/1‑149 (MQ=255)                                                                                                 cagcAATGCCCGCGCCAGCCCCCCACATATTGCTTGCATGCTGGCAATCG‑CAGCCGGCCCAGGTGCATGCGGTGGGGTGAAGCCAGAGCGCCCACCCAGACGAGCCAGGAACTCTCCAGACGGTTTCTCGGTGATACCGACGGTTCGGt                                       <  2:700877/146‑1 (MQ=255)                                                                                                  aTCAATGCTCGAGCCTGCCCCGTACATATTACTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCAACAAGATCTACGGAGGGTTGCTCGGTGATACCGACGGTTCGGTc                                      >  1:409050/1‑149 (MQ=255)                                                                                                  aTCAATGCTCGAGCCTGCCCCGGACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTc                                      >  1:75564/1‑149 (MQ=255)                                                                                                          tCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTa                              >  1:114722/1‑149 (MQ=255)                                                                                                                 tGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                       <  1:616689/149‑1 (MQ=255)                                                                                                                   ccccTTACATATTGATTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAgt                     <  1:21016/149‑1 (MQ=255)                                                                                                                   ccccGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAgt                     >  1:64421/1‑149 (MQ=255)                                                                                                                   ccccGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAgt                     >  1:64346/1‑149 (MQ=255)                                                                                                                                  ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGGCGGTTCGGTTGCCCTGTACATTAGAGTGTCCACGTAGTGGGc      >  1:145955/1‑149 (MQ=255)                                                                                                                                  ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGCCGGTTCGGTGGCGCTGTGAATTAGGGTGTGCACGTAGTgggg      >  2:74873/1‑148 (MQ=255)                                                                                                                                  ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGc      >  1:6306/1‑149 (MQ=255)                                                                                                                                  ggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGc      >  1:560915/4‑149 (MQ=255)                                                                                                                                  ggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGc      >  1:714404/4‑149 (MQ=255)                                                                                                                                     cctGCTGGCAATTG‑CAGCATGTCCAGGTTCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATAAGAGTGTCCACGTAGTGGGCAGa   <  1:434709/147‑1 (MQ=255)                                                                                                                                     cATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGGGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGa   >  2:270746/1‑149 (MQ=255)                                                                                                                                      aTGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGAt  >  2:650232/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                 CTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGAT  >  NC_000913/1585011‑1585292

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