Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,811,682	G→A	14.3%	M213I (ATG→ATA)	dsbG →	protein disulfide isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,811,682	0	G	A	14.3%	61.9 / 4.2	28	M213I (ATG→ATA)	dsbG	protein disulfide isomerase
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base A (2/2);  total (14/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 = 5.58e-01

GGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTGGGAAGATTCTGGGCAAGCGCTAGTGCATGTGGTGGGGCACAG  >  NC_002947/4811536‑4811821                                                                                                                                                   |                                                                                                                                            ggCAAGTCGGCGGCGCTGCTGACGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGCACATGGCGTTGATGGaa                                                                                                                                          <  2:439331/150‑1 (MQ=255)       tCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATgg                                                                                                                                            >  1:474781/1‑142 (MQ=255)       tCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATgg                                                                                                                                            <  2:474781/142‑1 (MQ=255)        gggcggcGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATgg                                                                                                                                   <  2:826205/149‑1 (MQ=255)                         gcgAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGCAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCg                                                                                                                  <  2:547075/150‑1 (MQ=255)                                 cccGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGCTGGAAGAGATGGGGTTGCAGGCGACTCCGGCCAgttt                                                                                                          >  1:271566/1‑150 (MQ=255)                                                  agcagcaTGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGAATACGGCCATTTTCTATCAGGATgttcaat                                                                                         <  2:189102/150‑7 (MQ=255)                                                        aTGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAgg                                                                                        <  2:4795/145‑1 (MQ=255)                                                        aTGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAgg                                                                                        >  1:4795/1‑145 (MQ=255)                                                        aTGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAAcc                                                                                   >  2:207891/1‑150 (MQ=255)                                                              cctccGGCGAGGCGACGACGGTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCaaa                                                                             <  2:239538/147‑1 (MQ=255)                                                                    gCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGc                                                                                            >  1:454404/1‑129 (MQ=255)                                                                    gCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGc                                                                                            <  2:454404/129‑1 (MQ=255)                                                                           cAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGc                                                                >  2:497479/1‑150 (MQ=255)                                                                                   tGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATAGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCCCGGc                                                        >  1:723728/1‑150 (MQ=255)                                                                                          aCTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTa                                                                                                       <  1:667858/96‑1 (MQ=255)                                                                                          aCTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTa                                                                                                       >  2:667858/1‑96 (MQ=255)                                                                                              gAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCtt                                             <  1:652412/150‑1 (MQ=255)                                                                                                gCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTgg                                           >  2:258066/1‑150 (MQ=255)                                                                                                gCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATAGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTgg                                           <  2:723728/150‑1 (MQ=255)                                                                                                    gTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGCAACAGCCCTTAGATGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCCCGGCCGGAGATGCTTGGGAAg                                       >  1:480562/1‑150 (MQ=255)                                                                                                          gagGCTGTGCAACATAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAg                                                                                         <  1:362499/94‑1 (MQ=255)                                                                                                          gagGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAg                                                                                         >  2:362499/1‑94 (MQ=255)                                                                                                                           gCTGGCGGCGAACATGGCGTTGATAGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGga                                                    <  1:159374/114‑1 (MQ=255)                                                                                                                           gCTGGCGGCGAACATGGCGTTGATAGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGga                                                    >  2:159374/1‑114 (MQ=255)                                                                                                                             tGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTGGGAAGATTCTGGGCAAGCGCTAGTGCAtgt              >  2:690855/1‑150 (MQ=255)                                                                                                                             tGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTGGGAAGATTCTGGGCAAGCGCTAGTGCAtgt              <  1:690855/150‑1 (MQ=255)                                                                                                                                   cgAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTGGGAAGATTCTGGGCAAGCGCTAGTGCATGTGGTggg        <  1:497479/150‑1 (MQ=255)                                                                                                                                    gAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTGGGAAGATTCTGGGCAAGCGCTAGTGCATGTGGTgggg       >  2:35732/1‑150 (MQ=255)                                                                                                                                         tGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTGGGAAGATTCTGGGCAAGCGCTAGTGCATGTGGTGGGGCACAg  >  1:252974/1‑150 (MQ=255)                                                                                                                                                   ggAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTGGGaa                                        <  2:508909/102‑1 (MQ=255)                                                                                                                                                   ggAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTGGGaa                                        >  1:508909/1‑102 (MQ=255)                                                                                                                                                   |                                                                                                                                            GGCAAGTCGGCGGCGCTGCTGGCGGCGAAAGACCCGGTCAAGGCCCTGCAGCAGCATGAAAAGGCCGGCAAGGCCAGCACGCTGAAGGCACTCGAGCAGGTGCCAGAGGCTGTGCAACAGAAGCTGGCGGCGAACATGGCGTTGATGGAAGAGATGGGGTTGCAGGCGACTCCGGCCATTTTCTATCAGGATGAGCAAGGCAACCTGCAAAGCCAGCAAGGGGCACCACGGCCGGAGATGCTTGGGAAGATTCTGGGCAAGCGCTAGTGCATGTGGTGGGGCACAG  >  NC_002947/4811536‑4811821

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: ‑