Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,270,077	C→T	7.7%	D155D (GAC→GAT)	PP_1110 →	serine acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,270,077	0	C	T	7.7%	147.5 / 2.2	52	D155D (GAC→GAT)	PP_1110	serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (21/27);  new base T (2/2);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

TCAATGCCTGGCCGCCATGATCTCCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAGTGTGCTGGGCGCAATTCATGTCGGTGACGGTGC  >  NC_002947/1269936‑1270223                                                                                                                                               |                                                                                                                                                   tCAATGCCTGGCCGCCATGATCTCCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGAt                                                                                                                                                 >  1:9294/1‑145 (MQ=255) tCAATGCCTGGCCGCCATGATCTCCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGAt                                                                                                                                                 <  2:9294/145‑1 (MQ=255) tCAATGCCTGGCCGCCATGATCTCCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGc                                                                                                                                              >  3:240480/1‑148 (MQ=255)        cTGGCCGCCATGATCTCCAGGCGCCGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGACATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGc                                                                                                                                       >  7:200718/1‑149 (MQ=255)        cTGGCCGCCATGATCTCCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAATGTCGACCATTGGTGACGATTGCTACGTGCt                                                                                                                                      >  4:64392/1‑149 (MQ=255)         tGGCCGCCATGATCTCCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCtt                                                                                                                                     >  8:101302/1‑149 (MQ=255)           gccgccATGATCTCCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTg                                                                                                                                    <  3:151465/148‑1 (MQ=255)               ccATGATCTCCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTggcggc                                                                                                                               <  2:109429/149‑1 (MQ=255)                    aTCTCCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAAc                                                                                                                          >  2:48368/1‑149 (MQ=255)                      ctcCACCCG‑CGGAAAACTGCTTTCGGGGGCGGAGCTCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAAccc                                                                                                                        <  4:258242/149‑1 (MQ=255)                       tcCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCt                                                                                                                       >  6:141333/1‑149 (MQ=255)                         cAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTgg                                                                                                                     <  2:142554/149‑1 (MQ=255)                                cgGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGccc                                                                                                               <  5:268093/149‑1 (MQ=255)                                   aaaaaTGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCcgcg                                                                                                            <  4:240480/149‑1 (MQ=255)                                    aaaaTGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCgg                                                                                                           <  1:285059/149‑1 (MQ=255)                                      aaTGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATAACGGCATGGGCACCGTCATTGGTGAAACGTCGACAATTGGTGACGATTGCTACGTGGTTGGCGGCGTAACCCTGGGCGCCCGCGGTa                                                                                                         >  8:14913/1‑149 (MQ=255)                                        tGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTAt                                                                                                        >  6:22566/1‑148 (MQ=255)                                        tGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATc                                                                                                       <  5:144375/149‑1 (MQ=255)                                          cTTTCGGGGGCGGAGATCCATTTCCGATGCCGCTTTGGTGCGCGCTTCATCATCGTTCACGGCATGGGCACCGTCATTTGTTAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAg                                                                                                     <  2:325391/149‑1 (MQ=255)                                             tCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGt                                                                                                    <  3:32478/147‑1 (MQ=255)                                             tCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGACATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGt                                                                                                    <  8:200718/147‑1 (MQ=255)                                                 gggCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGATGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAAc                                                                                              >  2:242699/1‑149 (MQ=255)                                                 gggCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGATGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAAc                                                                                              >  4:157951/1‑149 (MQ=255)                                                     ggAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCt                                                                                          >  4:121781/1‑149 (MQ=255)                                                       agaTCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCg                                                                                        >  3:294196/1‑149 (MQ=255)                                                                tCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAATGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAg                                                                                <  3:64392/148‑1 (MQ=255)                                                                     aGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATcc                                                                          <  4:294196/149‑1 (MQ=255)                                                                       ccGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGa                                                                        >  2:316618/1‑149 (MQ=255)                                                                          cATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAG‑CCCAGGCATCCGACGtt                                                                    <  8:8676/149‑1 (MQ=255)                                                                           aTTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGATGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGtt                                                                    <  1:242699/149‑1 (MQ=255)                                                                           aTTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGATGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGtt                                                                    <  3:157951/149‑1 (MQ=255)                                                                           aTTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGtt                                                                    <  3:226527/149‑1 (MQ=255)                                                                           aTTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGtt                                                                    <  5:141333/149‑1 (MQ=255)                                                                              ggTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGGGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTgg                                                                  <  2:261622/148‑1 (MQ=255)                                                                                          atcatcGATCAGGGCTTGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGGCCGCGGTTTCATTGACTACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCa                                                     <  7:14913/149‑1 (MQ=255)                                                                                             atcGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTgg                                                                  >  5:125377/1‑133 (MQ=255)                                                                                             atcGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTgg                                                                  <  6:125377/133‑1 (MQ=255)                                                                                                 aTCACGGCATGGGCATCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGgc                                              <  4:235446/149‑1 (MQ=255)                                                                                                 aTCACGGCATGGGCATCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGgc                                              <  4:248295/149‑1 (MQ=255)                                                                                                           gggCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGc                                      >  2:39482/1‑147 (MQ=255)                                                                                                            ggCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCa                                    <  5:22566/148‑1 (MQ=255)                                                                                                               aCCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAgtgt                                >  4:127817/1‑149 (MQ=255)                                                                                                                     aTTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGGGCTTTCGCCAGTGTGCTggg                          >  8:196962/1‑149 (MQ=255)                                                                                                                     aTTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGGCAGTGTGCTggg                          >  4:9707/1‑149 (MQ=255)                                                                                                                          tGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAGTGTGCTGGGCGCaa                     >  1:15582/1‑149 (MQ=255)                                                                                                                             aaCGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAGTGTGCTGGGCGCAAtt                   >  2:35495/1‑148 (MQ=255)                                                                                                                                 tCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAGTGTGCTGGGCGCAATTCATGt              <  1:48368/149‑1 (MQ=255)                                                                                                                                       aTTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAGTGTGCTGGGCGCAATTCATGTCGGTGa        <  3:127817/149‑1 (MQ=255)                                                                                                                                        ttGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAGTGTGCTGGGCGCAATTCATGTCGGTGAc       >  2:122955/1‑149 (MQ=255)                                                                                                                                           gTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAGTGTGCTGGGCGCAATTCATGTCGGTGACGGt    <  7:101302/149‑1 (MQ=255)                                                                                                                                           gTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAGTGTGCTGGGCGCAATTCATGTCGGTGACGGt    >  7:193668/1‑149 (MQ=255)                                                                                                                                              aCGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAGTGTGCTGGGCGCAATTCATGTCGGTGACGGTGc  <  1:35495/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   TCAATGCCTGGCCGCCATGATCTCCAGGCG‑CGGAAAAATGCTTTCGGGGGCGGAGATCCATTTCCGCAGCCGCATTGGTGCGCGCTTCATCATCGATCACGGCATGGGCACCGTCATTGGTGAAACGTCGACCATTGGTGACGATTGCTACGTGCTTGGCGGCGTAACCCTGGGCGCCCGCGGTATCAGTGACAACCCCTCGAGCCCCAGGCATCCGACGTTGGGTAACCGGGTGCAGATCGGCGCTTTCGCCAGTGTGCTGGGCGCAATTCATGTCGGTGACGGTGC  >  NC_002947/1269936‑1270223

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