Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,405,327	A→G	9.1%	P78P (CCT→CCC)	PP_3886 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,405,327	0	A	G	9.1%	116.7 / 2.9	44	P78P (CCT→CCC)	PP_3886	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/19);  new base G (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGG  >  NC_002947/4405179‑4405471                                                                                                                                                     |                                                                                                                                                  ccAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGa                                                                                                                                                   <  7:65650/149‑1 (MQ=255)  cAGCGATCTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCTTGACGTAATTTGGAGGGACTGGGCATAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGTCATGTGCCGa                                                                                                                                                   <  7:121851/148‑1 (MQ=255)      gAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCtt                                                                                                                                             <  5:268768/149‑1 (MQ=255)      gAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCtt                                                                                                                                             <  7:12094/149‑1 (MQ=255)            ccTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGG‑GCtttttcg                                                                                                                                        <  2:179499/148‑6 (MQ=255)             cTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACTGAATTTGGAGGGACTGGGCTGAAATCTGAGCATCTAACGTCTCGACGCTTTCGATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATgg                                                                                                                                      <  5:252214/149‑1 (MQ=255)             cTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATgg                                                                                                                                      <  2:72538/149‑1 (MQ=255)             cTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATgg                                                                                                                                      >  7:37541/1‑149 (MQ=255)                ggAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTg                                                                                                                                   <  6:66997/149‑1 (MQ=255)                     aGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGAt                                                                                                                              <  6:37263/149‑1 (MQ=255)                     aGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGAt                                                                                                                              >  3:286273/1‑149 (MQ=255)                      gggCGACCCGACTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATg                                                                                                                             <  2:159251/149‑1 (MQ=255)                       ggCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGTGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGG‑Gc                                                                                                                                               >  6:330913/1‑130 (MQ=255)                       ggCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGTGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGG‑Gc                                                                                                                                               <  5:330913/130‑1 (MQ=255)                          gACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCaa                                                                                                                         >  8:9514/1‑149 (MQ=255)                           aCCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAAc                                                                                                                        <  7:9514/149‑1 (MQ=255)                              cGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCt                                                                                                                     >  3:300980/1‑149 (MQ=255)                                        cgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGc                                                                                                           <  3:114096/149‑1 (MQ=255)                                        cgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGc                                                                                                           >  4:114096/1‑149 (MQ=255)                                             aCCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGTGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTc                                                                                                       <  1:61854/149‑1 (MQ=255)                                              cccGTGACGGAATTTGCAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGt                                                                                                     >  1:302510/1‑149 (MQ=255)                                              accGTGACGGTATTTGGAGGGACTGGGCAGAAATCTTAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGt                                                                                                     <  4:286273/148‑1 (MQ=255)                                                 gTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGgc                                                                                                  >  2:3387/1‑149 (MQ=255)                                                      ggAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCACCCCTGTCACGTCGCTGTCGAGGCGCTTc                                                                                             >  8:237639/1‑149 (MQ=255)                                                      ggAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTc                                                                                             >  2:24787/1‑149 (MQ=255)                                                      ggAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCATGGAATGGGTGGAGATGGCAAACCTGTCAGGTCGCTGTCGTGGCGCTTc                                                                                             >  5:69736/1‑149 (MQ=255)                                                          tttGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                                         <  8:37541/149‑1 (MQ=255)                                                                                 tGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTAcc                                                                  >  2:57223/1‑149 (MQ=255)                                                                                          aaCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTc                                                                                    <  8:177239/122‑1 (MQ=255)                                                                                          aaCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTc                                                                                    >  7:177239/1‑122 (MQ=255)                                                                                          aaCGTCTCGAAGCTATCCATGACTGGCTACATTTAACGAAAAATCACCGGCATGTGCCGAG‑GCTTGGAATGTGTGGAGATGGAAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCATCCAGTCCCCACGCTGACTGTTAAACCTGCACGt                                                         >  6:181725/1‑149 (MQ=255)                                                                                           aCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGtt                                                        >  3:255281/1‑149 (MQ=255)                                                                                                         tCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCttt                                           <  1:24787/148‑1 (MQ=255)                                                                                                            aTGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatc                                       <  4:300980/149‑1 (MQ=255)                                                                                                            aTGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGat                                        >  6:266208/1‑148 (MQ=255)                                                                                                              gACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatcat                                     >  5:44699/1‑149 (MQ=255)                                                                                                                     tCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCaa                               >  4:100498/1‑148 (MQ=255)                                                                                                                     tCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCaa                               >  1:63715/1‑148 (MQ=255)                                                                                                                      ccATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCaaa                              >  2:204599/1‑148 (MQ=255)                                                                                                                        aTTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAAc                           <  4:255281/149‑1 (MQ=255)                                                                                                                               gAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg                    <  2:63715/149‑1 (MQ=255)                                                                                                                                   aaGCCTCGGCACATGCCGGG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGcc                  >  4:161714/1‑147 (MQ=255)                                                                                                                                       cccGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCGGcca            >  2:267606/1‑147 (MQ=255)                                                                                                                                       cccGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCACCAAAAACGACAAAGCCCACCTg            >  7:156451/1‑149 (MQ=255)                                                                                                                                        ccGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTgg           <  1:57223/149‑1 (MQ=255)                                                                                                                                                  ccGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCggg  <  1:131101/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAG‑GCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGGG  >  NC_002947/4405179‑4405471

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