Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,510	G→C	E445Q (GAG→CAA)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,510	0	G	C	91.2%	85.4 / 3.9	34	E445Q (GAG→CAA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base C (19/12);  total (21/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 = 2.33e-01

GTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTT  >  NC_000913/4639271‑4639748                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               gTATTTATTTGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGGGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCCTGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTAAGTATc                                                                                                                                                                                                                                     >  1:390815/1‑251 (MQ=255)      tatttGCAGGCGGTATTGAAAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGa                                                                                                                                                                                                                                 >  1:234528/1‑250 (MQ=255)        tttGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAAt                                                                                                                                                                                                                              >  1:246786/1‑251 (MQ=255)           gcgtgtggTATTGAAAGGCGGGCGCAACAGCATGTTGTGATCCCTCGGGCTGTTGGTGCTGGAGGGTGTTATGTTGGGCCTGCTCAACTCGGCCGATAGCGAGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATAAAGACTGGTATGCGTTTTCACTGCCGAATATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATaa                                                                                                                                                                                                                            <  2:163111/244‑1 (MQ=255)                   tATTGACAGGCTGGCGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGGTCAATTACGTATCTGGAAATAAGGTTGaaa                                                                                                                                                                                                                    >  2:114320/1‑250 (MQ=255)                              tggcgCAACAGAATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGCCTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTATCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGTACAATTACGTATATGGAAATAAGGTTGAAAAATAAAAACggg                                                                                                                                                                                                        >  2:50287/5‑250 (MQ=255)                                  gCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGCAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATaaaa                                                                                                                                                                                                             >  2:426654/1‑242 (MQ=255)                                  gCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATaaaa                                                                                                                                                                                                             <  1:426654/242‑1 (MQ=255)                                       ggCATGTTGTTTACCCACGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCGGCCGATAGCGCGCGGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGAATCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAg                                                                                                                                                                                               <  2:40980/250‑1 (MQ=255)                                                tttACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGtttt                                                                                                                                                                                      >  1:159240/1‑251 (MQ=255)                                                    cccTCGCGCTGTGGTGGATGGATGGTGTGATGTGGGGACTGCTCAACCCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGtttttttt                                                                                                                                                                                  <  1:101359/251‑1 (MQ=255)                                                     cctcggtgagtttgtgttGGATGGGTTGAGGTGGTGACTGTGCAACTCTGCCGATTGCTGGCTGTTGTTGGGAACCAGGGTG‑TGTGTGGTGGTGTTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAAGTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGgttttttga                                                                                                                                                                                 <  2:441729/234‑3 (MQ=255)                                                             tgttgtTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTaaaaa                                                                                                                                                                                                >  1:353107/1‑228 (MQ=255)                                                             tgttgtTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTaaaaa                                                                                                                                                                                                <  2:353107/228‑1 (MQ=255)                                                                    gCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCg                                                                                                                                                                  >  1:158824/1‑251 (MQ=255)                                                                               tgATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTcc                                                                                                                                                       >  2:97757/1‑251 (MQ=255)                                                                                       ggACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCAcc                                                                                                                                               >  2:245464/1‑251 (MQ=255)                                                                                                          aTAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTc                                                                                                                             >  1:289328/1‑250 (MQ=255)                                                                                                          aTAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGACACCCTTc                                                                                                                             >  2:370101/1‑250 (MQ=255)                                                                                                                          gggAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTTGCTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAATTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCga                                                                                                            <  1:286148/251‑1 (MQ=255)                                                                                                                                cAGTGTGGTG‑GGGGGGGCGCTGGCCGGCATGATGTTTGTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGAttt                                                                                                      <  2:208651/251‑1 (MQ=255)                                                                                                                                                             atggtgTTTGTGCCTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGTTCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAACCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGTTGGTGGCGATGATTTCCGGCGTATCCGGGGTAGATTCGAAATg                                                                          >  2:290076/5‑251 (MQ=255)                                                                                                                                                             atgatgTTTGGGAAACATAATAACGAATGGGATGCGTTTTCACTGCCGAAAAGGAAAGCCAGTAAAGAAAGTAAAATGGACGATCAATTACGTAACTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATg                                                                          <  2:480687/251‑1 (MQ=255)                                                                                                                                                                atgTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTGCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGttt                                                                       <  2:246786/251‑1 (MQ=255)                                                                                                                                                                                  atCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGt                                                     >  1:374085/1‑251 (MQ=255)                                                                                                                                                                                    cGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGGGGAAAAGCCGATTATTCTTCCAGTTCACCCCAGAGCCGATAACCTTCCCCGTGAATGGGGGCGTTGATTTCCGCCGTACCCGCCGTAGATTCGAATGGTTTAGGAATCCGCCGGATCgccc                                                   >  1:396217/1‑248 (MQ=255)                                                                                                                                                                                           tATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTAc                                            >  1:6052/1‑251 (MQ=255)                                                                                                                                                                                              gCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAg                                          <  2:158824/250‑1 (MQ=255)                                                                                                                                                                                                    tCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTc                                   >  1:115596/1‑251 (MQ=255)                                                                                                                                                                                                        tGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGCCGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAAGAATACGGCGGCGCGTCACGGCTACAGTAAGGTCgggg                               >  2:418837/1‑248 (MQ=255)                                                                                                                                                                                                            gAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACCGCTACAGTACGGGCGTGTGGtt                           >  2:484733/1‑251 (MQ=255)                                                                                                                                                                                                                 tGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGTTAAAAAACGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATTTTCCAGCTCACCACACAAGCGTTAACCTTCACCGTGGATGGTGGCGCTGATTTCCGGCGTATCCGGAGTAGTTTCGAATTGTTTACGAATACGGCGGATCGTCTCGTCCAAAGTATGGACGTGTGGTTTCAGc                      >  1:499163/1‑251 (MQ=255)                                                                                                                                                                                                                                aagTTACAATGGACGATCAATGACGTATCTGGAAATAAGGGTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTGTTGTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGCAGCGAGATCCTTCACCGTGGATGGTGGCGATGAGTGCCGGCGTATTCGGCGGAGATTCGAAATGTTTACGAATACGGCGGATCGACACGTCTTCAGTACGGTCGTGCGGTTTCGGCGTACTGCCGGTCAtt       <  2:57454/251‑1 (MQ=255)                                                                                                                                                                                                                                     aCAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTttctt  <  1:185456/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               GTATTTATTTGCAAGCGGTATTGAAAGGTTGGTGCAACAGCATGTTGTTTACCCTCGCGCTGTTGTTGCTGGATGGTGTGATGTGGGGACTGCTCAACTCTGCCGATAGCGCGCTGTTGTTGGGAACCAGTGTGCTG‑GTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTT  >  NC_000913/4639271‑4639748

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