Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,006,247	A→G	10.8%	A370A (GCA→GCG)	ydcJ →	metalloprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,006,247	0	A	G	10.8%	116.9 / 5.2	46	A370A (GCA→GCG)	ydcJ	metalloprotein
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (3/2);  total (23/23)
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.37e-01

GCGCTCCGGCCGAAGCCAATGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTTCCAATCCAACCTGGGCGACGATGC  >  NC_002947/6006105‑6006384                                                                                                                                               |                                                                                                                                          gcgcTCCGGCCGAAGCCAATGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGa                                                                                                                                      <  3:38050/148‑1 (MQ=255)  cgcTCCGGCCGAAGCCAATGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGAc                                                                                                                                     >  1:242474/1‑148 (MQ=255)    cTCCGGCCGAAGCCAATGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAg                                                                                                                                  <  2:285696/149‑1 (MQ=255)     tCCGGCCGAAGCCAATGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAg                                                                                                                                  >  3:71407/1‑148 (MQ=255)       cGGCCGAAGCCAATGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGTCGCGTGCGGCCAggcc                                                                                                                               <  2:213807/149‑3 (MQ=255)           cGAAGCCAATGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGgcc                                                                                                                           <  3:195480/149‑1 (MQ=255)                 cAATGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGtcctcgta                                                                                                                      >  2:15934/1‑143 (MQ=255)                  aaTGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGac                                                                                                                    <  5:195059/149‑1 (MQ=255)                    tGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGacca                                                                                                                  <  2:242474/149‑1 (MQ=255)                      ccGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGaccacc                                                                                                                >  5:175062/1‑149 (MQ=255)                       cGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCAccc                                                                                                               >  2:294649/1‑149 (MQ=255)                          gCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTg                                                                                                             <  2:147053/148‑1 (MQ=255)                             ttACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACTCGACCAGGAAGGCCGCCCGaa                                                                                                                    >  6:251722/1‑137 (MQ=255)                             ttACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGac                                                                                                                    <  5:251722/138‑1 (MQ=255)                               aCATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAgg                                                                                                        <  5:25819/148‑1 (MQ=255)                                 aTGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCtt                                                                                                     >  8:41056/1‑149 (MQ=255)                                   ggCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGa                                                                                                   >  8:241563/1‑149 (MQ=255)                                    gcgcTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGAt                                                                                                  <  7:285348/149‑1 (MQ=255)                                          ctACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGAcgc                                                                                            >  8:2607/1‑149 (MQ=255)                                                gCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCa                                                                                      >  4:223147/1‑149 (MQ=255)                                                          ccGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTc                                                                            <  3:223147/149‑1 (MQ=255)                                                            gAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGa                                                                          <  4:71407/149‑1 (MQ=255)                                                                tCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGgcg                                                                      >  7:122692/1‑149 (MQ=255)                                                                  ccGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGCCCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGCCGCGGGCCACGTGCACTTCGAGGCGCt                                                                    >  4:290878/1‑149 (MQ=255)                                                                  ccGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCt                                                                    >  6:139194/1‑149 (MQ=255)                                                                              gcgcAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGa                                                        >  1:228094/1‑149 (MQ=255)                                                                                  aaaTGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTc                                                    <  5:139194/149‑1 (MQ=255)                                                                                    aTGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCt                                                  <  1:14574/149‑1 (MQ=255)                                                                                                ggCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCgg                                             <  4:1258/142‑1 (MQ=255)                                                                                                ggCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCgg                                             >  3:1258/1‑142 (MQ=255)                                                                                                 gCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTgcg                                     <  7:241563/149‑1 (MQ=255)                                                                                                         tatTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGa                              >  3:193699/1‑148 (MQ=255)                                                                                                             tCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGCTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTtc                         >  1:250230/1‑149 (MQ=255)                                                                                                              ccGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTtc                         <  2:133940/148‑1 (MQ=255)                                                                                                              ccGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTtc                         <  8:170825/148‑1 (MQ=255)                                                                                                                     ttCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTtccaatccaa                 >  6:2985/1‑149 (MQ=255)                                                                                                                          cACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTTCCAATCCAACCTgg            <  7:2607/149‑1 (MQ=255)                                                                                                                           aCCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTTCCAATCCAACCTggg           <  4:193699/149‑1 (MQ=255)                                                                                                                              gagaAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTTCCAATCCAACCTGGGcga        >  5:242555/1‑149 (MQ=255)                                                                                                                              gagaAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCACGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTTCCAATCCAACCTGGGcga        <  8:122692/149‑1 (MQ=255)                                                                                                                                 gccttccTGGTCGCGTGCGGCCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTtccaa                      >  2:301664/6‑132 (MQ=255)                                                                                                                                   aGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTTCCAATCCAAc                >  2:43860/1‑136 (MQ=255)                                                                                                                                   aGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTTCCAATCCAAc                <  1:43860/136‑1 (MQ=255)                                                                                                                                    ttccTGGTCGCGTGCGGCCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTTCCAATCCAACCTGGGCGACGATgc  >  1:189092/3‑149 (MQ=255)                                                                                                                                       cTGGTCGCGTGCGGCCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTtccaa                      <  6:292292/126‑1 (MQ=255)                                                                                                                                       cTGGTCGCGTGCGGCCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTtccaa                      >  5:292292/1‑126 (MQ=255)                                                                                                                                               |                                                                                                                                          GCGCTCCGGCCGAAGCCAATGCCGAGCGTTACATGGCGCTGCTACAGGCAAACTTTGCCGAGTTCCCGGATGACCTGGCGCAAATGCGCGAACAGGGCCTCGCATATTTCCGCTACTTCGCCACCGAGAAAGGCCTGGCCGCACGCGACCAGGAAGGCCGCCCGACCACCCTGCAAGGCTTGATCGACGCGGGCCACGTGCACTTCGAGGCGCTGGTGTACGAGGACTTCCTGCCGGTGAGTGCGGCGGGGATCTTCCAATCCAACCTGGGCGACGATGC  >  NC_002947/6006105‑6006384

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