Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	848,505	C→G	12.3%	R22P (CGG→CCG)	PP_0730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	848,505	0	C	G	12.3%	147.2 / 9.9	56	R22P (CGG→CCG)	PP_0730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base G (5/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAA  >  NC_002947/848366‑848645                                                                                                                                            |                                                                                                                                             gccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCt                                                                                                                                      <  1:196493/148‑1 (MQ=255)    gccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAACGGCTTGGCTGcca                                                                                                                                  <  1:238068/149‑2 (MQ=255)           acaTTATCGATACGCAAATAAGGATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAACGGCTTGGCTGccacg                                                                                                                                >  4:313256/3‑141 (MQ=255)           aGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTCACGTGGCAGCCAAACGGCTTGGCTGccacg                                                                                                                                <  3:313256/144‑4 (MQ=255)            gATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACAGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGcc                                                                                                                          <  1:141047/149‑1 (MQ=255)               tATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGc                                                                                                                                       >  2:51709/1‑133 (MQ=255)                        gCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCt                                                                                                               >  3:314542/1‑148 (MQ=255)                        gCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCt                                                                                                               >  6:66777/1‑148 (MQ=255)                        gCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTATGGCTGCCTCAAAAGCCAGCAGGCCGCtt                                                                                                              <  1:220356/149‑1 (MQ=255)                          aaaTAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTAc                                                                                                            >  1:151882/1‑149 (MQ=255)                                gCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCt                                                                                                      <  5:203248/149‑1 (MQ=255)                                 cATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTg                                                                                                     <  4:16076/149‑1 (MQ=255)                                   tCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTc                                                                                                   >  4:311628/1‑149 (MQ=255)                                      acaTATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGttt                                                                                                                                          >  3:5038/1‑107 (MQ=255)                                      acaTATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGttt                                                                                                                                          <  4:5038/107‑1 (MQ=255)                                      acaTATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTcgac                                                                                                >  1:202156/1‑149 (MQ=255)                                               aGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCt                                                                                                                                  <  7:103082/106‑1 (MQ=255)                                               aGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCt                                                                                                                                  >  8:103082/1‑106 (MQ=255)                                                    aGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGtt                                                                                  <  4:314542/149‑1 (MQ=255)                                                      gCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGttt                                                                                 >  5:7952/1‑148 (MQ=255)                                                        tCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAt                                                                              <  6:7952/149‑1 (MQ=255)                                                                      gtcgtcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGAt                                                                 >  8:42032/1‑148 (MQ=255)                                                                       tcgtcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTa                                                               >  5:313304/1‑149 (MQ=255)                                                                       tcgtcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTa                                                               >  4:99750/1‑149 (MQ=255)                                                                            cTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAgg                                                          <  6:313304/149‑1 (MQ=255)                                                                             ttGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAATGGAACCTCGGATTACCAgg                                                          >  3:325804/1‑148 (MQ=255)                                                                               gATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGAt                                                       >  7:225135/1‑149 (MQ=255)                                                                                 tCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATg                                                      >  7:174512/1‑148 (MQ=255)                                                                                 tCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATg                                                      <  7:42032/148‑1 (MQ=255)                                                                                  cGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGa                                                     >  6:243520/1‑148 (MQ=255)                                                                                   gagaTCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGaaa                                                   <  3:99750/149‑1 (MQ=255)                                                                                    agaTCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAg                                                  <  4:325804/149‑1 (MQ=255)                                                                                     gaTCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAg                                                  >  1:38046/1‑148 (MQ=255)                                                                                           cAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACgt                                                                                     <  5:377156/107‑1 (MQ=255)                                                                                           cAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACgt                                                                                     >  6:377156/1‑107 (MQ=255)                                                                                             ggAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTg                                         <  8:174512/149‑1 (MQ=255)                                                                                             ggAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTg                                         >  4:83702/1‑149 (MQ=255)                                                                                             ggAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTg                                         >  2:259204/1‑149 (MQ=255)                                                                                             ggAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTg                                         >  4:261116/1‑149 (MQ=255)                                                                                              gAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGc                                        >  3:140266/1‑149 (MQ=255)                                                                                                 cGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGa                                     <  4:140266/149‑1 (MQ=255)                                                                                                  gAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAg                                    <  3:221584/149‑1 (MQ=255)                                                                                                       gTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGa                               <  2:202156/149‑1 (MQ=255)                                                                                                          gcgcTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAAc                             <  5:66777/148‑1 (MQ=255)                                                                                                            gcTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAg                           >  8:338257/1‑148 (MQ=255)                                                                                                                    atgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTCTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                  <  8:59585/149‑1 (MQ=255)                                                                                                                       aaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGAATTGCGaagc                <  5:100328/148‑1 (MQ=255)                                                                                                                          gTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTTCGaagcaagc            >  7:162527/1‑149 (MQ=255)                                                                                                                          gTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaagcaagc            >  4:158110/1‑149 (MQ=255)                                                                                                                           ttttgaggCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCg           >  5:130118/7‑149 (MQ=255)                                                                                                                           ttttgaggCAGCCAAACGGCTTGGCTGCCACAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaagcaagc            >  7:114211/7‑148 (MQ=255)                                                                                                                           tGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCg           <  3:158110/149‑1 (MQ=255)                                                                                                                                aggCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGACCCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGACCAGCATTGCGAAGCAAGCGACCGt      >  7:278513/2‑149 (MQ=255)                                                                                                                                aggCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTCCAACCACGGATTACCAGGGATGAAAGTAAATAGTGGTGCGCCCGAACAGCATTGCGACGCAAGCGCCCGt      >  6:95760/2‑149 (MQ=255)                                                                                                                                  gCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTc     >  6:273619/1‑148 (MQ=255)                                                                                                                                    aGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGaa  <  8:225135/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                             GCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAA  >  NC_002947/848366‑848645

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