Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,381,086	G→A	11.1%	A568A (GCC→GCT)	PP_3857 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,381,086	0	G	A	11.1%	93.5 / 2.7	36	A568A (GCC→GCT)	PP_3857	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/13);  new base A (2/2);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGTAGACGGAAGCGTTCGGCGCGGTACTCGTATTCGCCAATCTTCACGTAGGCTTCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGCCTGATAGTCGCTGACCCATTCGCGCATCGATGCTGCGTACGC  >  NC_002947/4380938‑4381226                                                                                                                                                     |                                                                                                                                              cgcgTAGACGGAAGCGTTCGGCGCGGTACTCGTATTCGCCAATCTTCACGTAGGCTTCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTgg                                                                                                                                               <  3:161295/149‑1 (MQ=255)      aGACGGAAGCGTTCGGCGCGGTACTCGTATTCGCCAATCTTCACGTAGGCTTCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAg                                                                                                                                          >  6:283535/1‑149 (MQ=255)          ggAAGCGTTCGGCGCGGTACTCGTATTCGCCAATCTTCACGTAGGCTTCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGAATCTGCAGCTGGGCCAGCGa                                                                                                                                       >  6:40158/1‑148 (MQ=255)                gTTCGGCGCGGTACTCGTATTCGCCAATCTTCACGTAGGCTTCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAg                                                                                                                                >  1:55987/1‑149 (MQ=255)                                      gCCAATCTTCACGTAGGCTTCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGAGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGt                                                                                                         >  8:293890/1‑149 (MQ=255)                                      gCCAATCTTCACGTAGGCTTCCGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGAGCCAGCGACACGAAGTTCATCGAGCCGCACCCA‑Cgct                                                                                                         >  6:13685/1‑147 (MQ=255)                                            cTTCACGTAGGCTGCGGCGTACGGGTGCCAGGTGTGGTGGTAGCAGCACTTATATTTCGGCGCTACCATAACGAACTGCGGATCAGTTTCGTGCATCTGCTGGTGGGCCAGCGACTCGAAGTTCTTCGAGGCGGCCCCA‑CGGTTCATc                                                                                                    <  4:196662/148‑1 (MQ=255)                                                aCGTAGGCTTCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCa                                                                                                                            <  3:210910/121‑1 (MQ=255)                                                aCGTAGGCTTCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCa                                                                                                                            >  4:210910/1‑121 (MQ=255)                                                    aGGCTTCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAgg                                                                                           >  7:26782/1‑149 (MQ=255)                                                        ttCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAg                                                                                            >  2:253959/1‑144 (MQ=255)                                                        ttCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAg                                                                                            <  1:253959/144‑1 (MQ=255)                                                           ggCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTACGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCAcc                                                                                    <  5:283535/149‑1 (MQ=255)                                                              gTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAg                                                                                 <  2:55987/149‑1 (MQ=255)                                                               tACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAgg                                                                                >  2:195415/1‑149 (MQ=255)                                                                    gTGCCAGGTGTCTTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATg                                                                           <  5:40158/149‑1 (MQ=255)                                                                     tGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGGTGCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTAATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATgg                                                                          <  5:333807/149‑1 (MQ=255)                                                                       ccAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACa                                                                                             >  3:13178/1‑128 (MQ=255)                                                                       ccAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACa                                                                                             <  4:13178/128‑1 (MQ=255)                                                                       ccAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGa                                                                        >  6:12790/1‑149 (MQ=255)                                                                        cAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGc                                                                                                 <  5:76891/123‑1 (MQ=255)                                                                        cAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGc                                                                                                 >  6:76891/1‑123 (MQ=255)                                                                                gTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCTCCAAGGCAATGGGATCGCAGCAt                                                               >  7:281828/1‑149 (MQ=255)                                                                                            cTGATATTTCGTCGCTACCAGAACGAACTGCTGATCTGTTTCGTGCATCTGCAGCTTAGCCAGCTACACGAAGTTCATCTAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCTGCATCCATAGCGGTaa                                                   <  5:13685/149‑1 (MQ=255)                                                                                            cTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGAGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTaa                                                   <  7:293890/149‑1 (MQ=255)                                                                                               atatTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAAcc                                                 <  5:12790/148‑1 (MQ=255)                                                                                               atatTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCg                                                >  4:269772/1‑149 (MQ=255)                                                                                                      ggCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGc                                            >  4:154952/1‑146 (MQ=255)                                                                                                      ggCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGc                                            <  3:154952/146‑1 (MQ=255)                                                                                                           tACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGCCTGATa                                      >  5:53530/1‑147 (MQ=255)                                                                                                             ccAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGCCTGATa                                      <  6:53530/145‑1 (MQ=255)                                                                                                                    gaacTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGCCTGATAGTCGCTGAcc                            >  4:109724/1‑148 (MQ=255)                                                                                                                     aacTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGCCTGATAGTCGCTGACCCa                          >  1:56442/1‑149 (MQ=255)                                                                                                                             aTCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGCCTGATAGTCGCTGACCCATTCGCGCa                  >  1:127504/1‑149 (MQ=255)                                                                                                                                    tCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGCCTGATAGTCGCTGACCCATTCGCGCATCGAtg            <  1:195415/148‑1 (MQ=255)                                                                                                                                       tGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCACCGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGCCTGATAGTCGCTGACCCATTCGCGCATCGAtgctg         >  1:78322/1‑149 (MQ=255)                                                                                                                                             tGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGCCTGATAGTCGCTGACCCATTCGCGCATCGATGCTGCGTAcgc  >  5:161638/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                              CGCGTAGACGGAAGCGTTCGGCGCGGTACTCGTATTCGCCAATCTTCACGTAGGCTTCGGCGTACGGGTGCCAGGTGTCGTGGTAGCAGTACTGATATTTCGGCGCTACCAGAACGAACTGCGGATCAGTTTCGTGCATCTGCAGCTGGGCCAGCGACACGAAGTTCATCGAGGCGGCCCCA‑CGGTTCATCTGCGACAGGTAGCACCAAGGCAATGGGATCGCAGCATCCATAGCGGTAACCGCTGCCTGATAGTCGCTGACCCATTCGCGCATCGATGCTGCGTACGC  >  NC_002947/4380938‑4381226

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