Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,595,188	T→G	11.1%	intergenic (+47/‑94)	PP_5511 → / → PP_2273	hypothetical protein/DNA polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,595,188	0	T	G	11.1%	147.6 / 10.7	54	intergenic (+47/‑94)	PP_5511/PP_2273	hypothetical protein/DNA polymerase I
Reads supporting (aligned to +/- strand):  ref base T (22/26);  new base G (4/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

CTGGGCGAAGCCCTCAAGGGCCTGTACCACGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAGGAGCAACATTCCCAAGGCCGCTTCATCGGCGGC  >  NC_002947/2595040‑2595335                                                                                                                                                     |                                                                                                                                                    cTGGGCGAAGCCCTCAAGGGCCTGTACCACGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAAt                                                                                                                                                     <  8:182485/149‑1 (MQ=255)    ggCGAAGCCCTCAAGGGCCTGTACCACGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATgg                                                                                                                                                   <  7:65770/148‑1 (MQ=255)       gAAGCCCTCAAGGGCCTGTACCACGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTc                                                                                                                                               >  1:308702/1‑149 (MQ=255)        aaGCCCTCAAGGGCCTGTACCACGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTcc                                                                                                                                              <  6:37717/149‑1 (MQ=255)         aGCCCTCAAGGGCCTGTACCACGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCATTCGTTTCGAGGTCGAACAGAAGGGGTTctt                                                                                                                                             <  7:192059/149‑3 (MQ=38)                   ggCCTGTACCACGGCTTTTGTGATGCCCTGGCCGATACCGGCGCGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTg                                                                                                                                            >  1:242616/1‑140 (MQ=255)                   ggCCTGTACCACGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTg                                                                                                                                            <  2:242616/140‑1 (MQ=255)                         tACCACGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCg                                                                                                                               <  6:207724/147‑1 (MQ=255)                            cACGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATc                                                                                                                          >  4:344525/1‑149 (MQ=255)                             aCGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATc                                                                                                                          >  7:209039/1‑148 (MQ=255)                                gCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGAGTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCAtt                                                                                                                      <  1:110108/149‑1 (MQ=255)                                gCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCAtt                                                                                                                      <  8:209039/149‑1 (MQ=255)                                      gtgATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCtt                                                                                                                 <  3:344525/148‑1 (MQ=255)                                             ccTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCg                                                                                                          >  8:113038/1‑148 (MQ=255)                                                       aCCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCa                                                                                                                            >  3:184150/1‑120 (MQ=255)                                                       aCCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCa                                                                                                                            <  4:184150/120‑1 (MQ=255)                                                         cGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAAc                                                                                             >  4:37976/1‑149 (MQ=255)                                                             aCGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCGTTGTCATCGAAGACACGGCAACTggg                                                                                         <  2:308702/149‑1 (MQ=255)                                                               gagctTCCGAGAGCTGTGGCATACGCCTACTCCCTGCTTACTGGGGTACCAATGCATGAAGGCTTTGTGTTGGACCTCGCATCGAATGGGTTACTGGAGGAAGTCTCAGTGATGCATTGCCTTGTCATCGAAGACACGGCAACTGGGGa                                                                                       <  8:236127/146‑1 (MQ=255)                                                               gCTCTACCGAGAGCTGCGGCATACGCCGACCCCCTGCTGACTTGGTGACCAATGCAAGAAACACTTCTGTTCGACCTAGAAACGAATGGGTTCCTGGAGTGTCTCTCCGTCATCAATTGCCTTGTCATCGATGACACGGCAACTGGGGa                                                                                       <  8:271076/149‑1 (MQ=255)                                                                       gagaGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTgggg                                                                                        <  1:82782/140‑1 (MQ=255)                                                                       gagaGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTgggg                                                                                        >  2:82782/1‑140 (MQ=255)                                                                        agagCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGa                                                                               >  5:120300/1‑148 (MQ=255)                                                                            cTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTc                                                                          <  3:37976/149‑1 (MQ=255)                                                                             tgtgGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCATGGAGGCAGTCTCCGTAAACCATTGCCTTGTCAACGAAGCAACGCCAACTGGGGACGTTAAGAAATTcc                                                                         >  6:265885/1‑149 (MQ=255)                                                                                    tACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGa                                                                                                      <  6:218889/113‑1 (MQ=255)                                                                                    tACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGa                                                                                                      >  5:218889/1‑113 (MQ=255)                                                                                    tACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTgg                                                                  >  3:64925/1‑149 (MQ=255)                                                                                       gCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCt                                                               >  2:302146/1‑149 (MQ=255)                                                                                        ccGACCCCCTGCTTACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTc                                                              <  5:265885/149‑1 (MQ=255)                                                                                             ccccTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGc                                                         >  6:95003/1‑149 (MQ=255)                                                                                                   ctgACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATgg                                                    <  5:95003/148‑1 (MQ=255)                                                                                                    tgACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGa                                                  <  7:113038/149‑1 (MQ=255)                                                                                                     gACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAg                                                 >  6:153401/1‑149 (MQ=255)                                                                                                      aCTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGa                                                  <  6:120300/147‑1 (MQ=255)                                                                                                      aCTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGa                                                  <  1:302146/147‑1 (MQ=255)                                                                                                           gggACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATc                                                                                                           <  5:322785/85‑1 (MQ=255)                                                                                                           gggACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTACTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATc                                                                                                           >  6:322785/1‑85 (MQ=255)                                                                                                                 aaTGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATgg                                                    >  2:188718/1‑134 (MQ=255)                                                                                                                 aaTGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATgg                                                    <  1:188718/134‑1 (MQ=255)                                                                                                                  ctccaggAACCCATTCGTTTCGAGGTGGAACAGAAGGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTc                                      >  1:186539/7‑147 (MQ=37)                                                                                                                    ccaggAACCCATTCGTTTCGAGGTCGAACAGAAGGGGTTCTTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAg                                   >  1:12185/5‑148 (MQ=37)                                                                                                                    ccaggAACCCATTCGTTTCGAGGTCGAACAGAAGGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCa                                                       >  1:110885/5‑128 (MQ=37)                                                                                                                    ccaggAACCCATTCGTTTCGAGGTCGAACAGAAGGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCa                                                       <  2:110885/124‑1 (MQ=37)                                                                                                                      aaGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAGGAg                                >  1:295610/1‑149 (MQ=255)                                                                                                                       ggAACCCATTCGTTTCGAGGTCGAACAGAAGGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGATCCAGGAg                                >  7:159880/2‑148 (MQ=38)                                                                                                                              cTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGa                                                  >  2:39291/1‑123 (MQ=255)                                                                                                                              cTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGa                                                  <  1:39291/123‑1 (MQ=255)                                                                                                                                          cTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAGGAGCAACATTCCCAAGGCCGCt             >  4:333351/1‑148 (MQ=255)                                                                                                                                          cTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAGGAGCAACATTCCCAAGGCCGCt             <  3:333351/148‑1 (MQ=255)                                                                                                                                           tCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAGGAGCAACATTCCCAAGGCCGc              >  2:80300/1‑146 (MQ=255)                                                                                                                                           tCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAGGAGCAACATTCCCAAGGCCGc              <  1:80300/146‑1 (MQ=255)                                                                                                                                            cGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAGGAGCAACATTCCCAAGGCCGCTTCa          >  1:306653/1‑149 (MQ=255)                                                                                                                                              aaaCGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAGGAGCAACATTCCCAAGGCCGCTTCATc        <  4:64925/149‑1 (MQ=255)                                                                                                                                                    aTGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAGGAGCAACATTCCCAAGGCCGCTTCATCggcggc  <  5:153401/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CTGGGCGAAGCCCTCAAGGGCCTGTACCACGGCTTTTGTGATGCCCTGGCCGATACCGGCACGCTCTTCCGAGAGCTGTGGCATACGCCGACCCCCTGCTGACTGGGGGACCAATGCAAGAACCCCTTCTGTTCGACCTCGAAACGAATGGGTTCCTGGAGGCAGTCTCCGTCATCCATTGCCTTGTCATCGAAGACACGGCAACTGGGGACGTTAAGAAATTCCCGCCTGGCCTCATCGCCATGGGAGTCAAATGGCTCCAGGAGCAACATTCCCAAGGCCGCTTCATCGGCGGC  >  NC_002947/2595040‑2595335

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