Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,296,481	G→A	V25M (GTG→ATG)	dolP →	division and outer membrane stress‑associated lipid‑binding lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,296,481	0	G	A	97.1%	93.7 / ‑4.6	34	V25M (GTG→ATG)	dolP	division and outer membrane stress‑associated lipid‑binding lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (16/17);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAGTGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCAAACAGATTGCTATGGGCGTAGACGGTGCCAACGA  >  NC_000913/3296234‑3296725                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     gACCGGCAATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATgg                                                                                                                                                                                                                                                   >  2:272400/1‑251 (MQ=255)   ccGGCTATGACGGCGGCGAACGTGCATGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATGCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGt                                                                                                                                                                                                                                                 <  2:104048/251‑1 (MQ=255)      gCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTAcc                                                                                                                                                                                                                                              <  1:39926/251‑1 (MQ=255)                                                         tCCGCATTCCTTCGCATCGAAGTGCTCGCCTTCAGGAAATGCATATGCTGACGGTAAATTGCCTGAGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGATGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGCAGGCAccc                                                                                                                                                                                             >  2:143509/1‑249 (MQ=255)                                                                    tCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCGCCCAGGTGGACGAt                                                                                                                                                                                  <  1:267822/249‑1 (MQ=255)                                                                    tCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATgg                                                                                                                                                                                >  1:205513/1‑251 (MQ=255)                                                                             aGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTgg                                                                                                                                                                        >  1:74846/1‑250 (MQ=255)                                                                               tGCTCGCATTCAGGCAATGCATATGCTGACGGTAAATTGCCTGAGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATAAATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGTTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCCGTGTCGGCACCCAGGGGGACGATGGTACCCTGGa                                                                                                                                                                       >  2:286998/1‑249 (MQ=255)                                                                                       ttCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCAGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTCGACGATGGTACCCTGGAAGTGCGCGTg                                                                                                                                                             >  2:204023/1‑251 (MQ=255)                                                                                           ggAAATGCTTATGCTGACGGTAAATTGCCTGTGCGAGCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACAGGAAGGCAGTATCGCCAATCGCAGTCCTTATTGCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCGGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCGTGGAAGTGCGCGTGAAcc                                                                                                                                                         <  2:29145/251‑2 (MQ=255)                                                                                           ggAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCAGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGCATACATGAAGGCATTATCGCCAATCGCCGTCCTTATTTCCGCGATGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCAGCAACTGACCCACGCAGGGTAGGCACCCAGGTGGACGATGGGACCCTgagcgaggccgtgaac                                                                                                                                                          >  2:248031/1‑235 (MQ=255)                                                                                                    tatGCTGACGGTAAATTGCCTGTGCGATCTGCTCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCACTCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCCCCCAGGTGGACGATGGTACCCTGGACGTGCGCGTGAACAGCGCAtggt                                                                                                                                                >  2:253299/1‑248 (MQ=255)                                                                                                                  aTTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAg                                                                                                                                  <  1:204023/251‑1 (MQ=255)                                                                                                                     gcttgtgCGAACTGATCGATAACACTCTTTTCCCTCACCAGGAGGATTAAGGAGAAGACTTGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGAtt                                                                                                                               <  2:235849/248‑1 (MQ=255)                                                                                                                     gCCTGTGCAATCTGATCGATAACACGCTTTTCCCTCACCAGGATTATTAAGGAGAATACATGAAGGCAGTATCGCCAATCGCAGTCCTTATTGCCGCGCTGCTGGTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGAtt                                                                                                                               <  2:225526/251‑1 (MQ=255)                                                                                                                              atatGATCGATAACACGCTTTTCCCTCACAAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGATGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAGTGGGTACCGCTGCTGTGGGTACCAAAGTCGCAACTGACACACGCAGTGTCGGCACCAAGGTGGAAGAAGGGACCCAGGACGTGCAGGATAACGGCGCATTGTCGAAAGACGaccgggtaaagacagc                                                                                                                       >  2:106015/4‑235 (MQ=255)                                                                                                                                  gATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAgc                                                                                                                    <  1:273915/249‑1 (MQ=255)                                                                                                                                     cGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGTCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTGGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCa                                                                                                               <  2:239448/251‑1 (MQ=255)                                                                                                                                         aaaacGCTTTTCCCTCACAAGGATGATTAAGGAGAATACATGAAGGCCTTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGCCGATGGTACCCTTGAAGAGCGCGTGATAAGCGCATTGGCGAACGAAGAACAGACTAAGCATGAAGAGCGCCTTa                                                                                                            >  2:202007/4‑250 (MQ=255)                                                                                                                                         aaCACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTaa                                                                                                           <  1:225186/251‑1 (MQ=255)                                                                                                                                           cacGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATg                                                                                                         <  1:271137/251‑1 (MQ=255)                                                                                                                                           cacGCTTTTCCCCCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCCCCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATg                                                                                                         <  2:214573/251‑1 (MQ=255)                                                                                                                                                          ccAGGATGAGTAAGGAGAATACATGAAGGCATTATCGCAAATCGCGGTCCTTATGGCCGTGCTGCTGGTGCAAGTTGGTGTTGCCGCTGCCGTAAGGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCTCGCAGTGTCGGCAACCAGGTGGACGATGGTACCCTTGAAGTGCGCGTTAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCTTATCAgg                                                                                          <  2:13341/251‑1 (MQ=255)                                                                                                                                                           cagtatgatTAAGGAGAAGACATGAAGGCATTATCGCCTATCGCAGTCCTTAGTTCCGCGTTGCTGTTGCAAGGTGGTGTTGCCGCTTCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCTCCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAggg                                                                                         <  2:12422/247‑1 (MQ=255)                                                                                                                                                              gatgatTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCa                                                                                       >  1:7317/1‑250 (MQ=255)                                                                                                                                                                              aCATGAAGGCATTATCGCCAATCGCAGTCCTTATTGCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTgg                                                                       <  1:261350/250‑1 (MQ=255)                                                                                                                                                                                     ggCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCAc                                                               <  1:155006/251‑1 (MQ=255)                                                                                                                                                                                      gCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCAcc                                                              >  1:102203/1‑251 (MQ=255)                                                                                                                                                                                          tATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTAGTGTTGCCGGTGCCGTAATGGGTACCGCTGCTGTGGGTACCAACGCCGCAACTGACCCACGCAGGGTCGGCACCCAGGCGGCCGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACCGATTAAGAAAGAAGCGCGCATTAATGTATCGGCATATGCGGGAAAAGTGCTGCTGGTTGGGCAGtacgccaa                                                           >  2:116306/1‑243 (MQ=255)                                                                                                                                                                                                          ccTTAATTCCACGCTGGTTTTGCAAGGTTTTGTTGCCGCGGCCGGAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCGCCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCt                                           <  2:288192/250‑1 (MQ=255)                                                                                                                                                                                                                      gctgctGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCaaaa                                 >  2:288235/1‑247 (MQ=255)                                                                                                                                                                                                                                           tGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGCACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGATGGACGATGGGACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGCAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCAAACTGATTGCTATGGGCGGAGACGGTg         >  1:99022/1‑251 (MQ=255)                                                                                                                                                                                                                                            gCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCCGTGTCGGCACCCAGGTGGACGATGGTACCCTGGCAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTAACCACATGCTGAACTTTCGGCTCGCGCCAAACAGTTTGCTAAGGGCgaagccgggtc        >  2:232248/1‑241 (MQ=255)                                                                                                                                                                                                                                                    cGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCAAACAGATTGCTATGGGCGTCGACGGTGCCAACGa  >  1:147773/1‑249 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     GACCGGCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAGTGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCAAACAGATTGCTATGGGCGTAGACGGTGCCAACGA  >  NC_000913/3296234‑3296725

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