Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 61,461 G→A Q602* (CAG→TAG)  rapA ← RNA polymerase remodeling/recycling factor ATPase, RNA polymerase‑associated, ATP‑dependent RNA translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091361,4610GA95.7% 60.3 / NA 23Q602* (CAG→TAG) rapARNA polymerase remodeling/recycling factor ATPase, RNA polymerase‑associated, ATP‑dependent RNA translocase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (12/10);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATACCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGAACCGATTTCTGAGCACAGCAGTACCTGTGCGCCGGT  >  NC_000913/61326‑61602
                                                                                                                                       |                                                                                                                                             
ataCGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATACCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGATa                                                                                                                                        <  2:199026/143‑1 (MQ=255)
            aaaTAGTGCCTCCGGTCCGTCAGGTGTGCTCAAATGCATCCAGACCCTCGTGGTACCACCCAACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGATACGATCCAgacga                                                                                                                            <  1:90443/143‑1 (MQ=255)
                         ggTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATACCAGCGCACCAGCACAGACTGAGAGGTTTTCTCCAGATAAGGAACATGGATCTGAATAACGTGGGCCTAGCCGGTACGATCCAGACGACGAATACgctgct                                                                                                               >  2:485665/1‑143 (MQ=255)
                           tCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATACCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGATCCGATCCAGACGACCAATACTCTGCTcc                                                                                                             >  1:63218/1‑143 (MQ=255)
                                        cTCAACTGCATCCAGACCTTCGTGATACCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGt                                                                                                <  2:267978/143‑1 (MQ=255)
                                                 aTCCAGGCCGGCGTGATACCATGGCACCAGCACCGTCTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGTCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCa                                                                                       <  2:215652/143‑1 (MQ=255)
                                                   ccAGACCTTCGTGATACCAGCGAACCAGCACCGACTGAGCGGGTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAgg                                                                                     <  2:322731/143‑1 (MQ=255)
                                                    cAGACCTTCGTGATACCATCGCGCCAGCACCGACTGTGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATCTCGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGt                                                                                    <  2:279387/143‑1 (MQ=255)
                                                     aGACCTTCGTGATACCAGCGCACCAGCACCGACTGAGCGGGTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGCTACGATCCAGACGACCAATACGCTGCCGCAGGCGATCCGGGCTGAATGGCACGTc                                                                                   >  1:2107/1‑143 (MQ=255)
                                                         cTTCGTGATACCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAAc                                                                               >  1:145361/1‑143 (MQ=255)
                                                                            cAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGa                                                            >  1:340354/1‑143 (MQ=255)
                                                                                 ccGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGATACGATCCAGACGGCCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTGAAACATCACCATGTGGCTGGCGAACTgg                                                       >  1:181693/1‑143 (MQ=255)
                                                                                   tgtttggattttttCTCCAGATAAGCCACATGGATCTGAATATCGTGCGCCTAACCGATACGATCCAGACGACCAATACGCTTCTCCAGTAGATCCTGTTTGAATTGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGaa                                                     <  2:62926/133‑1 (MQ=255)
                                                                                            gTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGGACTGGACGTTACGTcc                                            >  1:275770/1‑143 (MQ=255)
                                                                                              tttCTCCAGATAACGCACAAGGGTCGGAAGATCGCGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCtt                                          <  1:455174/143‑1 (MQ=255)
                                                                                                   ccAGATAAGGCACATGGATCTGAATATCGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGaa                                     >  1:370187/1‑143 (MQ=255)
                                                                                                            ccgcaTGGCTCTGGTTATCGTGCGCCTAGCCGCTACTGTCCAGACGTCCTATACGCTGCTCCAGTACGTGCGAGTTGAAGGGCAATTTAATCTGCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGAACCGATTTCt                            <  2:422877/140‑1 (MQ=255)
                                                                                                                     tCTGAATATCGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGCTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACCGGAAGTTACGGCCTTCAGAACCGATTCCTGAGCAccgc                   >  2:243042/1‑143 (MQ=255)
                                                                                                                      cTGAATATCGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGAACCGATTTCTGAGCAcagca                  <  1:186145/143‑1 (MQ=255)
                                                                                                                       tGAATATCGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGAACCGATTTCTGAGCAcagcag                 >  2:476900/1‑143 (MQ=255)
                                                                                                                              cGTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGAACCGATTTCTGAGCACAGCAGTACCtgt          >  1:161521/1‑143 (MQ=255)
                                                                                                                               gTGCGCCTAGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCCGCACCGATTTCTGCGCTCCGCACTACAtgtg         >  1:177953/1‑143 (MQ=255)
                                                                                                                                      tGGCCGCTACGCGCCAGTCGCCCAATCCGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGAACCGATTTCTGAGCACAGCAGTACCTGTGCGCCGgt  <  2:251696/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                             
ATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATACCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACATCACCATGTGGCTGGCGAACTGGAAGTTACGTCCTTCAGAACCGATTTCTGAGCACAGCAGTACCTGTGCGCCGGT  >  NC_000913/61326‑61602

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