Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,520,084	G→A	32.9%	intergenic (+20/‑178)	curA → / → mcbR	curcumin/dihydrocurcumin reductase, NADPH‑dependent/colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,520,084	0	G	A	32.9%	26.2 / 4.2	36	intergenic (+20/‑178)	curA/mcbR	curcumin/dihydrocurcumin reductase, NADPH‑dependent/colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled
Reads supporting (aligned to +/- strand):  ref base G (8/14);  new base A (7/3);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGAC  >  NC_000913/1519863‑1520299                                                                                                                                                                                                                                    |                                                                                                                                                                                                                         gCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGGTGATTAAA‑GTTTTAA‑AGGGGGC‑GTAAGACGCCa                                                                                                                                                                                                                    >  1:347670/1‑228 (MQ=255)         ttattaTCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGGGATGGTTAAA‑GTTTCAA‑CGCGGGC‑GGAAACCGCCATaaaaaac                                                                                                                                                                                                            >  1:656761/1‑226 (MQ=255)                 gCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑ATTGGGCGAAAGAGGATAAAATCCACTACCGCGAAGAAATTAATGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑AGGCGGG‑GTAAACCGGCATaaaaagaccgggaa                                                                                                                                                                                                     >  1:304876/1‑218 (MQ=255)                 gCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCACAGTGGAGATCAGCGGGGCGGGGGATGCGTAAA‑AATTTAA‑AGGCGGG‑GGAAGCGGCCCCaaaaaaaaagagaaa                                                                                                                                                                                                    >  2:368648/1‑218 (MQ=255)                            tGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGTTTAAA‑GTTTCAA‑CGGCGGC‑GTAACGCGCCATAAAAATCCAGGCGAAAAGCTtttttt                                                                                                                                                                                         >  1:12501/1‑228 (MQ=255)                            tGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGGGATGATTAAA‑GTTTAAA‑AGGGGGG‑GTAAACCGCCCTAAAAATCCCAGCAAAAAACTTAtttt                                                                                                                                                                                         >  1:472045/1‑228 (MQ=255)                                 aCCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGGGATGTTTAAA‑GTTTCAA‑CGGCAGC‑GGAAGCCGCCATAAAAAGCCTAGTCACAAGGTTAtttttacct                                                                                                                                                                                    >  1:229592/1‑223 (MQ=255)                                 aCCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGGGATGGTTAAA‑GTTTCCA‑CGGCGGCGGTAAACCGCCCAAAAAATCAGGGCAAAAACTTTTTTTTCAgt                                                                                                                                                                                     >  2:411308/1‑227 (MQ=255)                                 aCCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGGGATCGGCGTGGCGGGTGAGGGTTAAA‑GTTTAAA‑AGGCGGG‑GTAAGGCGCCCTAAAAATTCAGGCAAAAAACTTATTTTTAGCt                                                                                                                                                                                    >  1:252361/1‑228 (MQ=255)                                      aTCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGGGGTGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAACCGCCCTAAAAATCCACGCGAAAAGCTTATTTTCAGCTttaat                                                                                                                                                                               >  1:351581/1‑228 (MQ=255)                                            gggattCAGGGGGAAAAGGGGCAAAGGGGTAAAAGAGGAAAAAACCCCCTCCCGGGAAAAAATTCCGTACGGTTTA‑GAGATGGCCC‑CACAGCCGTTTATCGGCCTGCTGAAGGGTAAAAACTCCGGCAAAGTGGTGACCCGCGTGGCGGGTGATGATTAAA‑GTTTCAC‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaact                                                                                                                                                                          <  1:344224/224‑1 (MQ=255)                                               tttCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGGTTTCAA‑‑GGCGGC‑GGAAACCCGCACAAAAATCCACGACAAAAGTCTACATTCAACTTTCAATAAATAAcc                                                                                                                                                                      >  1:67376/1‑228 (MQ=255)                                                ttgagggggAGAGGGGGC‑AGGGGGTGAAAGAGGAAAAAACCCCACGGCGAGAAGAAATGGCGGGCGGTTGA‑GATAGGGACC‑AACAGGCGTTTACGGGCCTGGGGAAGTAAAAAAATTCCGGCAAATGGGTGCCCCGCGTGGGGGGTAAGGATTAAA‑GTTTCAA‑CGGCGGC‑TTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACca                                                                                                                                                                     <  2:521049/222‑1 (MQ=255)                                                                       ggTGAAAGGGGATAAATGCCACGACCGCAAAGAAATTCGTGACGGTTAA‑GGAAAGGCGC‑AACAGGCTTTTAAGGGCTGGTAGATGGATAAAAACTGAGGAAAAGTGTTAATCGGCGTGGGGGTGGATGAGTAAA‑GTTGCAA‑CGGCGGC‑TTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCATCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                                                                                                                                                <  1:320717/227‑1 (MQ=255)                                                                                     aaaCCCACTCCCGCGAAAAATTTTCGGAGGGTTTA‑AAGAAGGCGC‑CACAGGCGTTTATCGGCCTGGTGAAGGGTAAAACCTTCGGCAAAGGGGTGACCCGCGTGGCGGTTTATGAAAAAA‑GTTTAAA‑CGGCGGC‑TTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                                 <  2:550497/228‑1 (MQ=255)                                                                                      acccccccaCCGGGAAAAAATTTCTGACGGTTAA‑GAGAATCCGC‑CACGGACTTTTACGGGCCTGCGAAAGGATAAAAACTTGGGCAAGGTGGTCACCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑ATAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                                 <  1:57660/219‑1 (MQ=255)                                                                                      aaTCCGCTACCGAAAAAAAATTATTG‑TGGTTTAGGAGATTCCAC‑GAGAGACGTTTACCGGCGTGCTGAAGTGTAAAATCTACGGCAATGTGGTGACCCGCGTGGCGGGTGATGATAAAA‑TTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg                                                                                                                                 <  2:387509/227‑1 (MQ=255)                                                                                          cACTACCGAGAAAAAATTACTGACGGTTTA‑GGAAAGGCGC‑CACAGACGTTTATGGGCCTGGTAAAGGGTAAAAGCTTCGGAAAGGTGGTGATCCGGGTGGCGGGTGAGGATAAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGtt                                                                                                                            <  1:36057/228‑1 (MQ=255)                                                                                                  cgAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGGGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCCTAAAAATCCAGGGAAAAAAGTTATTTTTAGCTTTTATTAACTAACCACATCATTGTTTTGCATCTAATAAAACCAAAAGtatggggg                                                                                                                             >  2:36057/1‑212 (MQ=255)                                                                                                    acgaatatACTGACGTTATA‑GAGAAGTGGA‑CCCAGTAGTTGAGCGACCGGCTGAAGGTTAAAATCTGCGGGAAAGTGGTACTCCGCTGGGGGGGCGATGATTAAA‑TTTTAAA‑GAGCGGC‑ATAGACCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGattatt                                                                                                                  <  1:672892/221‑1 (MQ=255)                                                                                                    aagaaTTTCCTGATGTTTTA‑AGGAAT‑CGCACACGGGCGTTTACCGGCCTGCTAAAGGGAAAAAACTCAGAAAAAGTGGTATACCGCGTGGCGGGTGAAGATTAAA‑GTTTAAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAATTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGattatt                                                                                                                  <  1:368648/228‑1 (MQ=255)                                                                                                            aCTGACGTTTAA‑GAGAATGCGC‑CACGGACGTTTAGCGGCAGGTTGAAGTGTAAAAACTTCGAAAAGGTGGTAATCCGCGTGGCGGGTAATGATTAAA‑GTTTCAC‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTaa                                                                                                          <  1:158553/228‑1 (MQ=255)                                                                                                                aCGGTGTA‑GAGAAGGCGC‑CACAGACGTTTATGGGCGTGTTAAAGGGTAAAATCTTCAGAAAAGTGGTGATCCGGGTGGCGGGTGATGATTAAA‑TTTTCAA‑CGGCGGC‑TAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATaca                                                                                                      <  1:266526/228‑1 (MQ=255)                                                                                                                         gaaaaTGCGC‑CACAGGCTTTAAGCGGCGTGGTAAAGGATAAAAACTGCGGCAAGGTGGTGACCCGGGGGGCGGGGGAAGA‑TAAA‑ATTTCAAGCGGCGGC‑GAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGa                                                                                              <  1:564405/225‑1 (MQ=255)                                                                                                                                                                     actgcgcaaaaGTGGTGATCCGAGGGGCGGGTATTGTTTAAA‑TTCTCAA‑CGGCGGG‑TAAATCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGa                                                  <  1:301026/220‑1 (MQ=255)                                                                                                                                                                      ttgcGGAAAAGGGGTGACCCGCGGGGCGGTGAATGATTAAG‑TTTTCAA‑GGGCGGC‑GAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGa                                                  <  1:352664/224‑1 (MQ=255)                                                                                                                                                                         cGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTTTGCTGCTAATAACAACAAAGGGGAGTGGTTTATGATTATTTACTTTAATAACCCAGGTGAATTCCTTTTGCCATGCAGGCAGGGGTGGAC‑AGAAAACACCCGAAg                                                >  1:353137/1‑226 (MQ=255)                                                                                                                                                                          ggCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCCTAAAAATCCAGGGAAAAAAGTTATTTTTAGCTTTTATTTAACAACCACATCCTTTGTTTTCTGCCAATAAAAACCAAAGTGAGTGGGTTATGATTTTTTACTTTTATATACCAGGGGCATTCTTTTTGCCATGTATGTAGGGTTGGGC‑aaaaaaacaccgcaa                                                >  2:411924/1‑216 (MQ=255)                                                                                                                                                                                 tggtgATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACAAACCACATCATTGTTCTGCTGCTAATAAAAACAAAAGTGAGTGGGTTATTATTATTCACTTTAATACACCAGGGGAATTCCTTTTGCCATGCAGGCAGGGGTGGACAATAAAACAC‑CGAAGAAATTtaaaa                                      >  1:258757/1‑227 (MQ=255)                                                                                                                                                                                              ggggggTAATAATTAAA‑GTTCAAG‑CGGGGGT‑GAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACggg                          <  1:411308/224‑1 (MQ=255)                                                                                                                                                                                              ggCGGGTTATGATTAAA‑GTTACAA‑CGGGGGG‑GAAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACggg                          <  1:411924/227‑1 (MQ=255)                                                                                                                                                                                                      atttatAAA‑TTTTCAA‑CGGCGGC‑ATAAACCGCCATAAAAATCCAGGCAAAAAGCTTATTTACAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGaa                 <  2:353137/223‑1 (MQ=255)                                                                                                                                                                                                           tAAA‑TTTCCAA‑AGGCGGC‑GTAGTCCGGCATAAATATCCAGGCAAAAAGCTTATTTTCATCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACAt             <  1:87492/227‑1 (MQ=255)                                                                                                                                                                                                                    cAA‑CGGCGGC‑GTAAACCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCGTTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTTTTCACTTTATAACACCAGGGGAATTCCCTTCGCCATtcag                                                                            >  2:38884/1‑157 (MQ=255)                                                                                                                                                                                                                       ctgcgggc‑tTAAACCGCCATAAAAATCCAGGAAAAAAGTTTATTTTCAGCTTTAAGTAACTAACCACATCATTGTTCTTCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGAc  <  2:252361/223‑1 (MQ=255)                                                                                                                                                                                                                        acctggtaatAACCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAAATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGAc  <  1:282201/219‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                         GCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGC‑AATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTA‑GAGAATGCGC‑CACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAA‑GTTTCAA‑CGGCGGC‑GTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATGATTATTCACTTTAATACACCAGGTGAATTCCTTCTGCCATGCAGGCAGGGTTGGAC‑AGAAAACACGCGAAGGAATTAAAATGCCTGGAACGGGAAAAATGAAACATGTCAGTTTGAC  >  NC_000913/1519863‑1520299

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