Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 456,490 T→C 10.1% intergenic (+59/‑187) tig → / → clpP peptidyl‑prolyl cis/trans isomerase (trigger factor)/proteolytic subunit of ClpA‑ClpP and ClpX‑ClpP ATP‑dependent serine proteases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913456,4900TC10.1% 182.7 / 9.3 59intergenic (+59/‑187)tig/clpPpeptidyl‑prolyl cis/trans isomerase (trigger factor)/proteolytic subunit of ClpA‑ClpP and ClpX‑ClpP ATP‑dependent serine proteases
Reads supporting (aligned to +/- strand):  ref base T (26/27);  new base C (3/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AAAGTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACTAGCTGATAATCCGTCCATA  >  NC_000913/456378‑456626
                                                                                                                 |                                                                                                                                        
aaagTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGa                                                                                                                >  3:560017/1‑139 (MQ=255)
aaagTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGa                                                                                                                <  4:560017/139‑1 (MQ=255)
   gTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGt                                                                                                                                 <  3:1133877/119‑1 (MQ=255)
   gTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGt                                                                                                                                 >  4:1133877/1‑119 (MQ=255)
           aaaaGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCttttt                                                                                                                                      <  2:235673/106‑1 (MQ=255)
           aaaaGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCttttt                                                                                                                                      >  1:235673/1‑106 (MQ=255)
              ttcagtcACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg                                                                                         <  1:476534/142‑1 (MQ=255)
              ataaaCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTg                                                                                                                                  >  2:359621/3‑107 (MQ=255)
              aGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTg                                                                                                                                  <  1:359621/107‑1 (MQ=255)
              aGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGc                                                                                                                     <  3:52739/120‑1 (MQ=255)
              aGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGc                                                                                                                     >  4:52739/1‑120 (MQ=255)
              aGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg                                                                                         >  2:759506/1‑148 (MQ=255)
                aaaCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTc                                                                                           <  4:462748/144‑1 (MQ=255)
                aaaCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTc                                                                                           >  3:462748/1‑144 (MQ=255)
                  aCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGa                                                                                                                            <  4:507079/109‑1 (MQ=255)
                  aCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGa                                                                                                                            >  3:507079/1‑109 (MQ=255)
                  aCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGtt                                                                                     <  1:759506/148‑1 (MQ=255)
                                   gatgaACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTaaca                                                                            <  1:554150/140‑1 (MQ=255)
                                   gatgaACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTaaca                                                                            >  2:554150/1‑140 (MQ=255)
                                        aCCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg                                                                                         <  1:633182/122‑1 (MQ=255)
                                        aCCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg                                                                                         >  2:633182/1‑122 (MQ=255)
                                         ccAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCttttt                                                                                                                                      >  1:1106693/1‑76 (MQ=255)
                                         ccAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCttttt                                                                                                                                      <  2:1106693/76‑1 (MQ=255)
                                          cagcagGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg                                                                                         <  4:1128964/120‑1 (MQ=255)
                                          cagcagGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg                                                                                         >  3:1128964/1‑120 (MQ=255)
                                                   tAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGgtg                                                                                       <  4:636979/113‑1 (MQ=255)
                                                   tAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGgtg                                                                                       >  3:636979/1‑113 (MQ=255)
                                                   tAATTTACGCAGCATAACGCGCTAAATCCGCTCAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTc                                                                                                                                <  3:740738/72‑1 (MQ=255)
                                                   tAATTTACGCAGCATAACGCGCTAAATCCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTc                                                                                                                                >  4:740738/1‑72 (MQ=255)
                                                     aTTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGc                                                                                                                     <  3:362777/81‑1 (MQ=255)
                                                     aTTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGc                                                                                                                     >  4:362777/1‑81 (MQ=255)
                                                            cagcaTAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGAt                                            >  1:480397/1‑148 (MQ=255)
                                                            cagcaTAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTc                                                                                                                                >  2:974921/1‑63 (MQ=255)
                                                            cagcaTAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTc                                                                                                                                <  1:974921/63‑1 (MQ=255)
                                                                    cgcgCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAg                                                                                                   >  4:426176/1‑84 (MQ=255)
                                                                    cgcgCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAg                                                                                                   <  3:426176/84‑1 (MQ=255)
                                                                          aaaTTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGaa                                                      >  4:350470/1‑123 (MQ=255)
                                                                          aaaTTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGaa                                                      <  3:350470/123‑1 (MQ=255)
                                                                              tCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCAGTGTTAGCGTAACAACAAAAGGTTGTTATGCTTGAAata                                                   <  1:328625/122‑1 (MQ=255)
                                                                              tCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATCTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAata                                                   >  2:328625/1‑122 (MQ=255)
                                                                                 cacaAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAg                          >  2:581943/1‑144 (MQ=255)
                                                                                 cacaAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAg                          <  1:581943/144‑1 (MQ=255)
                                                                                   caAAGGCCCGTCACCGCCAGGTGGTGGGCTtttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttc                                                            <  3:924722‑M2/108‑79 (MQ=255)
                                                                                   caAAGGCCCGTCACCGCCAGGTGGTGGGCTtttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttc                                                            >  4:924722‑M2/1‑30 (MQ=255)
                                                                                   caAAGGCCCGTCACCGCCAGGTGGTGGGCTtttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggtt                                                                                          >  3:769524‑M2/1‑30 (MQ=255)
                                                                                   caAAGGCCCGTCACCGCCAGGTGGTGGGCTtttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggtt                                                                                          <  4:769524‑M2/78‑49 (MQ=255)
                                                                                        gCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACTAGCTg                <  2:480397/148‑1 (MQ=255)
                                                                                         cccGTCACCGCCAGGTGGTGGGCCtttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttcagtcactttcgctttcgcc                                         >  2:82436‑M2/1‑24 (MQ=255)
                                                                                         cccGTCACCGCCAGGTGGTGGGCCtttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttcagtcactttcgctttcgcc                                         <  1:82436‑M2/121‑98 (MQ=255)
                                                                                          ccGTCACCGCCAGGTGGTGGGCCtttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttcagtcactttcgctttcgccagtacagcttcaacagcctgttcttcca             <  2:218670‑M2/148‑126 (MQ=255)
                                                                                           cGTCACCGCCAGGTGGTGGGCCtttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttct                                                                 <  2:5113‑M2/95‑74 (MQ=255)
                                                                                           cGTCACCGCCAGGTGGTGGGCCtttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttct                                                                 >  1:5113‑M2/1‑22 (MQ=255)
                                                                                                ccGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTc                                                                                           <  2:879819/64‑1 (MQ=255)
                                                                                                ccGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTc                                                                                           >  1:879819/1‑64 (MQ=255)
                                                                                                 cGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACt                     <  2:667297/133‑1 (MQ=255)
                                                                                                 cGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACt                     >  3:56768/1‑133 (MQ=255)
                                                                                                 cGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACt                     >  1:667297/1‑133 (MQ=255)
                                                                                                 cGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACt                     <  4:56768/133‑1 (MQ=255)
                                                                                                     aGGTGGTGGGCCtttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttcagtcactttcgctttcgccagtacagcttcaacagcctgttcttccagagcaacattg  >  1:218670‑M2/1‑12 (MQ=255)
                                                                                                      ggtggtGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAggg                        >  1:186930/1‑125 (MQ=255)
                                                                                                      ggtggtGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAggg                        <  2:186930/125‑1 (MQ=255)
                                                                                                               cTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCg                                                                                 >  1:186983/1‑59 (MQ=255)
                                                                                                               cTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCg                                                                                 <  2:186983/59‑1 (MQ=255)
                                                                                                                 |                                                                                                                                        
AAAGTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACTAGCTGATAATCCGTCCATA  >  NC_000913/456378‑456626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: