Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,490,095	G→A	intergenic (+424/‑807)	aldA → / → cybB	aldehyde dehydrogenase A/superoxide oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,490,095	0	G	A	96.0%	69.6 / NA	26	intergenic (+424/‑807)	aldA/cybB	aldehyde dehydrogenase A/superoxide oxidase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (15/9);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTA‑GTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTTATAAACGATAGTTTTCATTTC  >  NC_000913/1489854‑1490341                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        gaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTCTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGat                                                                                                                                                                                                                                                >  1:324314/1‑251 (MQ=255)   aaTGGCTGCAAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATAACTCGTTATTGGAGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATCCCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGGGTTGCGCATGTCCTTTCAGTTTGCCGCGCCGTTCCGGGATCACCATACCACTGGCTTTTGCCTACCCCGTCGTGATGGGAATGatat                                                                                                                                                                                                                                              >  2:189057/1‑251 (MQ=255)          gCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGc                                                                                                                                                                                                                                       <  1:51650/251‑1 (MQ=255)                      ggAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGa                                                                                                                                                                                                                           <  1:207559/251‑1 (MQ=255)                                   tttgctacGGTATAATCAAATGTCTCGGTATTGGTGGTTGCGTGTTTAAGTGCGTTATGCTCCTCTGCTGCAGTCACTTTTTTCCGAGATTAGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGAACAATGGCTTTTGCCGTCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCCCGTAAATc                                                                                                                                                                                                               <  2:117944/243‑1 (MQ=255)                                      tcATCGGTATAACCAAATGACGCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCGCTTCGGCAGTCACTTTTTGCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTGTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCt                                                                                                                                                                                                              <  2:165726/248‑1 (MQ=255)                                                                                              cACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCAt                                                                                                                                                   >  1:61517/1‑251 (MQ=255)                                                                                                  tctTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTg                                                                                                                                               >  1:291700/1‑251 (MQ=255)                                                                                                     tCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTg                                                                                                                                               <  1:256589/248‑1 (MQ=255)                                                                                                      cGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGATCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTACGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCAACCAGTGACTGGGTGCCAGTATAGGCATGAAAGGTCGTCATCGTGCCg                                                                                                                                            >  2:220645/1‑250 (MQ=255)                                                                                                      cGGCAGACACTTTTTTCGGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGGTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAACACGTAAATCTTTAACACGCGGGCCTTCCACCACTGAATGGGTGCCAGTATAGGCCTGTAGGGTCGTCAACGTGcat                                                                                                                                            >  2:269766/1‑248 (MQ=255)                                                                                                       ggCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGAAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGCATGGTCGTCATCGTGCCGAc                                                                                                                                          >  1:178083/1‑251 (MQ=255)                                                                                                                   tttCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGa                                                                                                                                                                                    >  1:176935/1‑197 (MQ=255)                                                                                                                   tttCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGa                                                                                                                                                                                    <  2:176935/197‑1 (MQ=255)                                                                                                                                          aGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGcc                                                                                                        >  1:33159/1‑250 (MQ=255)                                                                                                                                            tGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAACGATATTTTCTGCCGCTGCGCGTGAAGCACGTCAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCTCGCCGCCTTCAATACCGAAACTGTAATCCAAGGCTTTGACCaat                                                                                                     >  2:367607/1‑249 (MQ=255)                                                                                                                                                      tGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGgcgaacgac                                                                                           >  2:323205/1‑245 (MQ=255)                                                                                                                                                           tCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGCCTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGACACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGtt                                                                                       >  2:162561/1‑250 (MQ=255)                                                                                                                                                                            gcatgagatctCAGTTTTCCGATCAGTTCCGGCATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGAGTAGGGAATGATATTTTCTGACGCTGCGAGTGACGCACGTAAATCTTTACCACGCGGTCCATCCACCCGTGCCTGGGGGCCTGTATAGGGATGAATGGTCGTCATCGTGCCGACTTCTGTCCCGCAACTGTCAGACAAGGCTTTCGCCATTGGCGCAAGAC‑GTTAAGTGGTGCAtaaaccct                                                                     >  2:219261/11‑243 (MQ=255)                                                                                                                                                                                                        ccGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTA‑GTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCag                                         <  1:216447/250‑1 (MQ=255)                                                                                                                                                                                                         cggcaTCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCCGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTA‑GTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCagag                                       <  1:161180/247‑1 (MQ=255)                                                                                                                                                                                                                        aaTGGCTTTTGCCGCCCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTA‑GTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGt                        >  1:213157/1‑251 (MQ=255)                                                                                                                                                                                                                                     gccCCCGTAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTA‑GTGGTGCATGACGCCACGGAAACAATGGGGTCGTTGCCATCCAGAGTGTCGTCATTGACGTGATAAACGATAGt           >  1:137716/1‑251 (MQ=255)                                                                                                                                                                                                                                        cacGTAGTGTGGGGAATGATATGTTCTGCCGCTGCTCGTGAAGCACGTAAATCTGTACCACGCGGGCCATCCACCAGGCTTTGGGTGCCAGTATGGGCATGAATGGTCGTCAGCGTGCCGACTGCTATCCCGAAAATGTCATGCAAGGCTTGGGCCATCCGCGCAAGACAGTTA‑GTGGTGCATGACGCCACGGAAACAATGGTGTCGGTGCCATCCAGAGTGTCGTCATTGACGTTATAAACGATAGtttt        <  2:324314/249‑1 (MQ=255)                                                                                                                                                                                                                                            tAGTGTAGGGAATGATATTTTCTGCCGCTGCGCGTGCAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTA‑GTGGTGCATAACGCCACGGAACCAATGGGGTCGTTGCCATCCAGAGTGGCTTCATTGACGTTATAAACGTTAGTtttcatt    >  2:20047/1‑251 (MQ=255)                                                                                                                                                                                                                                              gtgtAGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTA‑GTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTTATAAACGATAGTtttcatttc  <  1:313573/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        GAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATATTTTCTGCCGCTGCGCGTGAAGCACGTAAATCTTTACCACGCGGGCCATCCACCAGTGACTGGGTGCCAGTATAGGCATGAATGGTCGTCATCGTGCCGACTTCTATCCCGAAACTGTCATGCAAGGCTTTGGCCATCGGCGCAAGACAGTTA‑GTGGTGCATGACGCCACGGAAACAATGGTGTCGTTGCCATCCAGAGTGTCGTCATTGACGTTATAAACGATAGTTTTCATTTC  >  NC_000913/1489854‑1490341

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