Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,697,314 C→A 7.3% intergenic (+24/‑53) 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_0029475,697,3140CA 7.3% 172.1 / 3.9 55intergenic (+24/‑53)hslV/hslUATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base A (2/2);  total (29/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 = 2.21e-01

CTGCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCT  >  NC_002947/5697168‑5697455
                                                                                                                                                   |                                                                                                                                             
ctgctgAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCg                                                                                                                                             <  7:36993/149‑1 (MQ=255)
      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                                                                                                   >  8:6541/1‑137 (MQ=255)
                                                      aTCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCg                                                                                                   <  7:6541/137‑1 (MQ=255)
                                                      aTCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGt                                                                                       >  5:180904/1‑149 (MQ=255)
                                                        cGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATc                                                                                           <  6:160423/143‑1 (MQ=255)
                                                        cGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATc                                                                                           >  5:160423/1‑143 (MQ=255)
                                                         gCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCa                                                                                    <  2:44080/149‑1 (MQ=255)
                                                                  aTCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCa                                                                                    <  1:152266/140‑1 (MQ=255)
                                                                  aTCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCa                                                                                    >  2:152266/1‑140 (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                                                                         <  5:269593/117‑1 (MQ=255)
                                                                                                    aggagCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc                                                                         >  6:269593/1‑117 (MQ=255)
                                                                                                        gCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACt                                                        <  8:209306/130‑1 (MQ=255)
                                                                                                        gCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGc                                     >  7:103287/1‑149 (MQ=255)
                                                                                                         cAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGc                                                                                                                         <  1:217943/64‑1 (MQ=255)
                                                                                                         cAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGc                                                                                                                         >  2:217943/1‑64 (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                                 <  5:41911/130‑1 (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)
                                                                                                                                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      >  7:102548/1‑146 (MQ=255)
                                                                                                                                           gcCCCTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTg      <  8:102548/146‑1 (MQ=255)
                                                                                                                                            ccccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCt  >  2:83914/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CTGCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCT‑GGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCT  >  NC_002947/5697168‑5697455

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