Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,411,853	C→T	7.1%	G28S (GGT→AGT)	PP_1232 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,411,853	0	C	T	7.1%	163.2 / 2.1	56	G28S (GGT→AGT)	PP_1232	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/29);  new base T (2/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

TTTCGGCCAGCTTGTACACGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTGACGA  >  NC_002947/1411706‑1411991                                                                                                                                                     |                                                                                                                                           tttCGGCCAGCTTGTACACGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑ccc                                                                                                                                           <  6:234411/149‑1 (MQ=255)   tCGGCCAGCTTGTACACGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCa                                                                                                                                          <  7:18147/148‑1 (MQ=255)   tCGGCCAGCTTGTACACGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCa                                                                                                                                          <  5:181707/148‑1 (MQ=255)   tCGGCCAGCTTGTACACGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAg                                                                                                                                         >  8:169368/1‑149 (MQ=255)       ccAGCTTGTACACGCTGGTTTCGAGATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGTCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGAt                                                                                                                                     >  8:88403/1‑149 (MQ=255)          gCTTGTACACGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CTGGGTGATgcc                                                                                                                                  <  3:298945/149‑3 (MQ=255)          gCTTGTACACGCCGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACATGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCAAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATgg                                                                                                                                   >  5:231351/1‑148 (MQ=255)           cTTGTACACGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATgg                                                                                                                                   >  8:260217/1‑147 (MQ=255)            tggTACACGCTGGCTTCGACATAGTCCTTGAGCTGGGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAg                                                                                                                                <  8:89393/147‑1 (MQ=255)              gTACACGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGCCATCA‑CTGGGTGATGCCAgt                                                                                                                               <  2:171315/148‑2 (MQ=255)                 cacGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGgtcgt                                                                                                                           <  7:19537/149‑1 (MQ=255)                   cGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCa                                                                                                                         <  8:248219/149‑1 (MQ=255)                    gCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAg                                                                                                                        >  8:136228/1‑149 (MQ=255)                        gTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCa                                                                                                                      <  2:59008/147‑1 (MQ=255)                                   gTCCTTGAGCTGGGGGGCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAgg                                                                                                          <  6:379244/148‑1 (MQ=255)                                       ttGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCAt                                                                                                     >  3:266009/1‑149 (MQ=255)                                           gCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATcagg                                                                                                 >  8:184108/1‑149 (MQ=255)                                              gtgGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATcaggcag                                                                                              >  8:268201/1‑149 (MQ=255)                                                  gTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCa                                                                                          <  7:260217/149‑1 (MQ=255)                                                    cATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGc                                                                                        <  5:352029/149‑1 (MQ=255)                                                        gAGCTGCTTGTCCGGGCCGCGCTGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGTACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTca                                                                                    <  8:6175/149‑1 (MQ=255)                                                                    cTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAgg                                                                          >  6:141045/1‑147 (MQ=255)                                                                    cTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAgg                                                                          <  5:141045/147‑1 (MQ=255)                                                                         cgcgCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCACGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAg                                                                    >  8:158659/1‑148 (MQ=255)                                                                              agcagGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGcc                                                                                                       >  1:227562/1‑108 (MQ=255)                                                                              agcagGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGcc                                                                                                       <  2:227562/108‑1 (MQ=255)                                                                               gcagGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTc                                                             <  7:158659/149‑1 (MQ=255)                                                                                 agGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTcat                                                           <  1:250898/149‑1 (MQ=255)                                                                                   gCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTcatg                                                          >  8:133470/1‑148 (MQ=255)                                                                                   gCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTcatg                                                          >  3:320832/1‑148 (MQ=255)                                                                                        gCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGaa                                                    <  7:88403/149‑1 (MQ=255)                                                                                              cgcGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTc                                                <  8:115049/145‑1 (MQ=255)                                                                                                cGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGt                                            <  7:169368/149‑1 (MQ=255)                                                                                                   cccAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTc                                           <  6:231351/147‑1 (MQ=255)                                                                                                        ttGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCACCCCAGTGATGGCAGGACGTCAGCATGGCCgg                                                                                                               <  1:238352/75‑1 (MQ=255)                                                                                                        ttGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGACGTCAGCATGGCCgg                                                                                                               >  2:238352/1‑74 (MQ=255)                                                                                                        ttGGAGGTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGCCAGCTCGTCAGCAAGGCAGGGAAGGGCCATCAGGCAGCCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCa                                     >  8:107120/1‑148 (MQ=255)                                                                                                          ggTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCaaa                                   <  7:136228/148‑1 (MQ=255)                                                                                                          ggTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCaaa                                   <  7:268201/148‑1 (MQ=255)                                                                                                           gTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCa                                     >  4:140696/1‑145 (MQ=255)                                                                                                              cTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATAc                              <  4:320832/149‑1 (MQ=255)                                                                                                              cTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATAc                              <  4:266009/149‑1 (MQ=255)                                                                                                                ccTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTcatgca                                                        >  1:45008/1‑121 (MQ=255)                                                                                                                ccTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTcatgca                                                        <  2:45008/121‑1 (MQ=255)                                                                                                                 cTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATAcc                             <  7:184108/147‑1 (MQ=255)                                                                                                                    ttGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTg                        >  1:232030/1‑149 (MQ=255)                                                                                                                        ggcgAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGc                    >  6:68590/1‑149 (MQ=255)                                                                                                                           gagaCGATCGCAGAACTGGCATCA‑CCCAGTGATTGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCggg                 <  7:151537/149‑1 (MQ=255)                                                                                                                             gaCGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGgcg                                                                       >  4:62461/1‑93 (MQ=255)                                                                                                                             gaCGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGgcg                                                                       <  3:62461/93‑1 (MQ=255)                                                                                                                                      aGAACTGGCATCA‑CCCAGTGATGGCAGTTCGTCAGCATGGACGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGAATGAAGCTCTCGTCGGATCAAATATAACTAATGTAGCGGGCTCACACGtg       >  7:130547/1‑147 (MQ=255)                                                                                                                                      aGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGtt       >  3:150413/1‑148 (MQ=255)                                                                                                                                       gAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTga     <  4:150413/149‑1 (MQ=255)                                                                                                                                          cTGCCATCA‑CTGGGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTgacga  >  6:244871/1‑149 (MQ=255)                                                                                                                                          cTGCCATCA‑CTGGGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTgacga  >  2:152212/1‑149 (MQ=255)                                                                                                                                                  a‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGa                                        <  3:95582/103‑1 (MQ=255)                                                                                                                                                  a‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGa                                        >  4:95582/1‑103 (MQ=255)                                                                                                                                                     |                                                                                                                                           TTTCGGCCAGCTTGTACACGCTGGTTTCGACATAGTCCTTGAGCTGTGGGTCATTGAGCTGATTGACCTGGCCGCGCAGCAGGCTGAGCCAGGCACGGCCCAGTTGGTGCTCCTGTTGCGGCGAGACGATCGCAGAACTGGCATCA‑CCCAGTGATGGCAGGTCGTCAGCATGGCCGGGAAGGGCCATCAGGCAGGCCAGCGTCAGCAGGGTAGGGCGCAGTAGATTCATGCATGAAGCTCGCGTCGGATCAAAGATACCTACTGTAGCGGGCTCACACGTTGACGA  >  NC_002947/1411706‑1411991

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