Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	280,200	A→C	8.0%	L17R (CTG→CGG)	sctS ←	sulfur compound ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	280,200	0	A	C	8.0%	132.1 / 5.1	50	L17R (CTG→CGG)	sctS	sulfur compound ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (20/26);  new base C (2/2);  total (22/28)
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.56e-01

GCGTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCG  >  NC_002947/280051‑280349                                                                                                                                                      |                                                                                                                                                       gcgTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTca                                                                                                                                                        <  1:1168911/150‑1 (MQ=255)         cgGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGgaag                                                                                                                                               >  2:935074‑M2/1‑146 (MQ=255)              gaaCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGgaagggcgc                                                                                                                                          >  2:728189‑M2/1‑141 (MQ=255)               aaCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGgaagggcgcg                                                                                                                                         >  2:130442‑M2/1‑140 (MQ=255)               aaCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGgaagggcgcg                                                                                                                                         <  1:422022‑M2/150‑11 (MQ=255)                aCGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGGGCAGCACCGTATAACCCGCGCCCTTCA‑GCAAgaaaggcgcgg                                                                                                                                        <  2:74924‑M2/150‑12 (MQ=255)                 cGACACGTAGACACGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGa                                                                                                                                       <  1:263489/150‑1 (MQ=255)                   acacGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGt                                                                                                                                     <  1:568717/150‑1 (MQ=255)                            cccGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACgc                                                                                                                            <  1:53744/150‑2 (MQ=255)                             ccGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGg                                                                                                                                                  <  2:97984‑M2/127‑2 (MQ=255)                             ccGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGg                                                                                                                                                  >  1:97984‑M2/1‑126 (MQ=255)                             ccGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAg                                                                                                                           <  2:1098581/150‑1 (MQ=255)                                   cagccagTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAg                                                                                                                     >  2:1145660/1‑150 (MQ=255)                                                            gACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACAt                                                                                            <  1:895841/150‑1 (MQ=255)                                                            gACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACAt                                                                                            >  2:271717/1‑150 (MQ=255)                                                             accgCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGGATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCGGCGATCATTTGGTAACATc                                                                                           >  1:316405/4‑150 (MQ=255)                                                                  cGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCAc                                                                                      <  1:469525/150‑1 (MQ=255)                                                                    cATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTc                                                                                                               <  2:116951/123‑1 (MQ=255)                                                                    cATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTc                                                                                                               >  1:116951/1‑123 (MQ=255)                                                                     aTCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGGCACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGa                                                                                   >  2:944322/1‑150 (MQ=255)                                                                                 gCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCgg                                                                       >  1:1229988/1‑150 (MQ=255)                                                                                    aaGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCCTGCTGaagggcgcgggttatacg                                                                                                                                 <  2:384431‑M2/90‑19 (MQ=255)                                                                                    aaGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCCTGCTGaagggcgcgggttatacg                                                                                                                                 >  1:384431‑M2/1‑72 (MQ=255)                                                                                      gCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGc                                                                  <  2:542500/150‑1 (MQ=255)                                                                                           gCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGc                                                             <  1:1145660/150‑1 (MQ=255)                                                                                               cAGACCGAAGAACAGGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCCATCATGTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGc                                                         <  1:474956/150‑1 (MQ=255)                                                                                                aGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCCTGCTGaagggcgcgggttatacggtgctgctcagt                                                                                                                     <  2:549910‑M2/90‑31 (MQ=255)                                                                                                aGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCCTGCTGaagggcgcgggttatacggtgctgctcagt                                                                                                                     >  1:549910‑M2/1‑60 (MQ=255)                                                                                                      aagaaCATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCAt                                                  >  1:586211/1‑150 (MQ=255)                                                                                                             tGCCGCCGGCACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCttt                                           <  1:944322/150‑1 (MQ=255)                                                                                                                     acacTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGt                                   >  2:385786/1‑150 (MQ=255)                                                                                                                         tGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTg                                                      <  1:586459/127‑1 (MQ=255)                                                                                                                         tGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTg                                                      >  2:586459/1‑127 (MQ=255)                                                                                                                         tGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGAt                               >  1:759643/1‑150 (MQ=255)                                                                                                                          gAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATg                              <  1:74924/150‑1 (MQ=255)                                                                                                                             cagcaCCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGcc                           <  2:422022/150‑1 (MQ=255)                                                                                                                               gcaCCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCtt                         >  1:436254/1‑150 (MQ=255)                                                                                                                                caCCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCttg                        >  2:1093770/1‑150 (MQ=255)                                                                                                                                caCCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCttg                        >  1:923747/1‑150 (MQ=255)                                                                                                                                       taACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATgg                 <  2:1229988/150‑1 (MQ=255)                                                                                                                                       taACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATgg                 <  1:714130/150‑1 (MQ=255)                                                                                                                                       taACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATgg                 <  1:1244228/150‑1 (MQ=255)                                                                                                                                         aCCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCt               <  1:271717/150‑1 (MQ=255)                                                                                                                                              cgcCCTTCA‑GCAGgaagggcgcgggttatacggtgctgctcagtgtcggcggcatgttcttcggtctggtgctgggctttgccctggcgctgatgcggctgtcgaagattctgccgctgaactggctggcgcgggtctacgtgtcgttct          <  2:199686‑M2/150‑138 (MQ=255)                                                                                                                                               gcCCTTCA‑GCAGgaagggcgcgggttatacggtgctgctcagtgtcggcggcatgttcttcggtctggtgctgggctttgccctggcgctgatgcggctgtcgaagattctgccgc                                           >  2:960582‑M2/1‑12 (MQ=255)                                                                                                                                               gcCCTTCA‑GCAGgaagggcgcgggttatacggtgctgctcagtgtcggcggcatgttcttcggtctggtgctgggctttgccctggcgctgatgcggctgtcgaagattctgccgc                                           <  1:960582‑M2/116‑105 (MQ=255)                                                                                                                                                 ccTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGc                                                             <  2:703367/96‑1 (MQ=255)                                                                                                                                                 ccTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGc                                                             >  1:703367/1‑96 (MQ=255)                                                                                                                                                   ttCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGtt     >  2:504661/1‑150 (MQ=255)                                                                                                                                                    tCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTc    <  2:436254/150‑1 (MQ=255)                                                                                                                                                    tCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTc    <  1:1093770/150‑1 (MQ=255)                                                                                                                                                      a‑gcagGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCg  <  1:504661/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                       GCGTGCCGCGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCG  >  NC_002947/280051‑280349

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