Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,661,637	G→A	11.3%	S103S (TCG→TCA)	PP_1456 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,661,637	0	G	A	11.3%	94.9 / 2.8	39	S103S (TCG→TCA)	PP_1456	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (19/14);  new base A (2/2);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

TGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGGTGTAACCGCGGTGTGGCTGGCCATGGCATTTCGCGCCCGGCGGGCCGCCG  >  NC_002947/1661489‑1661785                                                                                                                                                     |                                                                                                                                                     tGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCg                                                                                                                                                      <  5:206943/149‑1 (MQ=255) tGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCg                                                                                                                                                      <  7:208434/149‑1 (MQ=255)      gtggCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCt                                                                                                                                                 <  7:27728/149‑1 (MQ=255)               gTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATccc                                                                                                                                        <  2:304448/149‑1 (MQ=255)                     tCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGCCCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGa                                                                                                                                  >  8:116180/1‑149 (MQ=255)                     tCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGa                                                                                                                                  <  2:25126/149‑1 (MQ=255)                       cTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGAt                                                                                                                                 <  8:59336/148‑1 (MQ=255)                        tATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATAc                                                                                                                               >  7:50879/1‑149 (MQ=255)                         aTGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATAcc                                                                                                                              >  1:32540/1‑149 (MQ=255)                            tCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTc                                                                                                                           <  8:50879/149‑1 (MQ=255)                                      gCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCg                                                                                                                 <  2:32540/149‑1 (MQ=255)                                         ggCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCa                                                                                                              >  2:160108/1‑149 (MQ=255)                                         ggCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCa                                                                                                              >  6:64740/1‑149 (MQ=255)                                          gCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGc                                                                                                               >  1:279394/1‑147 (MQ=255)                                              tCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATc                                                                                                         <  5:64740/149‑1 (MQ=255)                                                       ccAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCgg                                                                                                                                                     >  2:71829/1‑96 (MQ=255)                                                       ccAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCgg                                                                                                                                                     <  1:71829/96‑1 (MQ=255)                                                        cAGCGCCACTGCGTTCTTCATTTCCGACCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATACCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGa                                                                                               >  1:14741/1‑149 (MQ=255)                                                          gcgcCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCAGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGAcc                                                                                             >  8:164592/1‑149 (MQ=255)                                                          gcgcCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCAGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGAcc                                                                                             >  8:160113/1‑149 (MQ=255)                                                              cACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCgg                                                                                         >  6:313516/1‑149 (MQ=255)                                                                     ttcttcATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCTGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGTCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTTAt                                                                                  <  2:279394/149‑1 (MQ=255)                                                                           aTTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCAGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCa                                                                             <  7:164592/148‑1 (MQ=255)                                                                           aTTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCAGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCa                                                                             <  7:160113/148‑1 (MQ=255)                                                                                  aGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATg                                                                     <  3:108382/149‑1 (MQ=255)                                                                                  aGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATg                                                                     >  5:326260/1‑149 (MQ=255)                                                                                   gCTGATCGACTGGCTGGTGTTCAGCATCAGCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGGATCCCGGTGGATAGCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATgg                                                                    <  5:234159/149‑1 (MQ=255)                                                                                       aTCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTc                                                                <  1:287031/149‑1 (MQ=255)                                                                                       aTCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTc                                                                >  5:142203/1‑149 (MQ=255)                                                                                           aCTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGaaa                                                            >  5:334514/1‑149 (MQ=255)                                                                                                ctggTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGATGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGc                                                       <  6:326260/149‑1 (MQ=255)                                                                                                 tggTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCTGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGAt                                                                      <  1:165746/133‑1 (MQ=255)                                                                                                 tggTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCTGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGAt                                                                      >  2:165746/1‑133 (MQ=255)                                                                                                   gtgtTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCg                                                     <  6:334514/148‑1 (MQ=255)                                                                                                      ttCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGgtg                                                  >  1:215039/1‑148 (MQ=255)                                                                                                       tCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGgtgt                                                 >  2:34345/1‑148 (MQ=255)                                                                                                          gCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCCCCGGTGTAAcc                                             >  8:11764/1‑149 (MQ=255)                                                                                                                  cGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGGTGTAACCGCGGtgtg                                     >  4:162322/1‑149 (MQ=255)                                                                                                                                    cGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGGTGTAACCGCGGTGTGGCTGGCCAt                            >  5:242789/1‑140 (MQ=255)                                                                                                                                    cGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGGTGTAACCGCGGTGTGGCTGGCCAt                            <  6:242789/140‑1 (MQ=255)                                                                                                                                    cGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGGTGTAACCGCGGTGTGGCTGGCCATGGCATTTcg                   >  4:47985/1‑149 (MQ=255)                                                                                                                                       cTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGGTGTAACCGCGGTGTGGCTGGCCATGGCATTTCgcgc                >  4:190601/1‑149 (MQ=255)                                                                                                                                                    cGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGGTGTAACCGCGGTGTGGCTGGCCATGGCATTTCGCGCCCGGCGGgccgcc   >  7:299078/1‑149 (MQ=255)                                                                                                                                                     ggCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGGTGTAACCGCGGTGTGGCTGGCCATGGCATTTCGCGCCCGGCGGGccgccg  <  3:162322/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGTGGCTCAGTTCGGCGCTGAGCATCCCGGTGGATACCTTCATTTTCTTCGGCATGATCGGCGCGCTGACCCCGGCAGTGATCGGCACCGCGATGGCCTCGAAATTCGCCGGTGTAACCGCGGTGTGGCTGGCCATGGCATTTCGCGCCCGGCGGGCCGCCG  >  NC_002947/1661489‑1661785

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