Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,210,754	T→G	7.4%	intergenic (+31/+38)	PP_3692 → / ← PP_5621	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,210,754	0	T	G	7.4%	161.0 / 4.2	54	intergenic (+31/+38)	PP_3692/PP_5621	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/23);  new base G (2/2);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GTGTTTGAGCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGC  >  NC_002947/4210607‑4210902                                                                                                                                                    |                                                                                                                                                     gtgtTTGAGCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTc                                                                                                                                                     <  2:76191/149‑1 (MQ=255)         gCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCgg                                                                                                                                             >  5:270332/1‑149 (MQ=255)              cccAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGATAAATGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGAATAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAAACGCGTCAAACTCGCTTCTg                                                                                                                                        >  4:234683/1‑149 (MQ=255)               ccAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTgg                                                                                                                                       >  5:232874/1‑149 (MQ=255)               ccAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTgg                                                                                                                                       >  5:38280/1‑149 (MQ=255)                     ggaggaGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTg                                                                                                                                        >  5:254474/1‑142 (MQ=255)                     ggaggaGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTg                                                                                                                                        <  6:254474/142‑1 (MQ=255)                      gaggagCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGc                                                                                                                                 >  7:84465/1‑148 (MQ=255)                             ttCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCt                                                                                                                                         >  6:211641/1‑133 (MQ=255)                             ttCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCt                                                                                                                                         <  5:211641/133‑1 (MQ=255)                               cccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCACCCTCGGTTCTGGCTTAGCCaaa                                                                                                                             <  7:47375/143‑1 (MQ=255)                               cccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCaaa                                                                                                                             >  8:47375/1‑143 (MQ=255)                                 cagcagTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGgcagc                                                                                                                     <  3:107041/149‑1 (MQ=255)                                           ttGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCCCAGAACCGAGGTTGACTCGGTTCTgagcccg                                                                                                                                  >  2:111835/1‑120 (MQ=38)                                           ttGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCCCAGAACCGAGGTTGACTCGGTTCTgagcccg                                                                                                                                  <  1:111835/126‑7 (MQ=38)                                                 cTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACg                                                                                                      >  5:313530/1‑148 (MQ=255)                                                       gTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACATCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAgcggc                                                                                               >  1:146672/1‑149 (MQ=255)                                                          ttCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGAt                                                                                            >  2:236148/1‑149 (MQ=255)                                                           tCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGAt                                                                                            <  6:232874/148‑1 (MQ=255)                                                             gAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTcc                                                                                         >  2:268596/1‑149 (MQ=255)                                                              aGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCg                                                                                        <  4:78222/149‑1 (MQ=255)                                                                 cTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGAt                                                                                      >  1:20586/1‑148 (MQ=255)                                                                  tGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGAtt                                                                                     >  4:289346/1‑148 (MQ=255)                                                                    aCGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTa                                                                                    <  6:313530/147‑1 (MQ=255)                                                                    aCGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTa                                                                                    <  1:268596/147‑1 (MQ=255)                                                                             cactgcaTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCaa                                                                           <  4:151249/147‑1 (MQ=255)                                                                                          gATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGgca                                                                                                                       <  4:107553/90‑1 (MQ=255)                                                                                          gATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGgca                                                                                                                       >  3:107553/1‑90 (MQ=255)                                                                                               aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCg                                                                                                                                              >  2:25354/1‑62 (MQ=255)                                                                                               aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGAt                                                                    <  5:199267/136‑1 (MQ=255)                                                                                               aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGAt                                                                    >  6:199267/1‑136 (MQ=255)                                                                                               aCGAATTCTCAGCACGCAAGTAGAAGCATACCGTGAGCGGGCTCAGAACCGAGTCAACCTCg                                                                                                                                              <  1:25354/62‑1 (MQ=255)                                                                                                cGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCgg                                                         >  1:126359/1‑146 (MQ=255)                                                                                                cGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCgg                                                         <  2:126359/146‑1 (MQ=255)                                                                                                 gAATTCTCAGAACGCAAGTAGAAGCAAACGGTGAGCGGGCTCAGAACCGAGTAAACCTCGGTTCTGCCTTAGCCAAACTGGCAGCATAGTCACCAAAACGTAGAGGCGATTCAGATTACCATCCCAATACCGATCAAGCTACCGGGGAt                                                     >  2:274354/1‑149 (MQ=255)                                                                                                  aaTTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGTGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCTATg                                                    <  8:84465/149‑1 (MQ=255)                                                                                                   aTTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGa                                                   <  6:38280/149‑1 (MQ=255)                                                                                                    ttCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAg                                                  <  3:289346/149‑1 (MQ=255)                                                                                                    ttCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAg                                                  <  2:146672/149‑1 (MQ=255)                                                                                                     tctcAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGTACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGc                                                 >  1:272596/1‑149 (MQ=255)                                                                                                      ctcAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCt                                                >  6:260479/1‑149 (MQ=255)                                                                                                       tcAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTc                                               >  2:98167/1‑149 (MQ=255)                                                                                                        cAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTcc                                              >  5:61931/1‑149 (MQ=255)                                                                                                           cACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCACACTGGCAGCATACTCAGCAAAACGTAGCGGCGATTCCGAGTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCg                                            >  8:181798/1‑148 (MQ=255)                                                                                                           cACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCATCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCg                                            >  8:177132/1‑148 (MQ=255)                                                                                                             cGCAAGTAGAAGCATACTGTGAGCTGGCTCAGAACCTAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGAATACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGAt                                         <  3:234683/149‑1 (MQ=255)                                                                                                              gCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATg                                        <  2:20586/149‑1 (MQ=255)                                                                                                                   tAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGt                                   >  5:229869/1‑149 (MQ=255)                                                                                                                         cATACGGTGAGCGGGCTCAGTACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTc                             <  2:272596/149‑1 (MQ=255)                                                                                                                          aTACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTcc                            <  1:98167/149‑1 (MQ=255)                                                                                                                          aTACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTcc                            <  1:236148/149‑1 (MQ=255)                                                                                                                           tACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTccc                                             >  7:9632/1‑131 (MQ=255)                                                                                                                           tACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTccc                                             <  8:9632/131‑1 (MQ=255)                                                                                                                                                  aGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTg      >  5:89850/1‑147 (MQ=255)                                                                                                                                                    tCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGc  >  8:172361/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTGTTTGAGCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGC  >  NC_002947/4210607‑4210902

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