Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,879,894	T→G	11.4%	intergenic (+78/+58)	PP_1686 → / ← PP_1687	glutathione peroxidase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,879,894	0	T	G	11.4%	115.5 / 9.0	44	intergenic (+78/+58)	PP_1686/PP_1687	glutathione peroxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base G (2/3);  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 = 4.32e-01

CTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATCGCCAACT  >  NC_002947/1879751‑1880040                                                                                                                                                |                                                                                                                                                   cTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTcc                                                                                                                                               <  6:158876/149‑1 (MQ=255)    aCCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTg                                                                                                                                             >  8:108372/1‑148 (MQ=255)     ccAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTg                                                                                                                                             >  7:6346/1‑147 (MQ=255)       aaaCCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCAACTTCAGGAGAACGGCCTCCTGAAGt                                                                                                                                         <  5:138382/149‑1 (MQ=39)        aaCCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCAACTTCAGGAGAACGGCCTCCTGAAGtc                                                                                                                                        <  3:113567/149‑2 (MQ=39)          ccGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGc                                                                                                                                      <  7:108372/149‑1 (MQ=255)            ggATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCAGTTCTCCTGAAGTTGCgg                                                                                                                                    <  5:325530/149‑1 (MQ=255)                    ccGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatc                                                                                                                                >  1:313061/1‑145 (MQ=255)                     cgcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaagat                                                                                                                           <  2:249840/149‑1 (MQ=255)                     cgcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatc                                                                                                                                >  2:88902/1‑144 (MQ=255)                      gcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaaga                                                                                                                            <  5:148125/147‑1 (MQ=255)                      gcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaaga                                                                                                                            >  6:148125/1‑147 (MQ=255)                        gTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaagatcaa                                                                                                                        <  4:150037/149‑1 (MQ=255)                                    tATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtg                                                                                                            >  2:11481/1‑149 (MQ=255)                                         aaaaGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggctgg                                                                                                       >  7:260323/1‑149 (MQ=255)                                         aaaaGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggctgg                                                                                                       >  2:59120/1‑149 (MQ=255)                                         aaaaGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggctgg                                                                                                       <  1:59120/149‑1 (MQ=255)                                                      cTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACtt                                                                                          >  4:45219/1‑149 (MQ=255)                                                      cTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACt                                                                                           >  4:103057/1‑148 (MQ=255)                                                         acaGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCa                                                                                        <  1:284091/148‑1 (MQ=255)                                                                  aCCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGtgcatg                                                                               <  8:260323/148‑1 (MQ=255)                                                                      ggCTGCGCTTGCTGAGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACgg                                                                          <  3:103057/149‑1 (MQ=255)                                                                           cgcTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGtgaa                                                                      <  2:313061/148‑4 (MQ=255)                                                                           cgcTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTaa                                                                       <  3:164430/147‑1 (MQ=255)                                                                                  tgcGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGccacc                                                               >  7:5865/1‑148 (MQ=255)                                                                                    cggcgTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAAc                                                            <  1:11481/149‑1 (MQ=255)                                                                                        gTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGca                                                        <  3:45219/149‑1 (MQ=255)                                                                                         ttGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcac                                                       <  8:6346/149‑1 (MQ=255)                                                                                            cTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccac                                                    >  8:183300/1‑149 (MQ=255)                                                                                            cTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccac                                                    >  7:38888/1‑149 (MQ=255)                                                                                            cTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccac                                                    >  7:84965/1‑149 (MQ=255)                                                                                                aaaGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccaccac                                                 <  1:88902/148‑1 (MQ=255)                                                                                                     gcgAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCAccggccg                                           <  7:183300/149‑1 (MQ=255)                                                                                                                    cGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGcacaca                            >  8:218703/1‑149 (MQ=255)                                                                                                                            cccgcaaCTTCAAGAGAACGGCCTCCTGAAGTCGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATg                    >  3:81085/7‑149 (MQ=255)                                                                                                                               gcaaCTTCAGGAGAACGGCCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGac                                        <  4:148302/123‑1 (MQ=255)                                                                                                                               gcaaCTTCAGGAGAACGGCCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGac                                        >  3:148302/4‑126 (MQ=255)                                                                                                                                    ttCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCgg             >  4:182760/1‑148 (MQ=255)                                                                                                                                            gCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGAccacca                                   <  4:269437/118‑1 (MQ=255)                                                                                                                                            gCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGAccacca                                   >  3:269437/1‑118 (MQ=255)                                                                                                                                            gCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATCGCCa     >  1:109766/1‑148 (MQ=255)                                                                                                                                               gTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCAccggccg                                           >  4:34489/1‑107 (MQ=255)                                                                                                                                               gTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCAccggccg                                           <  3:34489/107‑1 (MQ=255)                                                                                                                                               gTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTCGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCCGACCACCAGAACACACACAAATGCCACAGGCATCGCCAACt  >  8:81258/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                   CTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATCGCCAACT  >  NC_002947/1879751‑1880040

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