Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	280,203	A→T	8.0%	L16Q (CTG→CAG)	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,203	0	A	T	8.0%	148.9 / 5.3	50	L16Q (CTG→CAG)	sctS	sulfur compound ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (20/26);  new base T (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 = 7.06e-01

CGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCGGC  >  NC_002947/280059‑280351                                                                                                                                                  |                                                                                                                                                     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                                                        >  2:586459/1‑127 (MQ=255)                                                                                                                 tGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTg                                                        <  1:586459/127‑1 (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:1244228/150‑1 (MQ=255)                                                                                                                               taACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATgg                   <  1:714130/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                                                               >  1:703367/1‑96 (MQ=255)                                                                                                                                         ccTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGc                                                               <  2:703367/96‑1 (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)                                                                                                                                                gcagGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCgg   >  1:562502/1‑150 (MQ=255)                                                                                                                                                 cagGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCCACAGTTCCGGc  >  1:696414/1‑150 (MQ=255)                                                                                                                                                 cagGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCGGc  >  2:591629/1‑150 (MQ=255)                                                                                                                                                  agGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCgg   <  1:448573/148‑1 (MQ=255)                                                                                                                                                  agGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCgg   >  2:448573/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CGGAAGAACGACACGTAGACCCGCGCCAGCCAGTTCAGCGGCAGAATCTTCGACAGCCGCATCAGCGCCAGGGCAAAGCCCAGCACCAGACCGAAGAACATGCCGCCGACACTGAGCAGCACCGTATAACCCGCGCCCTTCA‑GCAGGAAGGGCGCGGAGTCGACAACGAGTTGCAGGCTTTCAGCGATCATTTGGTGACATCGGCACCGAAGTACTTCTCGGACAGCTTGGCCAGCGTGCCATCAGCTTTGAGCTTGTCGATGGCCTTGTTGATGGCTGCCAACAGTTCCGGC  >  NC_002947/280059‑280351

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