Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,938,896	A→G	7.5%	L68S (TTG→TCG)	PP_1738 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,938,896	0	A	G	7.5%	154.7 / 2.7	53	L68S (TTG→TCG)	PP_1738	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base G (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

ATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTC  >  NC_002947/1938754‑1939023                                                                                                                                                |                                                                                                                                aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAg                                                                                                                          >  1:613396/1‑150 (MQ=255) aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATAAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAg                                                                                                                          >  1:615297/1‑150 (MQ=255)  ttGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTAGGGCAACGTGAACTTGCGCCCGGA‑GCGCGAGTTCac                                                                                                                          >  1:874951/1‑148 (MQ=255)  ttGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTAGGGCAACGTGAACTTGCGC‑CCGGAGCGCGAGTTCac                                                                                                                          <  2:874951/149‑2 (MQ=255)    gATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCt                                                                                                                       >  1:983513/1‑150 (MQ=255)      tgctgctgcTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGt                                                                                                                     <  1:1112470/150‑1 (MQ=255)      tgctgctgcTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGt                                                                                                                     >  2:1094535/1‑150 (MQ=255)        ctgctgctgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAg                                                                                                                   <  1:980403/150‑1 (MQ=255)           ctgctgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCg                                                                                                                >  1:1021316/1‑150 (MQ=255)             gctgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGa                                                                                                              <  1:212397/150‑1 (MQ=255)               tgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCACAGAAGTTCTCGCGACACGTTCATTCGCGCCACGC‑GGTCAATTGCAGAGGGTAGTCGGACt                                                                                                            >  2:140405/1‑150 (MQ=255)                    tcatcgcCGGCACCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGc                                                                                                                        <  2:498336/129‑1 (MQ=255)                    ggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTTGCGC‑CCGGAGCGCGAGTTCacg                                                                                                                         <  1:535276/132‑3 (MQ=255)                    ggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTTGCGC‑CCGGAGCGCGAGTTCacg                                                                                                                         >  2:535276/1‑130 (MQ=255)                         gcCGGTGCCGGCGATGAAGAAGATCAGAGACCCTTCGTTATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTc                                                                                                  <  2:1104719/150‑1 (MQ=255)                            ggTGCCGGCGATGATTAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCATAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTTTAGTCGGACTGAGTGATTTCAtc                                                                                               <  2:562803/150‑1 (MQ=255)                                         atgaAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGc                                                                                  >  1:723254/1‑150 (MQ=255)                                           gaAGATCAGCGGCGCTTCTTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTc                                                                                >  2:1022496/1‑150 (MQ=255)                                              gATCAGCGGCGCTTCGTGATCCTGGCGTGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAgg                                                                             <  2:965765/150‑1 (MQ=255)                                                         cTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcg                                                                  <  2:1021316/150‑1 (MQ=255)                                                               gATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTACGgcggcg                                                            <  2:983513/150‑1 (MQ=255)                                                                 tCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtg                                                          >  2:486659/1‑150 (MQ=255)                                                                    tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtc                                                                                               <  1:766135/110‑1 (MQ=255)                                                                    tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtc                                                                                               >  2:766135/1‑110 (MQ=255)                                                                    tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTc                                                                    <  2:1040129/137‑1 (MQ=255)                                                                    tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTc                                                                    >  1:1040129/1‑137 (MQ=255)                                                                    tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtggtg                                                       <  2:615297/150‑1 (MQ=255)                                                                    tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtggtg                                                       <  2:613396/150‑1 (MQ=255)                                                                    tGGCGGGCCAGGCGGTACTTGAGCTAGGTGACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGACTGTAGTCTGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtggtg                                                       >  1:174684/1‑150 (MQ=255)                                                                               gCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAt                                                  >  1:498551/1‑144 (MQ=255)                                                                               gCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAt                                                  <  2:498551/144‑1 (MQ=255)                                                                                      ttGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAAc                                     >  1:218381/1‑150 (MQ=255)                                                                                         aGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACggg                                  <  1:786080/150‑1 (MQ=255)                                                                                         aGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACggg                                  <  2:551956/150‑1 (MQ=255)                                                                                           cttcttCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGtt                                >  1:673280/1‑150 (MQ=255)                                                                                              cttcACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGt                             >  2:427231/1‑150 (MQ=255)                                                                                                  aCTGCCGAGAAGTTGTCGGGCAACGTGAACTCGCTCTGCGG‑GCGCAAGTTCATGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCTCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAaga                         <  2:178860/150‑1 (MQ=255)                                                                                                      tccAGGAGTTGGCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTAAGGCTATAGTCGCACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGt                     >  2:343313/2‑150 (MQ=255)                                                                                                       ccAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGc                                                                                     <  1:948034/85‑1 (MQ=255)                                                                                                       ccAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGc                                                                                     >  2:948034/1‑85 (MQ=255)                                                                                                       ccAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTa                    <  2:723254/150‑1 (MQ=255)                                                                                                                   gggCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGc        <  1:427231/150‑1 (MQ=255)                                                                                                                     gCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGCTTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAg      >  2:15236/1‑150 (MQ=255)                                                                                                                     gCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAg      >  2:335329/1‑150 (MQ=255)                                                                                                                       aaCGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAagag                        >  2:673120/1‑130 (MQ=255)                                                                                                                       aaCGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAagag                        <  1:673120/130‑1 (MQ=255)                                                                                                                         cGTGACCTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTc  <  1:1094535/150‑1 (MQ=255)                                                                                                                            gAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAtggtg                                              <  2:320999/103‑1 (MQ=255)                                                                                                                            gAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAtggtg                                              >  1:320999/1‑103 (MQ=255)                                                                                                                                tCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAagag                        <  2:567485/121‑1 (MQ=255)                                                                                                                                tCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAagag                        >  1:567485/1‑121 (MQ=255)                                                                                                                                      tCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTg                               >  1:161006/1‑108 (MQ=255)                                                                                                                                      tCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTg                               <  2:161006/108‑1 (MQ=255)                                                                                                                                                aaGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAaga                         >  1:1133444/1‑105 (MQ=255)                                                                                                                                                aaGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAaga                         <  2:1133444/105‑1 (MQ=255)                                                                                                                                                |                                                                                                                                ATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTC  >  NC_002947/1938754‑1939023

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑