Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	367,367	T→G	16.0%	intergenic (‑53/+21)	PP_5442 ← / ← PP_0306	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	367,367	0	T	G	16.0%	87.0 / 15.5	44	intergenic (‑53/+21)	PP_5442/PP_0306	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/16);  new base G (3/4);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

CCCATGTATCGAGTGTCAGCCAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGAAGCGCGATAGACGCGGGTAGAGG  >  NC_002947/367230‑367509                                                                                                                                           |                                                                                                                                               cccATGTATCGAGTGTCAGCCAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGAGCGCGCTGGTAccc                                                                                                                                       <  3:284230/147‑1 (MQ=255)   cATGTATCGAGTGTCAGCCAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGcc                                                                                                                                   >  5:202213/1‑149 (MQ=255)   cATGTATCGAGTGTCAGCCAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGAGCGCGCTGGTACCCg                                                                                                                                      <  6:229868/146‑1 (MQ=255)   cATGTATCGAGTGTCAGCCAGGGCGCAGGCTTGATGCTTAGCAGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGAGCGCGCTGGTACCCg                                                                                                                                      >  5:229868/1‑146 (MQ=255)    aTGTATCGAGTGTCAGCCAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCt                                                                                                                                  >  8:99517/1‑149 (MQ=255)         tCGAGTGTCAGCCAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCt                                                                                                                              >  7:256738/1‑148 (MQ=255)         tCGAGTGTCAGCCAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCt                                                                                                                              >  5:139028/1‑148 (MQ=255)                   gCCAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATTCATATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGc                                                                                                                   <  1:88983/149‑1 (MQ=255)                    ccAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGAGCGCGCTGGTAccc                                                                                                                                       <  4:21699/128‑1 (MQ=255)                    ccAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGAGCGCGCTGGTAccc                                                                                                                                       >  3:21699/1‑128 (MQ=255)                        ggCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCgaga                                                                                                               >  2:44053/1‑148 (MQ=255)                         gcgcAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGaaa                                                                                                             <  6:139028/149‑1 (MQ=255)                         gcgcAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGaaa                                                                                                             >  5:190981/1‑149 (MQ=255)                         gcgcAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGaaa                                                                                                             <  3:71028/149‑1 (MQ=255)                                cTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGcc                                                                                                       >  2:161989/1‑148 (MQ=255)                                cTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGcc                                                                                                       >  4:15699/1‑148 (MQ=255)                                   gATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAg                                                                                                    <  7:99517/148‑1 (MQ=255)                                   gATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGAGCGCGCTGGTACCCgg                                                                                                                                     >  6:142291/1‑115 (MQ=255)                                   gATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGAGCGCGCTGGTACCCgg                                                                                                                                     <  5:142291/115‑1 (MQ=255)                                       cTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGc                                                                                               <  2:147050/149‑1 (MQ=255)                                          aGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCt                                                                                             <  1:44053/148‑1 (MQ=255)                                            ccGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTg                                                                                          <  6:190981/149‑1 (MQ=255)                                                   tcGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGcc                                                                                                       <  7:116597/129‑1 (MQ=255)                                                   tcGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGcc                                                                                                       >  8:116597/1‑129 (MQ=255)                                                   tcGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGAtcttc                                                                                    >  5:6499/1‑148 (MQ=255)                                                   tcGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATAGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGAtcttct                                                                                   >  1:229431/1‑149 (MQ=255)                                                               gtgtCCAGGCTCGTATCGTCCATGTGGCATTGCACGATTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACAAGGCCTCGCTAATCAGATAGCGAGAAAAGCGCCGAGGTCGCCTTTGATATTATAATTCACAAgt                                                                        >  1:216299/1‑148 (MQ=255)                                                                            tATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTc                                                          >  2:248555/1‑149 (MQ=255)                                                                             aTCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTc                                                          >  6:142891/1‑148 (MQ=255)                                                                                                          gatgCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCttgtt                             <  1:161989/149‑1 (MQ=255)                                                                                                          gatgCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCttgtt                             <  3:15699/149‑1 (MQ=255)                                                                                                               aGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAg                                                                         <  1:209386/100‑1 (MQ=255)                                                                                                               aGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAg                                                                         >  2:209386/1‑100 (MQ=255)                                                                                                               aGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGaa                        <  1:248555/149‑1 (MQ=255)                                                                                                               aGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGaa                        <  6:202213/149‑1 (MQ=255)                                                                                                               aGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGaa                        <  2:229431/149‑1 (MQ=255)                                                                                                                     ccAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGAAGcgcg                   >  4:50538/1‑148 (MQ=255)                                                                                                                      cAGAACCGGGTACCAGCGCTCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGAAGCGCGAt                 <  2:216299/149‑1 (MQ=255)                                                                                                                        gAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGAAGCGCGa                  >  8:178274/1‑146 (MQ=255)                                                                                                                              ggTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGAAGCGCGATAGACGCgg         >  6:71594/1‑149 (MQ=255)                                                                                                                                   cAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGAtcttc                                                                                    <  1:121715/69‑1 (MQ=255)                                                                                                                                   cAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGAtcttc                                                                                    >  2:121715/1‑69 (MQ=255)                                                                                                                                    aGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGAAGCGCGATAGACGCGGGTagag   >  7:302120/1‑149 (MQ=255)                                                                                                                                     gcgcgcTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGAAGCGCGATAGACGCAGGTAGAgg  >  1:93144/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                               CCCATGTATCGAGTGTCAGCCAGGGCGCAGGCTTGATGCTTAGCCGCTCCTCGAATATTCGGGTGTCCAGGCTCGTATCGTCCATGTGGCATTGCTC‑CTTGATGGATGCAGATGACCAGAACCGGGTACCAGCGCGCTCTGGTACCCGGCCTCGCTAATCAGATAGCGAGAAACGCGCCGAGGTCGCCTTTGATCTTCTCATTCACAAGGTCACGTGCAGCGTCGAAGTGACGAACGCTGAACGTGCCTTGTTGCGAAGCGCGATAGACGCGGGTAGAGG  >  NC_002947/367230‑367509

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