Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,867,551	T→G	M1M (ATG→CTG) †	rpoS ←	RNA polymerase sigma factor RpoS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,867,551	0	T	G	97.0%	104.8 / ‑6.4	33	M1M (ATG→CTG) †	rpoS	RNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (18/14);  minor base C (1/0);  total (19/14)
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.00e+00

AAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATC  >  NC_000913/2867315‑2867796                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      aaaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCTTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAATCCTCAACTCCGTTCTCATCAAATTCCTCATCTTCATTTAAATCAGGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTAc                                                                                                                                                                                                                                         <  1:20623/251‑1 (MQ=255)   aaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGACGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTATCCTCAAATTCCGCCTCTTCATTTAAAACATGACCTTTCAGCGAATTCTGACTCACAAGGTGGCTCCTAccc                                                                                                                                                                                                                                       >  2:339943/1‑251 (MQ=255)     tAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACACATCTCACAAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGAGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTACGTTCTCATCAAATTCCGCATCTTCCTTTAAATCATGAACTTACAGCGTATTCTGACTCAGAAGGTGGCTCCTACACGt                                                                                                                                                                                                                                     >  1:250948/1‑251 (MQ=255)            cttcttcGGCCGTGAACAGTGGTGAATAACCAATCTAACCCAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCTGCCAAATCGTTATCACTGGACTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACATCAACTCCGTTCTCATCAAATTCCGACTCTTCCTTTAGATACTACACTTACAGCGTATTCTGACTCAGAAGGTGGCTCATACCCGTGATCCCt                                                                                                                                                                                                                              >  2:83890/1‑251 (MQ=255)              tcttcGGACGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGGCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATACCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAATGCCTTTTCGTCAAAAACCTCAAATCCGTTCTCATCAAATTCCGCATCTTCATTTAATTCATGAACTTTCCGCGTATTCTGACTCAGACGGTTGATCCTCACCGTGATCCCTTg                                                                                                                                                                                                                            >  2:218616/1‑251 (MQ=255)               cttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTg                                                                                                                                                                                                                            >  1:215379/1‑250 (MQ=255)                      cGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAAc                                                                                                                                                                                                                     >  1:227673/1‑250 (MQ=255)                        ttAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACCCACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTACGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCGTGCACTTTCCGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAtt                                                                                                                                                                                                                  >  1:223575/1‑251 (MQ=255)                                ggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAACCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCACGCaa                                                                                                                                                                                                           >  2:126721/1‑250 (MQ=255)                                     aTAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCGTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCt                                                                                                                                                                                                     <  2:258910/251‑1 (MQ=255)                                                            gCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCg                                                                                                                                                                              <  1:126721/251‑1 (MQ=255)                                                                   cccttCCATCACACGCTGTTTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGAAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATa                                                                                                                                                                       <  2:226584/247‑1 (MQ=255)                                                                     cgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAc                                                                                                                                                                     <  1:339943/251‑1 (MQ=255)                                                                               acGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGtt                                                                                                                                                           >  1:27093/1‑251 (MQ=255)                                                                                      gtggatcacTGCGATAACCGTTACTCTTCGGCCAAATCGTTATCACAGGGATCCTGTTCTACTAAGGCCTTTTCGTTAAAAACATCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAACTCATCACATTTCCGAGTATTCTGACTCAGAAGGATGCTCCTACCCGTGATCCCTTTACGGAACCTTCAAGCAACAACCTGGTTCCGCAGATTTCTCGCTGCGGCAAATAACGTAACGGGATTAAGGAt                                                                                                                                                    >  2:210827/9‑251 (MQ=255)                                                                                        ggttCCCTGCTATAACAGTTCCTCTCCGGCCAAACCGTTATCCCTGGTTTCCTGTACGACAAAGTCCTTTTTGGAAAAATCCTCAACGCCGTTCTCATCAAAAGCCGCATCTTCATTTAAACCATGAGCTTTCAGCGTATTCTGACACAGAAGGTGGCTCCTCCCCGTGAGCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCTCCGATTTATCACTTCGGCATATATCGCAGCGTCTTGACGGAtt                                                                                                                                                   <  2:35931/247‑1 (MQ=255)                                                                                            cccTGCGATACCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCACAACTCCGTTCGCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCTTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCTTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTcccc                                                                                                                                              <  1:211449/251‑1 (MQ=255)                                                                                                   aTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAAc                                                                                                                                       >  1:234561/1‑251 (MQ=255)                                                                                                                        aaTCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgtg                                                                                                                  <  1:54622/251‑1 (MQ=255)                                                                                                                                             gTTCTACTAAGGCCTTTTCGTCAAAATCCTCAACTCCGTTCTCATCAAATCCCGCTTCTTCATTTAAATCATGACCTTTCAGCGTATTGTGACTCAGAAGGTGGCTCCTGACCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCTGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGcg                                                                                               <  2:88113/249‑2 (MQ=255)                                                                                                                                                                    aaaaCCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGcccc                                                                      <  1:218616/251‑1 (MQ=255)                                                                                                                                                                        ccTCAACTCCGTTCTCATCAAATTCCGCACCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCt                                                                  <  1:255388/251‑1 (MQ=255)                                                                                                                                                                                 cGTTCTCATCAAATTCCTCCTCTCTATTTAAATCATGAACTTTCAGCGTATGCTGCCCCAGATGGTGGCTCCTCCCCGGGTTCCTTTGACGGAACATTAAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTCCGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCTGGTTGACCTGGTTCCGGTGCTACCAATGGTCGCTATTTTTTGCCCCGCCTTAActtctt                                                         <  2:227673/251‑1 (MQ=255)                                                                                                                                                                                   ttCTCATCAAATTCCGCATCTTCATTTACATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAGAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCAATGGTCGCTATTTTTTGCCCCGCCTTACCTTCttgt                                                       >  2:76248/1‑251 (MQ=255)                                                                                                                                                                                    tctcATCAAATTCCGCATCTTCATTTCAATCATGTACTTTCACCGTATTCTGCCACAGTATGAGGCTCCTCCCCGTCCTCACTTGACGGACGATTCAAGCAAACGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACACAGCGGGTTACCGCATTTCCCCTTGTAACGAATTTCAAAATCCAAGCTTTTTGAACTGGTTCTGGTCCTACCATTGGTCGCTATTTTTTGCTCAGCCTTAActtcttctc                                                      >  2:34716/1‑248 (MQ=255)                                                                                                                                                                                       catcaAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCCAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTgttgtt                                                   >  1:342144/1‑251 (MQ=255)                                                                                                                                                                                                    accTTCATTTAAACCAGTAACTTTCATCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAATCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTAt                                                                                <  2:264666/207‑1 (MQ=255)                                                                                                                                                                                                    aTCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTAt                                                                                >  1:264666/1‑209 (MQ=255)                                                                                                                                                                                                         cATTTAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTCCCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATtgtg                                  <  1:7668/250‑1 (MQ=255)                                                                                                                                                                                                            ttAAATCATGAACTTTCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                              >  1:269028/1‑251 (MQ=255)                                                                                                                                                                                                                     gAACTTTCAGCGTATTCTGACTCAGACGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGATTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCAATTTTTTGCCCCGCCTTACCTTCTTGTTGTTCCCGGACCAGCCTTGTGTCGTTTTGGGCGt                     >  2:215547/1‑251 (MQ=255)                                                                                                                                                                                                                         tttCAGCGTATTCTGACTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGc                 >  1:237279/1‑251 (MQ=255)                                                                                                                                                                                                                                         cTCAGAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAatcatc  >  1:339267/1‑250 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      AAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCCGGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATC  >  NC_000913/2867315‑2867796

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