Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,884,194 A→C intergenic (+88/‑142) ysdD → / → rpmH protein YsdD/50S ribosomal subunit protein L34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,884,1940AC98.0% 170.6 / NA 50intergenic (+88/‑142)ysdD/rpmHprotein YsdD/50S ribosomal subunit protein L34
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (22/27);  minor base T (1/0);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCAT  >  NC_000913/3884066‑3884304
                                                                                                                                |                                                                                                              
aTCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGaa                                                                                                            <  1:575794/133‑1 (MQ=255)
aTCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGaa                                                                                                            >  1:112912/1‑133 (MQ=255)
    aGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAagag                                                                                                        >  2:405766/1‑133 (MQ=255)
     ggACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAATGTGCGGAGAAGCGCGGAAAATCCGGAAAgaga                                                                                                       >  2:239533/1‑133 (MQ=255)
      gACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGATGCAGAATAGCGCGGTAAATACGGAAAGAGaa                                                                                                      >  2:217756/1‑133 (MQ=255)
          aTCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGa                                                                                                  >  2:67149/1‑133 (MQ=255)
           tCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGACAAGCGCGGTAAGTATGGAAAGGGAATTGAc                                                                                                 >  1:172622/1‑133 (MQ=255)
           tCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGAc                                                                                                 >  2:142739/1‑133 (MQ=255)
               tGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCg                                                                                             <  1:535032/133‑1 (MQ=255)
               tGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCg                                                                                             <  1:313250/133‑1 (MQ=255)
               tGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCg                                                                                             <  1:146206/133‑1 (MQ=255)
                gcgcTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCgg                                                                                            <  2:318580/133‑1 (MQ=255)
                 cgcTTTACCCTTCAGCCCGTATAAACCTCCACACGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGa                                                                                           <  2:172622/133‑1 (MQ=255)
                 cgcTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCg                                                                                             >  2:508157/1‑131 (MQ=255)
                 cgcTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCg                                                                                             <  1:508157/131‑1 (MQ=255)
                 cgcTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGa                                                                                           >  2:29800/1‑133 (MQ=255)
                 cgcTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGa                                                                                           <  2:366609/133‑1 (MQ=255)
                   cTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGa                                                                                           >  1:382896/1‑131 (MQ=255)
                   cTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGa                                                                                           <  2:382896/131‑1 (MQ=255)
                     ttACCCATCAGCCCGTGTAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAgtgt                                                                                       <  1:64800/133‑1 (MQ=255)
                     ttACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAgtgt                                                                                       >  2:83416/1‑133 (MQ=255)
                     ttACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAgtgt                                                                                       <  1:142739/133‑1 (MQ=255)
                      tACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAgtg                                                                                        >  1:379103/1‑131 (MQ=255)
                      tACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAgtg                                                                                        <  2:379103/131‑1 (MQ=255)
                      tACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTa                                                                                      >  2:573537/1‑133 (MQ=255)
                        cccATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACa                                                                                    <  2:501383/133‑1 (MQ=255)
                        cccATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACa                                                                                    <  1:405766/133‑1 (MQ=255)
                         ccATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACaa                                                                                   <  1:544047/133‑1 (MQ=255)
                         ccATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACaa                                                                                   <  2:125032/133‑1 (MQ=255)
                          cATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAat                                                                                  <  1:29800/133‑1 (MQ=255)
                                       aaTCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCgg                                                                     <  1:573537/133‑1 (MQ=255)
                                         tcctccACCCGACGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATcc                                                                       >  2:355401/1‑129 (MQ=255)
                                         tcctccACCCGACGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATcc                                                                       <  1:355401/129‑1 (MQ=255)
                                          cctccACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGAACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGAGTACAATTATTACAATCCGGCct                                                                  >  2:50728/1‑133 (MQ=255)
                                           ctccACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCctc                                                                 <  1:83416/133‑1 (MQ=255)
                                                cccGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTaa                                                            >  1:510356/1‑133 (MQ=255)
                                                    gcgcgcCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTCCAATTATAACAGTCTGGCCTCATTAATCAc                                                        >  1:285594/1‑133 (MQ=255)
                                                          cATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATgg                                                  <  1:67149/133‑1 (MQ=255)
                                                                        cTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATc                                    >  1:63929/1‑133 (MQ=255)
                                                                              tgAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCa                              >  2:12584/1‑133 (MQ=255)
                                                                              tgAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCATCAGTGTCCGTCGTTTCa                              >  1:286349/1‑133 (MQ=255)
                                                                                     gTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTcg                       >  1:343172/1‑133 (MQ=255)
                                                                                      tACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGACGTTGGTGTCCACCGTTTTATTTTTctg                      >  2:282908/1‑131 (MQ=255)
                                                                                       aCATTTTCACTGCGATAATGTGGGGAAAAGGGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTcggc                     <  1:217756/133‑1 (MQ=255)
                                                                                        cATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTcggcg                    <  2:510356/133‑1 (MQ=255)
                                                                                         aTTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTcggcgg                   <  1:506548/133‑1 (MQ=255)
                                                                                            ttCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGata                <  1:282908/133‑1 (MQ=255)
                                                                                             tCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGatat               >  1:42259/1‑133 (MQ=255)
                                                                                                    gAAAAGGTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATaa        <  1:614185/133‑1 (MQ=255)
                                                                                                          gTGCGGAAAAGCGCGGTAAATACGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCAt  <  1:12584/133‑1 (MQ=255)
                                                                                                                                |                                                                                                              
ATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCCCTGCGAAAAGGTGCGGAAAAGCGCGGTAAATAAGGAAAGAGAATTGACTCCGGAGTGTACAATTATTACAATCCGGCCTCTTTAATCACCCATGGCTTCGGTGTCCATCGTTTCATTTTTCGGCGGATATCCAATAAAGCCAT  >  NC_000913/3884066‑3884304

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