Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	848,507	T→C	12.7%	K21K (AAA→AAG)	PP_0730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	848,507	0	T	C	12.7%	137.0 / 9.2	55	K21K (AAA→AAG)	PP_0730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base C (5/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCA  >  NC_002947/848366‑848654                                                                                                                                              |                                                                                                                                                    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                                                                                                                        >  6:66777/1‑148 (MQ=255)                        gCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCt                                                                                                                        >  3:314542/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                                                                                              >  6:377156/1‑107 (MQ=255)                                                                                           cAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACgt                                                                                              <  5:377156/107‑1 (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)                                                                                                                                             gTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCa  <  3:311628/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    GCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCA  >  NC_002947/848366‑848654

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