Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,152,876 G→A P85S (CCG→TCG)  yegL ← VMA domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,152,8760GA96.3% 158.2 / ‑1.7 54P85S (CCG→TCG) yegLVMA domain protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (31/21);  total (32/22)
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.15e-01

GAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGTTCATCACGAAAGGTAACCAATCCGGCGTTA  >  NC_000913/2152748‑2153009
                                                                                                                                |                                                                                                                                     
gAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTa                                                                                                                             >  2:112776/1‑139 (MQ=255)
            gCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTa                                                                                                                 <  1:112776/139‑1 (MQ=255)
               ataataCGAAATACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAg                                                                                                              >  1:608274/1‑139 (MQ=255)
                     cGAAATACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGgttgt                                                                                                        <  1:527745/139‑1 (MQ=255)
                       aaaTACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTc                                                                                                      <  1:885587/139‑1 (MQ=255)
                        aaTACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTcc                                                                                                     >  2:374126/1‑139 (MQ=255)
                          tACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCcac                                                                                                   >  1:764210/1‑139 (MQ=255)
                             cATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATg                                                                                                >  1:304928/1‑139 (MQ=255)
                                       aTACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCCGCTCAAGGGGTGTTCCACATGCCCCGCGCCg                                                                                      >  2:693283/1‑139 (MQ=255)
                                           tCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACg                                                                                  >  1:223184/1‑139 (MQ=255)
                                                  tCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAAt                                                                           >  1:805653/1‑139 (MQ=255)
                                                  tCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGCCGAAAAAAAATTAGCGGCGCTGGTAAACGGTTGTTCCACATGCACCGGGCCGAACGTCACAAt                                                                           >  1:636842/1‑139 (MQ=255)
                                                    cGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGc                                                                         >  1:539158/1‑139 (MQ=255)
                                                         ctctACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGt                                                                    >  1:795211/1‑139 (MQ=255)
                                                         ctctACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGt                                                                    >  2:398372/1‑139 (MQ=255)
                                                            tACCATATCCAGGGCTTTGGTAATGGCAGCACCCATTGGTGTATCGCCCTGGGCAAACAGGCTGGGCGAAAAAAAATTAACGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTcc                                                                 >  2:639227/1‑139 (MQ=255)
                                                                 tatCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCActct                                                            >  1:642052/1‑139 (MQ=255)
                                                                   tCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCtttt                                                         >  2:35922/1‑139 (MQ=255)
                                                                    ccAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTca                                                                               >  2:289441/1‑117 (MQ=255)
                                                                    ccAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTca                                                                               <  1:289441/117‑1 (MQ=255)
                                                                      aGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTa                                                       >  1:262147/1‑139 (MQ=255)
                                                                        ggCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAgc                                                                                                                           <  1:148204/69‑1 (MQ=255)
                                                                        ggCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAgc                                                                                                                           >  2:148204/1‑69 (MQ=255)
                                                                        ggCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGcc                                                                                       >  1:1014468/1‑105 (MQ=255)
                                                                        ggCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGcc                                                                                       <  2:1014468/105‑1 (MQ=255)
                                                                        ggCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTaaa                                                     >  2:842071/1‑139 (MQ=255)
                                                                        ggCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTaaa                                                     >  1:25600/1‑139 (MQ=255)
                                                                          cTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGc                                                   <  2:820029/139‑1 (MQ=255)
                                                                            ttGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAg                                                <  1:35922/139‑1 (MQ=255)
                                                                             tGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAg                                                <  1:748075/139‑1 (MQ=255)
                                                                               gTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCg                                              <  1:52229/139‑1 (MQ=255)
                                                                                    ggCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCg                                              <  2:45261/134‑1 (MQ=255)
                                                                                    ggCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCg                                              >  1:45261/1‑134 (MQ=255)
                                                                                      ccGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCagc                                       >  1:974826/1‑139 (MQ=255)
                                                                                         cACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAg                                                <  2:885577/127‑1 (MQ=255)
                                                                                         cACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAg                                                >  1:885577/1‑127 (MQ=255)
                                                                                            ccATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAg                                 <  2:533203/139‑1 (MQ=255)
                                                                                            ccATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAg                                 <  1:374126/139‑1 (MQ=255)
                                                                                             cATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGt                                <  2:65367/139‑1 (MQ=255)
                                                                                                        cGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGTtc                              >  1:728020/1‑130 (MQ=255)
                                                                                                        cGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGTtc                              <  2:728020/130‑1 (MQ=255)
                                                                                                        cGCCCTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGTTCATCACGaaa                     <  2:539158/139‑1 (MQ=255)
                                                                                                          cccTGGGCAAACAGGATGGGCGAAAAAAAATTAgcggcg                                                                                                                       <  2:205983/39‑1 (MQ=255)
                                                                                                          cccTGGGCAAACAGGATGGGCGAAAAAAAATTAgcggcg                                                                                                                       >  1:205983/1‑39 (MQ=255)
                                                                                                          cccTGGGCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACTTCACAATGCCGAGTTCCACTCTTTTTAGAGCCAGCGGATCAGCAAGCAGTTCATCACGAAAgg                   >  1:637124/1‑139 (MQ=255)
                                                                                                              gggCAAACAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGTTCATCACGAAAGGTAAc               >  2:317862/1‑139 (MQ=255)
                                                                                                                    aCAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCAc                                                               <  2:886546/85‑1 (MQ=255)
                                                                                                                    aCAGGATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCAc                                                               >  1:886546/1‑85 (MQ=255)
                                                                                                                        gATGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGTTCATCACGAAAGGTAACCAATCCGGCg     <  2:642052/139‑1 (MQ=255)
                                                                                                                         aTGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGcc                                                                        <  1:430451/71‑1 (MQ=255)
                                                                                                                         aTGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGcc                                                                        >  2:430451/1‑71 (MQ=255)
                                                                                                                         aTGGGCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGTTCATCACGAAAGGTAACCAATCCGGCGt    >  1:606901/1‑139 (MQ=255)
                                                                                                                           gggCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGTTCATCACGAAAGGTAACCAATCCGGCGTTa  >  1:618559/1‑139 (MQ=255)
                                                                                                                           gggCGAAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGTTCATCACGAAAGGTAACCAATCCGGCGTTa  <  2:618559/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                     
GAAAATCCACGGGCGATAATACGAAATACCATTTGCACGATACTCACGTTTCCGCTCCTCTACCATATCCAGGGCTTTGGTAATGGCCGCACCCATTGGTGTATCGCCCTGGGCAAACAGGATGGGCGGAAAAAAATTAGCGGCGCTGGTAAAGGGTTGTTCCACATGCACCGGGCCGAACGTCACAATGCCGAGTTCCACTCTTTTTAAAGCCAGCGGATCAGCAAGCAGTTCATCACGAAAGGTAACCAATCCGGCGTTA  >  NC_000913/2152748‑2153009

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