Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,418,522	T→C	9.6%	intergenic (‑111/+48)	prtN ← / ← PP_3032	transcriptional regulator PrtN/DNA‑binding protein Roi‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,418,522	0	T	C	9.6%	108.6 / 3.4	43	intergenic (‑111/+48)	prtN/PP_3032	transcriptional regulator PrtN/DNA‑binding protein Roi‑like protein
Reads supporting (aligned to +/- strand):  ref base T (20/18);  new base C (2/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

CGCAGCGCGATTTGATTGGTATTGCTCATCAGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGGGGATCCCGGCCTGCTTCACGCGTGTCGATTCGC  >  NC_002947/3418383‑3418670                                                                                                                                            |                                                                                                                                                     cgcAGCGCGATTTGATTGGTATTGCTCATCAGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGctcct                                                                                                                                             <  6:112519/149‑1 (MQ=255)  gctgcgcGAGTTGATTGGTTTTGCTCATCAGGCTTCGGGCTCTCGCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAACGCTGCGctccag                                                                                                                                            <  4:208938/146‑3 (MQ=255)  gcAGCGCGATTTGATTGGTATTGCTCATCAGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAACGCTGCGctccag                                                                                                                                            <  6:190263/149‑3 (MQ=255)               aTTGGTATTGCTCATCAGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAAc                                                                                                                               >  8:39552/1‑149 (MQ=255)                  ggTATTGGTGATCAGGCTTCGGGCTCTGCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTg                                                                                                                             <  4:226131/148‑1 (MQ=255)                   gTATTGCTCATCAGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTgg                                                                                                                            <  8:285186/148‑1 (MQ=255)                   gTATTGCTCATCAGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTgg                                                                                                                            <  5:241878/148‑1 (MQ=255)                   gTATTGCTCATCAGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTg                                                                                                                             >  4:91184/1‑147 (MQ=255)                        gCTCATCAGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATc                                                                                                                      >  2:70569/1‑149 (MQ=255)                          tcatcaGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTg                                                                                                                    >  6:178699/1‑149 (MQ=255)                            atcaGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTggg                                                                                                                  >  2:225593/1‑149 (MQ=255)                            atcaGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTgg                                                                                                                   <  2:222067/148‑1 (MQ=255)                                   ttCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGGACTGGAGGATCTGGGCCGCGTa                                                                                                           >  5:281143/1‑149 (MQ=255)                                      gggCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACt                                                                                                         <  2:156608/148‑1 (MQ=255)                                      cgttcttcCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGCATCTGGGCCGCGTACTg                                                                                                        <  1:308673/143‑1 (MQ=255)                                             cccTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTg                                                                                                        >  2:308673/1‑142 (MQ=255)                                                          gCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCCCCATTCGTCAGGAACCTCCTGGAGCGCAGCGATGCGCTCCTCACGGGTTTTGAACTGCAGGATCTCGGCCGCGTACTGCCTAGGCCGCATAGCGATTg                                                                                    >  7:129624/1‑149 (MQ=255)                                                             tGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGt                                                                                                            <  6:65622/122‑1 (MQ=255)                                                                    tGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGACCCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCggtg                                                                          <  5:229381/149‑1 (MQ=255)                                                                          cAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCg                                                                    <  5:178699/149‑1 (MQ=255)                                                                                aTCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGt                                                               <  8:193235/148‑1 (MQ=255)                                                                                   tcGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGATTGTCTTCCGGTGGTGTCGGAAGTTg                                                             >  2:187072/1‑147 (MQ=255)                                                                                    cGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCGCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAg                                                          <  2:212177/149‑1 (MQ=255)                                                                                     gCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTg                                                             <  3:91184/145‑1 (MQ=255)                                                                                     gCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAgg                                                         >  4:74618/1‑149 (MQ=255)                                                                                           cTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGGACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCgag                                                    <  6:281143/148‑1 (MQ=255)                                                                                            ttGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGAc                                                  >  4:158110/1‑149 (MQ=255)                                                                                              gCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGt                                                >  6:139642/1‑149 (MQ=255)                                                                                                  aCCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTg                                            >  4:112457/1‑149 (MQ=255)                                                                                                   ccAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCggtg                                                                          >  8:122692/1‑118 (MQ=255)                                                                                                   ccAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCggtg                                                                          <  7:122692/118‑1 (MQ=255)                                                                                                    cAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGCTCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAACTTGCAGGTCGAGACGTTCTGc                                           >  7:229528/1‑148 (MQ=255)                                                                                                     aGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTa                                         <  1:225593/149‑1 (MQ=255)                                                                                                                   ttcAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGGGGATCCCg                           >  2:142831/1‑149 (MQ=255)                                                                                                                      aGGAACCTCCTGGAGCGCATCGTTGCGCTCCTCACGTGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGCGGATCCCGGcc                        <  8:129624/149‑1 (MQ=255)                                                                                                                      aGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGGGGATCCCGGcc                        >  6:122107/1‑149 (MQ=255)                                                                                                                            ctcctGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGGGGATCCCg                           >  4:6319/1‑140 (MQ=255)                                                                                                                            ctcctGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGGGGATCCCg                           <  3:6319/140‑1 (MQ=255)                                                                                                                              tgaggAGCGCAACGCTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGGGGATCCCGGCCTGCTTCAc                >  5:320940/4‑149 (MQ=255)                                                                                                                                 ggAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGGGGATCCCGGCCTGCTTCAcgcg             >  7:132290/1‑149 (MQ=255)                                                                                                                                  gAGCGCAACGCTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGGGGATCCCGGCCTGCTTCACGCgt            >  2:257167/1‑149 (MQ=255)                                                                                                                                      gcAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGCCGTTCCGCTAGCCAGGGGATCCCGGCCTGCTTCACGCGTGTc         >  2:308823/1‑148 (MQ=255)                                                                                                                                      gcAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCAGCTAGCCAGGGGATCCCGGCCTGCTTCACGCGTGTc         >  2:304072/1‑148 (MQ=255)                                                                                                                                            ttGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCgg                                                                   >  4:282073/1‑84 (MQ=255)                                                                                                                                            ttGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGGGGATCCCGGCCTGCTTCACGCGTGTCGATTCGc  <  1:70569/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                     CGCAGCGCGATTTGATTGGTATTGCTCATCAGGCTTCGGGCTCTCCCTGAGCTTGTGGCGTGACGGATGGTGCCAAGTGATCTCGCAGTGCTTGCGTACCAGGTCGCGCCATTCTTCAGGAACCTCCTGGAGCGCAGCGTTGCGCTCCTCACGGGTTTTGAACTGGAGGATCTGGGCCGCGTACTGCCTAGGCCGCATAGCGCTTGTCTTCCGGTGGTGTCGGAAGTTGCAGGTCGAGACGTTCTGCTAGCCAGGGGATCCCGGCCTGCTTCACGCGTGTCGATTCGC  >  NC_002947/3418383‑3418670

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