Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP025268	4,270,732	T→C	*68* (TAA→TAG)	CXP41_RS22090 ←	protein YjbS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP025268	4,270,732	0	T	C	98.4%	201.0 / NA	63	*68* (TAA→TAG)	CXP41_RS22090	protein YjbS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (34/28);  minor base G (1/0);  total (35/28)
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

ACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGA  >  NZ_CP025268/4270605‑4270857                                                                                                                                |                                                                                                                              aCCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATaaa                                                                                                                   >  2:3980794/1‑140 (MQ=255) aCCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATaaa                                                                                                                   >  2:1479077/1‑140 (MQ=255) aCCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATaaa                                                                                                                   >  1:236648/1‑140 (MQ=255) aCCTTATATTGATTACCAAGATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCGTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATaaa                                                                                                                   >  1:1110397/1‑140 (MQ=255)            aTTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGGAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaaaca                                                                                                        >  2:3016701/1‑138 (MQ=255)            aTTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGaaaaa                                                                                                        >  2:944857/1‑140 (MQ=255)              tACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACt                                                                                                      >  1:45736/1‑140 (MQ=255)                 cAATATCAAATAAAAAGGTCCTGTACTGAGAATAATCACCTTTACCCTCCTCGGCAAAACAGCAAAGCGTGCGAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGCCAAACTTTa                                                                                                   >  2:4566795/1‑140 (MQ=255)                   atatCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCACAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACt                                                                                                 <  1:4565878/140‑1 (MQ=255)                   atatCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACt                                                                                                 <  1:1947354/140‑1 (MQ=255)                   atatCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACt                                                                                                 <  1:1046789/140‑1 (MQ=255)                   atatCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACt                                                                                                 <  1:3299287/140‑1 (MQ=255)                    tatCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACtt                                                                                                >  2:1284709/1‑140 (MQ=255)                      tCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGc                                                                                              >  2:4584954/1‑140 (MQ=255)                          aTAAAAAGGTCCTGTACTGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCCTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATa                                                                                          <  1:2835804/140‑1 (MQ=255)                          aTAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATa                                                                                          <  2:2408078/140‑1 (MQ=255)                             aaaaGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGaa                                                                                       <  1:1269896/140‑1 (MQ=255)                               aaGGTCCTGTAATGAGAATAATCACCCTTAACCTCCTCGGCAAAACAGCAAGGCGTGCAAATATACATAAAGGCTTTCTTACCCCTGATTGTGGTAAGTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACt                                                                                     >  1:3052528/1‑140 (MQ=255)                                  gTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGa                                                                                  <  1:2272376/140‑1 (MQ=255)                                  gTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGa                                                                                  <  1:3980794/140‑1 (MQ=255)                                     cTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATAc                                                                               >  2:2307709/1‑140 (MQ=255)                                     cTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATAc                                                                               >  1:1845260/1‑140 (MQ=255)                                            gagaATAATCACCTTTAACCTCCTCGGGAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTtcttct                                                                        >  2:3738745/1‑140 (MQ=255)                                             agaATAATCACCTTTAACCTCCTCGGCAAAACATCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACt                                                                                                 <  1:4065409/114‑1 (MQ=255)                                             agaATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACt                                                                                                 >  2:4065409/1‑114 (MQ=255)                                             agaATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTtctc                                                                       >  2:3272579/1‑140 (MQ=255)                                             agaATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTtctc                                                                       >  2:3325769/1‑140 (MQ=255)                                               aataatCACCTTTAACCTCCTCGTCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCa                                                                     >  2:4717701/1‑140 (MQ=255)                                               aataatCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCa                                                                     <  2:713580/140‑1 (MQ=255)                                                 taatCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGCACTTGATACTTCTTCTCCAgg                                                                   <  2:3959300/140‑1 (MQ=255)                                                  aatCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGt                                                                  >  1:3867553/1‑140 (MQ=255)                                                  aatCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGt                                                                  >  1:1202775/1‑140 (MQ=255)                                                    tCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCt                                                                >  2:4394525/1‑140 (MQ=255)                                                         tttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTaaaaa                                                           >  2:2962643/1‑140 (MQ=255)                                                         tttAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTaaaaa                                                           >  2:4595305/1‑140 (MQ=255)                                                                  ctcGGCAAAACAGCAAAGCGTGCAAATATACATAACGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGc                                                  >  1:3332012/1‑140 (MQ=255)                                                                  ctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGc                                                  >  2:3366920/1‑140 (MQ=255)                                                                  ctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGc                                                  >  1:541773/1‑140 (MQ=255)                                                                  ctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGc                                                  >  1:3373653/1‑140 (MQ=255)                                                                  ctcGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGc                                                  >  1:3330807/1‑140 (MQ=255)                                                                     ggCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGAt                                                                                 <  2:1384741/106‑1 (MQ=255)                                                                     ggCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGAt                                                                                 >  1:1384741/1‑106 (MQ=255)                                                                        aaaaCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAAtttt                                            <  1:1412775/140‑1 (MQ=255)                                                                           aCAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTaa                                         <  1:3094548/140‑1 (MQ=255)                                                                             aGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATCAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                       <  2:3237158/140‑1 (MQ=255)                                                                             aGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                       <  1:4394525/140‑1 (MQ=255)                                                                             aGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                       <  2:3161753/140‑1 (MQ=255)                                                                             aGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                       <  1:3854733/140‑1 (MQ=255)                                                                             aGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                       <  1:316187/140‑1 (MQ=255)                                                                             aGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                       <  1:2876198/140‑1 (MQ=255)                                                                                   gCGTGCAAATATACATAACGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACATGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTCGCACTTTTTAACTGTGttt                                 <  2:3289855/140‑1 (MQ=255)                                                                                      tGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTa                              >  1:2488354/1‑140 (MQ=255)                                                                                           atatACATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGtat                         <  1:3845356/140‑1 (MQ=255)                                                                                              taCATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATaaa                      <  1:456009/140‑1 (MQ=255)                                                                                                cATAAAGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATaaaaa                    >  1:189771/1‑140 (MQ=255)                                                                                                     aGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                       <  2:950312/116‑1 (MQ=255)                                                                                                     aGTCTTTCTTACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACt                                       >  1:950312/1‑116 (MQ=255)                                                                                                        ctttcttACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCtt            <  2:236648/140‑1 (MQ=255)                                                                                                        ctttcttACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAAGTTTTAACTGTGTTTCTATGTATAAAAAGAAATCtt            <  2:827874/140‑1 (MQ=255)                                                                                                           tcttACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTgg         >  1:565982/1‑140 (MQ=255)                                                                                                           tcttACCCCTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTgg         >  1:5022780/1‑140 (MQ=255)                                                                                                               aCCCCTGATTTTGTTAACTAAATCAATAAATACGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATa     >  2:3887871/1‑140 (MQ=255)                                                                                                                  ccTGATTTTGTTAACTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAaga  <  1:4717701/140‑1 (MQ=255)                                                                                                                                |                                                                                                                              ACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATCACCTTTAACCTCCTCGGCAAAACAGCAAAGCGTGCAAATATACATAAAGTCTTTCTTACCCCTGATTTTGTTAATTAAATCAATAAATGCGTGAAAAACTTTACTTGCAATAGAACTTGATACTTCTTCTCCAGGTCTAAAAAATCTTTTGCAATTTTTAACTGTGTTTCTATGTATAAAAAGAAATCTTTGGCATAAGA  >  NZ_CP025268/4270605‑4270857

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