Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	968,252	G→A	K544K (AAG→AAA)	msbA →	ATP‑binding lipopolysaccharide transport protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	968,252	0	G	A	97.0%	94.1 / ‑5.7	33	K544K (AAG→AAA)	msbA	ATP‑binding lipopolysaccharide transport protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (18/14);  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 = 9.92e-01

ATAAGATGGATAACGGTCTCGATACAGTGATTGGTG‑AAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCG  >  NC_000913/968011‑968493                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   aTAAGATGGATAACGGTCTCGGTACAGTGATTGGTG‑AACACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGGGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCGCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGCTACCGAATCCGATCGTGCGATTCAGGCGGAACTGGATGAGTTGCAGAAAAACCGCACCTCTATGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGa                                                                                                                                                                                                                                            >  2:117793/1‑250 (MQ=255)            aaCGGCCTCGATACAGTTAGTG‑TGTAAAACGGCGTGCTGCTCTCTGGCGGTCTGCGTCAGCGTATGTCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCGGGACGAAGCTACCTCGGTTCTGGATACCGAAGCCGAACGTGCGATTCAGGCGGCTCTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTTCCATTGAAAAAGCATCCGAAATCTTGgtcc                                                                                                                                                                                                                                 <  2:51721/250‑2 (MQ=255)               ggTCTCGATACAGTGATTGGTG‑AAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATGCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAA‑CCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAg                                                                                                                                                                                                                            <  1:15768/251‑1 (MQ=255)                        aCAGTGATTGGTG‑AAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTTCCATTGAAAAAGCAGACGAAATCGTGGTCGGCGAGGATGgtgt                                                                                                                                                                                                                    <  2:196424/251‑1 (MQ=255)                             gATTGGTG‑AAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGAGTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATtg                                                                                                                                                                                                               <  1:264942/251‑1 (MQ=255)                                         aCGGCGAGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGAACGAAGCTACCTCGGATCTGGATACCGAATCCGAACGTGCGAATCAAGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGt                                                                                                                                                                                                    >  1:142941/1‑251 (MQ=255)                                          cGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGt                                                                                                                                                                                                    >  1:81321/1‑250 (MQ=255)                                              gTGCAGCTCTCTCGCGGTAAGCGTCACCCTATTGCTATCGCATGAGACTTGTTGCGTGCTAGACCGCTTCCGGTTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGCATCAGGCCGCAATGGCTGAGTATCAGAAAAACCGTACCTCTCTGGTGATTGCCAACCGCTTGTCTCACATTGAAAAAACAGACGAAATCGTGGTCGTCGAGTATGGGGTCGTTGAGGAACGCggagcgc                                                                                                                                                                                                >  2:188099/1‑245 (MQ=255)                                                                  gcgtcagcgtATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCAcc                                                                                                                                                                           <  1:141963/251‑1 (MQ=255)                                                                     tcagcgtATTCCTATCGCTCGAGCCTTGTTGCGTGAAAGCTCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTCCCATTGTAAAAGCAGCCGAAATCGTGGTCGTCGAGGCTGGTGGCATTGTGGAACGCGGTACGCCTAACGACTTGCTTGAGCACcgcg                                                                                                                                                                        >  2:100505/1‑251 (MQ=255)                                                                       agcgtATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCgcggc                                                                                                                                                                      <  1:196379/251‑1 (MQ=255)                                                                                  aTCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACgcgc                                                                                                                                                            >  1:49842/1‑250 (MQ=255)                                                                                                          aGCCCGATTCTGATTCTGGACGAAGCTTCCTCGGCTCTGGATACCGAATCCGAACGTGCGAGTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCa                                                                                                                                   <  2:63828/251‑1 (MQ=255)                                                                                                              cGATTCTGATTCTCGACGAAGCTACATCGGCTCTAGTTACCGAATCCGAACGTGAGATACAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGATAAAGCCGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGAGGCGTTTACGCGCCACTTCACAAAATGCAGTTTGGCCAATg                                                                                                                                >  2:121838/1‑250 (MQ=255)                                                                                                                           ggACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGAGGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCGTGTCTACCATGGAAAAAGCAGACGAAAGCGTGGTGGTCGAGGATGGTGTCAGGGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGGTTACGCGGAACATCACATTATGCAGTTTgtcgatttagcgataaaat                                                                                                                     <  2:139654/248‑19 (MQ=255)                                                                                                                              cGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAAtctggtc                                                                                                               <  2:81321/251‑1 (MQ=255)                                                                                                                                  gCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGTTGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATctggtctggt                                                                                                           >  1:107164/1‑251 (MQ=255)                                                                                                                                    tACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTg                                                                                                          >  2:67728/1‑250 (MQ=255)                                                                                                                                      ccTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGCTGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGCTGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCACCTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAAt                                                                                                       >  2:222481/1‑251 (MQ=255)                                                                                                                                                ggATACCGAAGCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCGCGGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTtgtg                                                                                             <  1:117793/251‑1 (MQ=255)                                                                                                                                                    aCCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTggcg                                                                                          <  2:85703/250‑1 (MQ=255)                                                                                                                                                       gAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCCCAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCt                                                                                       >  2:3722/1‑250 (MQ=255)                                                                                                                                                               aCGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGCGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGCTGcc                                                                             >  1:227955/1‑251 (MQ=255)                                                                                                                                                                          aGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGCTGCCCCTCTCCTGGt                                                                  >  1:145023/1‑251 (MQ=255)                                                                                                                                                                              ggcACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGCTGCCACTCTCCTGGTTGt                                                               >  1:48023/1‑250 (MQ=255)                                                                                                                                                                                 aCTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGCTGCCACTCTCCTGGGTGTATgg                                                           >  1:70002/1‑251 (MQ=255)                                                                                                                                                                                       tGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTATGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGATGCCCTTCTCCTGGTTGTATGGCCTggg                                                     >  1:228397/1‑250 (MQ=255)                                                                                                                                                                                               aGAAAACCCGTACCGCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGAAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCTCCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGgcg                                             <  1:3722/251‑1 (MQ=255)                                                                                                                                                                                                   aaaCCGTACCTCTCGGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGAt                                         <  1:67728/251‑1 (MQ=255)                                                                                                                                                                                                            ctctctGGTGCTTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAACTCGTGGTCGTCGAGGTTGGTGTCGTTGTGGACCGCGGTACGCATACCGATTTGCTTGAGCACCGCGGCGTTTCCGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGCTGCCACTCGCCTGGTTATATGGCCTGGTGAGTGGCGCGATCCGTCTTTg                                >  1:178497/1‑251 (MQ=255)                                                                                                                                                                                                                             aCCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGCGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGTCGGCT‑ATTGACCCCACTCTCCTGGTTGTATCGTCTTGTGAGTGGGTCGATGCGGCTTTGcgttaaaagaaacact                >  2:88053/1‑235 (MQ=255)                                                                                                                                                                                                                                  ttGTCTATCATTGAGAAAGAAGACGAAATCGCGGCCGTTAAGGTTGGAGTCATTGTGGAAAGAGGTACGCATAACGCTTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCCGTTTGGCCAATGATCGAAAACATCTGGTCTGGAGAATCCCCTTTGTGGCGGCTCATT‑CTGCCACTCTCCTGGTTGTATGGCCTGGGGCGTGGTGGGTTCCGTCTTTGCTATAAACTaaaatttaagggg          >  2:103000/1‑242 (MQ=255)                                                                                                                                                                                                                                          cATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCg  <  2:70002/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   ATAAGATGGATAACGGTCTCGATACAGTGATTGGTG‑AAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCT‑ATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCG  >  NC_000913/968011‑968493

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