Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,697,318	T→G	7.6%	intergenic (+28/‑49)	hslV → / → hslU	ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,697,318	0	T	G	7.6%	163.0 / 4.8	55	intergenic (+28/‑49)	hslV/hslU	ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base T (28/22);  new base G (2/2);  total (30/25)
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.99e-01

AACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACC  >  NC_002947/5697174‑5697463                                                                                                                                                  |                                                                                                                                                  aaCAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAg                                                                                                                                               >  5:161411/1‑149 (MQ=255) aaCAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAg                                                                                                                                               >  2:2806/1‑149 (MQ=255)   cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGATCGCCAACCACAACCTGACCATCGAGTAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCTCCCCTTCCCGCTGAAggg                                                                                                                                             >  8:77848/1‑149 (MQ=255)       accgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCtt                                                                                                                                         <  6:161411/149‑1 (MQ=255)            ccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCa                                                                                                                                    <  5:52907/149‑1 (MQ=255)              tctcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGACGGCGACATCTGCGTGTTAACCAAACACAAACTGACAATCGAGGAGCAGGAACTGGCCGAATGATTCATCTGTTCT‑GGCTCCCATTCCAGGTGAAGGGGCTTTTCCACg                                                                                                                                  >  4:109031/1‑149 (MQ=255)                   gCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGAc                                                                                                                             >  1:5739/1‑149 (MQ=255)                       gcgAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCAc                                                                                                                         >  2:124959/1‑149 (MQ=255)                         gAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCt                                                                                                                       >  1:44080/1‑149 (MQ=255)                           aaTCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GCCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGc                                                                                                                     >  8:141745/1‑149 (MQ=255)                               gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTCGGCGCCCCTTCACCGGGAAGGGGcgcca                                                                                                                                       <  1:16350/128‑5 (MQ=255)                               gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCACCGGGAAGGGGcgcca                                                                                                                                       >  2:16350/1‑123 (MQ=255)                                       cGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAggg                                                                                                                                             <  5:52618/113‑1 (MQ=255)                                       cGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAggg                                                                                                                                             >  6:52618/1‑113 (MQ=255)                                              aaCATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATc                                                                                                   <  7:72057/148‑1 (MQ=255)                                                 aTCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCg                                                                                                           <  7:6541/137‑1 (MQ=255)                                                 aTCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCg                                                                                                           >  8:6541/1‑137 (MQ=255)                                                 aTCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGt                                                                                               >  5:180904/1‑149 (MQ=255)                                                   cGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATc                                                                                                   >  5:160423/1‑143 (MQ=255)                                                   cGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATc                                                                                                   <  6:160423/143‑1 (MQ=255)                                                    gCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCa                                                                                            <  2:44080/149‑1 (MQ=255)                                                             aTCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCa                                                                                            >  2:152266/1‑140 (MQ=255)                                                             aTCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCa                                                                                            <  1:152266/140‑1 (MQ=255)                                                                  cGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCga                                                                               <  8:30051/148‑1 (MQ=255)                                                                          ccaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCa                                                                                                                  >  8:143776/1‑105 (MQ=255)                                                                                    ccTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAAc                                                            >  6:22193/1‑149 (MQ=255)                                                                                         ccATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCca                                                       <  8:89463/149‑1 (MQ=255)                                                                                            tCGAGGAGCAGGACCTGGCCGACTGATTCAGCTCTTCT‑GGCGCCCCTTCCCGGTGAAGGTGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcat                                                    <  7:77848/149‑1 (MQ=255)                                                                                               aggagCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc                                                                                 >  6:269593/1‑117 (MQ=255)                                                                                               aggagCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc                                                                                 <  5:269593/117‑1 (MQ=255)                                                                                                   gCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACt                                                                <  8:209306/130‑1 (MQ=255)                                                                                                   gCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGc                                             >  7:103287/1‑149 (MQ=255)                                                                                                    cAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGc                                                                                                                                 >  2:217943/1‑64 (MQ=255)                                                                                                    cAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGc                                                                                                                                 <  1:217943/64‑1 (MQ=255)                                                                                                    cAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCAt                                                                                           <  4:197386/102‑1 (MQ=255)                                                                                                    cAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCAt                                                                                           >  3:197386/1‑102 (MQ=255)                                                                                                      ggACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAg                                          <  6:180904/149‑1 (MQ=255)                                                                                                               cGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGc                                  >  8:157725/1‑148 (MQ=255)                                                                                                                     aTTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcgc                           >  3:184820/1‑149 (MQ=255)                                                                                                                      ttCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcgcg                          >  2:149969/1‑149 (MQ=255)                                                                                                                          gCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAgg                                         >  6:45990/1‑130 (MQ=255)                                                                                                                          gCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAgg                                         >  6:41911/1‑130 (MQ=255)                                                                                                                          gCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAgg                                         <  5:45990/130‑1 (MQ=255)                                                                                                                          gCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAgg                                         <  5:41911/130‑1 (MQ=255)                                                                                                                           cTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTg                     <  2:5739/149‑1 (MQ=255)                                                                                                                             gTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTgg                    <  1:242456/148‑1 (MQ=255)                                                                                                                             gTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATAATTAAAATGAACCCCCCCGAGAAAGTACAAGAACAAAAACGCAAAATAAACGGACAGAACAACACCAAACGCGCCGTgg                    >  2:242456/1‑148 (MQ=255)                                                                                                                              ttCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTgg                    <  1:124959/147‑1 (MQ=255)                                                                                                                                    gcgcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTg              <  7:157725/148‑1 (MQ=255)                                                                                                                                      gcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGcc            >  1:101820/1‑148 (MQ=255)                                                                                                                                      gcCCCTTCCCGGCGAAGGGGCTGTTCCACGCTGACTCACCTGCCCAACGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGcc            <  2:101820/148‑1 (MQ=255)                                                                                                                                      gcCCCTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTg              <  8:102548/146‑1 (MQ=255)                                                                                                                                      gcCCCTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTg              >  7:102548/1‑146 (MQ=255)                                                                                                                                       ccccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCt          >  2:83914/1‑149 (MQ=255)                                                                                                                                               cGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAAcc  >  3:216490/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                  AACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCAACC  >  NC_002947/5697174‑5697463

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