Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,946,955 G→A A141V (GCC→GTC)  ruvC ← component of RuvABC resolvasome, endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,946,9550GA98.2% 184.8 / ‑4.2 54A141V (GCC→GTC) ruvCcomponent of RuvABC resolvasome, endonuclease
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (20/33);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

TAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGGCGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCT  >  NC_000913/1946856‑1947054
                                                                                                   |                                                                                                   
tAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATcgg                                                                                                     >  2:1981335/1‑100 (MQ=255)
  aCGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATcgacg                                                                                                   <  2:391581/100‑1 (MQ=255)
  aCGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATcgacg                                                                                                   <  1:1771244/100‑1 (MQ=255)
  aCGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATcgacg                                                                                                   >  1:1820970/1‑100 (MQ=255)
      aGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCAt                                                                                               <  2:351495/100‑1 (MQ=255)
        tCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATcc                                                                                             <  2:1311773/100‑1 (MQ=255)
          gCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGc                                                                                           >  1:1151876/1‑100 (MQ=255)
           cccTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGcc                                                                                          <  1:1304153/100‑1 (MQ=255)
           cccTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGcc                                                                                          <  1:158278/100‑1 (MQ=255)
            ccTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCt                                                                                         >  2:1110353/1‑100 (MQ=255)
             ctctCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCtg                                                                                        <  2:1306994/100‑1 (MQ=255)
               ctcGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCtgtg                                                                                      <  1:157767/100‑1 (MQ=255)
                      ggTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGc                                                                               >  1:459252/1‑100 (MQ=255)
                       gTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCg                                                                              >  2:676666/1‑100 (MQ=255)
                          cAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGc                                                                           <  1:4357/100‑1 (MQ=255)
                              cgcgATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGtt                                                                       >  2:1408307/1‑100 (MQ=255)
                                cgATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTc                                                                     <  1:269086/100‑1 (MQ=255)
                                cgATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTc                                                                     <  2:1965445/100‑1 (MQ=255)
                                 gATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCa                                                                    <  1:1188227/100‑1 (MQ=255)
                                  aatCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCa                                                                    <  2:1732754/97‑1 (MQ=255)
                                   ttCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGc                                                                  >  1:206514/1‑100 (MQ=255)
                                          aTCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGc                                                           <  1:542982/100‑1 (MQ=255)
                                            cTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCgg                                                         <  2:980762/100‑1 (MQ=255)
                                             tGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGa                                                        <  2:1368642/100‑1 (MQ=255)
                                              gCATCGCATTCTGTCTAACTTTTCAGTTGGTTATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGAc                                                       <  2:1939202/100‑1 (MQ=255)
                                              gCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGAc                                                       <  2:1820954/100‑1 (MQ=255)
                                              gCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGAc                                                       <  2:1512464/100‑1 (MQ=255)
                                               cATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGAcc                                                      <  1:1597218/100‑1 (MQ=255)
                                                aTCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCa                                                     <  2:1585692/100‑1 (MQ=255)
                                                  cGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCata                                                   >  2:1175203/1‑100 (MQ=255)
                                                  cGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCata                                                   <  1:187291/100‑1 (MQ=255)
                                                    cATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATAtg                                                 >  1:692900/1‑100 (MQ=255)
                                                      tccTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATAtgct                                               >  1:1378221/3‑100 (MQ=255)
                                                        cTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATAtgctgc                                             <  1:145664/100‑1 (MQ=255)
                                                         tGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCa                                            <  2:459251/100‑1 (MQ=255)
                                                             tAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTg                                        >  1:1365945/1‑100 (MQ=255)
                                                               aCGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGc                                      >  1:398984/1‑100 (MQ=255)
                                                                 gTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCtt                                    <  1:333207/100‑1 (MQ=255)
                                                                  tGGCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCttt                                   >  2:198534/1‑100 (MQ=255)
                                                                    gCAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCttttt                                 <  2:356218/100‑1 (MQ=255)
                                                                     cAGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCtttttt                                >  1:509709/1‑100 (MQ=255)
                                                                      aGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTc                               <  1:1823904/100‑1 (MQ=255)
                                                                      aGTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTc                               <  1:676668/100‑1 (MQ=255)
                                                                       gTGGGTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCCGTTTCAGCAAGGCGCGGAGCATATGCTGCAACTGGCTTTTTTCg                              >  2:1872567/1‑100 (MQ=255)
                                                                         gggTGATAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGc                            <  1:1175209/100‑1 (MQ=255)
                                                                              aTAGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGc                       >  2:240907/1‑100 (MQ=255)
                                                                                aGCAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCa                     <  2:135212/100‑1 (MQ=255)
                                                                                  cAATCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAAt                   >  1:398292/1‑100 (MQ=255)
                                                                                   aaTCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAAta                  <  1:1482867/100‑1 (MQ=255)
                                                                                   aaTCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAAta                  >  2:229515/1‑100 (MQ=255)
                                                                                   aaTCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAAta                  <  2:1100437/100‑1 (MQ=255)
                                                                                     tCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAAtacc                >  2:605352/1‑100 (MQ=255)
                                                                                     tCGCCAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAAtacc                <  1:401987/100‑1 (MQ=255)
                                                                                       gCCAGCGCATCGGCGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAAtaccta              >  1:83844/1‑100 (MQ=255)
                                                                                        ccAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAAtacctac             >  1:668622/1‑100 (MQ=255)
                                                                                         cAGCGCATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAAtacctacc            <  1:1110358/100‑1 (MQ=255)
                                                                                           gcgcATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAATACCTaccac          <  2:1151870/100‑1 (MQ=255)
                                                                                              cATCGACGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAATACCTACCACCGt       >  2:1821169/1‑100 (MQ=255)
                                                                                                   gcggcATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCt  <  1:1981353/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGCGATTCGCTCATCTGCATCGCATTCTGACTAACGTGGCAGTGGGTGATAGCAATCGCCAGCGCATCGGCGGCATCCGCCTGTGGATTAGCGGGCAGTTTCAGCAAGGTGCGGACCATATGCTGCACCTGGCTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCT  >  NC_000913/1946856‑1947054

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