Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,488,689	C→A	9.1%	R83S (CGT→AGT)	ybgI →	metal‑binding hydrolase/oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,488,689	0	C	A	9.1%	134.1 / 4.3	44	R83S (CGT→AGT)	ybgI	metal‑binding hydrolase/oxidase
Reads supporting (aligned to +/- strand):  ref base C (20/20);  new base A (2/2);  total (22/22)
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.49e-01

GGCCGCAGGTCAGCCGTATTGTCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCC  >  NC_002947/1488543‑1488835                                                                                                                                                     |                                                                                                                                                   ggccgCAGGTCAGCCGTATTGTCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGt                                                                                                                                                  >  6:242581/1‑149 (MQ=255)    cgCAGGTCAGCCGTATTGTCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGc                                                                                                                                               >  2:212739/1‑149 (MQ=255)      cAGGTCAGCCGTATTGTCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCt                                                                                                                                             <  5:148668/149‑1 (MQ=255)       aGGTCAGCCGTATTGTCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTg                                                                                                                                            <  6:9656/149‑1 (MQ=255)             gCCGTATTGTCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGAcc                                                                                                                                      >  6:102584/1‑149 (MQ=255)                gTATTGTCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTg                                                                                                                                   <  3:121916/149‑1 (MQ=255)                 tATTGTCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGc                                                                                                                                  >  6:36799/1‑149 (MQ=255)                   ttGTCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCt                                                                                                                                 <  4:24389/148‑1 (MQ=255)                      tCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGaa                                                                                                                              >  7:252510/1‑148 (MQ=255)                          tgttgtTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGaa                                                                                                                                          <  6:176082/129‑1 (MQ=255)                          tggtgTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGaa                                                                                                                                          >  5:176082/1‑132 (MQ=255)                                 aCCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAGTCAACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGatat                                                                                                                     <  1:212739/148‑1 (MQ=255)                                        gccagGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTg                                                                                                           <  5:242581/149‑1 (MQ=255)                                         ccagGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTgg                                                                                                          >  3:155805/1‑149 (MQ=255)                                           agGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCtt                                                                                                                                    <  7:82344/121‑1 (MQ=255)                                           agGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCtt                                                                                                                                    >  8:82344/1‑121 (MQ=255)                                            ggCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGt                                                                                                       >  6:169142/1‑149 (MQ=255)                                             gCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGAGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGtt                                                                                                      >  3:213141/1‑149 (MQ=255)                                             gCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGAGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGt                                                                                                       >  3:209544/1‑148 (MQ=255)                                                   cTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGTCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCt                                                                                                <  7:31643/149‑1 (MQ=255)                                                              ggTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACg                                                                                     >  1:8366/1‑149 (MQ=255)                                                                           gccgaTCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGt                                                                        <  4:155805/149‑1 (MQ=255)                                                                             cgaTCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGc                                                                                                                                  <  7:87643/89‑1 (MQ=255)                                                                             cgaTCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGc                                                                                                                                  >  8:87643/1‑89 (MQ=255)                                                                              gaTCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTggg                                                                     >  5:167187/1‑149 (MQ=255)                                                                                    gtgctggtgcACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTgg                                                                      <  8:252510/142‑1 (MQ=255)                                                                                    gtgctggtgcACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAAAAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTggg                                                                     <  5:260745/143‑1 (MQ=255)                                                                                       ctggtgcACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAAtgt                                                            <  3:223189/149‑1 (MQ=255)                                                                                          gtgcACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGc                                                          <  6:167187/148‑1 (MQ=255)                                                                                              accacGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGcaacaa                                                               >  5:12233/1‑139 (MQ=255)                                                                                              accacGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGcaacaa                                                               <  6:12233/139‑1 (MQ=255)                                                                                                             ggAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTg                                         >  1:36129/1‑146 (MQ=255)                                                                                                              gAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGaa                                       <  2:8366/147‑1 (MQ=255)                                                                                                              gAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGAGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGaa                                       <  4:213141/147‑1 (MQ=255)                                                                                                              gAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGAGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGaa                                       <  4:209544/147‑1 (MQ=255)                                                                                                                               gTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGc                      >  4:187865/1‑149 (MQ=255)                                                                                                                                  cATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGc                   <  5:102584/149‑1 (MQ=255)                                                                                                                                   aTCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTCACCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGTGCCGc                   >  1:181309/1‑148 (MQ=255)                                                                                                                                    tCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGc                   >  1:68677/1‑147 (MQ=255)                                                                                                                                             ttCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCg                            >  7:73191/1‑129 (MQ=255)                                                                                                                                             ttCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCg                            <  8:73191/129‑1 (MQ=255)                                                                                                                                                 ccAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAAcc    >  2:166784/1‑149 (MQ=255)                                                                                                                                                   aGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAAcccc  <  5:36799/149‑1 (MQ=255)                                                                                                                                                   aGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAAccc   >  5:1097/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGCCGCAGGTCAGCCGTATTGTCAGTGGTGTTACCGCCAGCCAGGCATTGCTGGATGCCGCGGTCGAGGCCGAGGCCGATCTGGTGCTGGTGCACCACGGTTATTTCTGGAAGGGCGAG‑‑AACCCGTGCATCACCGGTATTCGCCAGCGTCGCCTGAAGACCTTGCTGAATAACGATATAAGCCTGTTGGCGTTCCACCTGCCCTTGGACGTGCACCCGGACGTGGGCAACAATGTGCAGCTGGCGCGGCGACTGAAGATCACTGTCGAAGGGCCGCTAGACCCGGAAAACCCC  >  NC_002947/1488543‑1488835

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