Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,422,516	C→T	7.4%	A366V (GCG→GTG)	PP_1246 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,422,516	0	C	T	7.4%	155.1 / 2.2	54	A366V (GCG→GTG)	PP_1246	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

CCCGGCCGGCAGATGAAGGTGTACCTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTGGGGCAACGACGCGCAGAAAATCACCCGCAACCGCATCGGCA  >  NC_002947/1422374‑1422659                                                                                                                                               |                                                                                                                                                 cccGGCCGGCAGATGAAGGTGTACCTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑gg                                                                                                                                            >  3:174305/1‑148 (MQ=255)       cggcAGATGAAGGTGTACCTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCtt                                                                                                                                     <  7:118749/149‑1 (MQ=255)               gAAGGTGTACCTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGa                                                                                                                             >  1:332485/1‑149 (MQ=255)                aaGGTGTACCTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGAt                                                                                                                            <  1:5764/149‑1 (MQ=255)                   gtgtACCTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTc                                                                                                                                    <  5:229283/138‑1 (MQ=255)                   gtgtACCTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTc                                                                                                                                    >  6:229283/1‑138 (MQ=255)                      tACCTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGac                                                                                                                      >  4:259986/1‑149 (MQ=255)                      tACCTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGa                                                                                                                       <  2:332485/148‑1 (MQ=255)                        ccTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGacc                                                                                                                     <  8:291349/148‑1 (MQ=255)                         cTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGaccac                                                                                                                   >  8:279813/1‑149 (MQ=255)                          tGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGaccac                                                                                                                   >  1:87027/1‑148 (MQ=255)                              gTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGc                                                                                                              <  7:279813/149‑1 (MQ=255)                                 cACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCg                                                                                                            <  3:259986/148‑1 (MQ=255)                                  aCCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGac                                                                                                          >  1:263989/1‑149 (MQ=255)                                   cccGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGacc                                                                                                         >  5:177080/1‑149 (MQ=255)                                    ccGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGAc                                                                                                                                 <  6:156736/124‑1 (MQ=255)                                    ccGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGAc                                                                                                                                 >  5:156736/1‑124 (MQ=255)                                               aCGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGAtgtg                                                                                             >  3:297942/1‑149 (MQ=255)                                                cGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTgg                                                                                            <  5:227425/149‑1 (MQ=255)                                                       cTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCt                                                                                     <  2:263989/149‑1 (MQ=255)                                                          cGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGtt                                                                                  >  4:44006/1‑149 (MQ=255)                                                            cccGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCa                                                                                >  6:144884/1‑149 (MQ=255)                                                              cGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATc                                                                              <  6:177080/149‑1 (MQ=255)                                                              cGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATc                                                                              <  4:174305/149‑1 (MQ=255)                                                               gTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATc                                                                              >  8:61144/1‑148 (MQ=255)                                                                      agcGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGa                                                                      <  5:144884/149‑1 (MQ=255)                                                                           cccTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGa                                                                                                                       >  6:290939/1‑95 (MQ=255)                                                                           cccTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGa                                                                                                                       <  5:290939/95‑1 (MQ=255)                                                                             cTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCg                                                               >  2:294351/1‑149 (MQ=255)                                                                                actacGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGTGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGacca                                                                                                        >  6:322519/1‑105 (MQ=255)                                                                                actacGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGTGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGacca                                                                                                        >  6:316932/1‑105 (MQ=255)                                                                                actacGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGTGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGacca                                                                                                        <  5:322519/105‑1 (MQ=255)                                                                                actacGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGTGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGacca                                                                                                        <  5:316932/105‑1 (MQ=255)                                                                                  tacGTCGCGGTCACCCCGGACGACTCGTTAGCACTCGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTACGACCTGATCAAGACCACCTCGCCGACCACCACCGATCAGGCCAACCTGAACATCGACAAAGACATGTCGGCGAt                                                          >  2:232725/1‑149 (MQ=255)                                                                                     gTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGCCCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGAtctc                                                       >  5:356867/1‑149 (MQ=255)                                                                                      tCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCg                                                      <  4:297942/149‑1 (MQ=255)                                                                                                cccGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCTACGAAATATTCGACCTGATCAAGACCACATCGCCGACCCCCACCGATGTGGCAAACCTGTTAATCGTAAAAGACCTTTCGGCGAAATCGAACAACAta                                             >  1:271953/1‑148 (MQ=255)                                                                                                 ccGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTg                                           <  2:87027/149‑1 (MQ=255)                                                                                                  cGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCTATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTgg                                          <  1:232725/149‑1 (MQ=255)                                                                                                   ggACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTggg                                         <  7:61144/149‑1 (MQ=255)                                                                                                       gAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTgg                                                                                            <  3:317681/94‑1 (MQ=255)                                                                                                       gAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTgg                                                                                            >  4:317681/1‑94 (MQ=255)                                                                                                            gTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTgg                                          <  2:330812/139‑1 (MQ=255)                                                                                                            gTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTgg                                          >  1:330812/1‑139 (MQ=255)                                                                                                             tttGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCg                                                            <  6:25775/120‑1 (MQ=255)                                                                                                             tttGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCg                                                            >  5:25775/1‑120 (MQ=255)                                                                                                                  tctGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTc                                                                           >  3:101116/1‑100 (MQ=255)                                                                                                                  tctGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTc                                                                           <  4:101116/100‑1 (MQ=255)                                                                                                                  tctGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTGGGGCAACGACGCGCAg                           >  2:299720/1‑148 (MQ=255)                                                                                                                  tctGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTGGGGCAACGACGCGCAg                           >  2:217255/1‑148 (MQ=255)                                                                                                                           gACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTGGGGCAACGACGCGCAGAAAATCAccc                 <  1:217255/149‑1 (MQ=255)                                                                                                                              cAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTGGGGCAACGACGCGCAGAAAATCACCCGc               >  5:201213/1‑148 (MQ=255)                                                                                                                                  tACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTGGGGCAACGACGCGCAGAAAATCACCCGCAACCg          >  3:107291/1‑149 (MQ=255)                                                                                                                                          aGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTGGGGCAACGACGCGCAGAAATTCACCCGCAACCGCATCggga  <  6:356867/149‑3 (MQ=255)                                                                                                                                               |                                                                                                                                                 CCCGGCCGGCAGATGAAGGTGTACCTGACGTTCACCCGCTCCAGCTACGGCGTGCTCCGCCCGTTCAGCAGCGACCCTTACTACGTCGCGGTCACCCCGGACGAGTCGTTTGCTCTGGACGCGACCAAGTACAAGGAAGCTGCGCT‑GGAAATCTTCGACCTGATCAAGACCACCTCGCCGACCACCACCGATGTGGCCAACCTGTTCATCGTCAAAGACCTGTCGGCGATCTCGAACAACATCTGGGGCAACGACGCGCAGAAAATCACCCGCAACCGCATCGGCA  >  NC_002947/1422374‑1422659

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