Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,094,234 T→C T283A (ACA→GCA)  yggR ← putative PilT family AAA+ ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,094,2340TC96.2% 172.6 / ‑2.6 52T283A (ACA→GCA) yggRputative PilT family AAA+ ATPase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (23/27);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

CGTCACAAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGTGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGC  >  NC_000913/3094098‑3094354
                                                                                                                                        |                                                                                                                        
cGTCACAAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGgtgcgt                                                                                                                        >  1:759310/1‑139 (MQ=255)
cGTCACAAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGgtgcgt                                                                                                                        >  1:211457/1‑139 (MQ=255)
    acaAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGTTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAAt                                                                                                                    <  2:759310/139‑1 (MQ=255)
     caAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATc                                                                                                                   >  2:763266/1‑139 (MQ=255)
       aaCGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCCCCGCGGGTGCGTTAATCAg                                                                                                                 <  2:704493/139‑1 (MQ=255)
        aCGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGc                                                                                                                <  1:763266/139‑1 (MQ=255)
            ccTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCTACCGCGGGTGCGTTAATCAGCAAtt                                                                                                            >  1:334823/1‑139 (MQ=255)
             cTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTc                                                                                                           >  1:1268055/1‑139 (MQ=255)
                cccccACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCaaa                                                                                                        <  1:611139/139‑1 (MQ=255)
                          gTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcg                                                                                              >  1:836308/1‑139 (MQ=255)
                           tGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcgc                                                                                             >  1:1302795/1‑139 (MQ=255)
                             cTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGt                                                                                           >  2:458174/1‑139 (MQ=255)
                                 tAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGtcctt                                                                                       <  2:1302795/139‑1 (MQ=255)
                                       tgctgAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                  <  1:144532/106‑1 (MQ=255)
                                       tgctgAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                  >  2:144532/1‑106 (MQ=255)
                                       tgctgAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGtt                                                                                 >  2:1105134/1‑139 (MQ=255)
                                           gAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATc                                                                             <  2:334823/139‑1 (MQ=255)
                                           gAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATc                                                                             <  1:192579/139‑1 (MQ=255)
                                                        ccccACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCttt                                                                <  2:836308/139‑1 (MQ=255)
                                                         cccACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCttt                                                                >  1:12599/1‑138 (MQ=255)
                                                         cccACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCttt                                                                <  2:12599/138‑1 (MQ=255)
                                                                 ctgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                  <  1:274197/80‑1 (MQ=255)
                                                                 ctgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCa                                                                                                                  >  2:274197/1‑80 (MQ=255)
                                                                 ctgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGc                                                       >  1:1110790/1‑139 (MQ=255)
                                                                  tgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCa                                                                                                >  1:1085506/1‑97 (MQ=255)
                                                                  tgcCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCa                                                                                                <  2:1085506/97‑1 (MQ=255)
                                                                           tGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCa                                                                     <  1:1155236/115‑1 (MQ=255)
                                                                           tGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCa                                                                     >  2:1155236/1‑115 (MQ=255)
                                                                           tGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTa                                             <  2:1110790/139‑1 (MQ=255)
                                                                                  cATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAAt                                                                                                                    >  2:906154/1‑61 (MQ=255)
                                                                                  cATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAAt                                                                                                                    <  1:906154/61‑1 (MQ=255)
                                                                                   aTGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGTGTTAATCAGCAATTCAAATAGCGCCACGCGtcct                                                                                        >  1:527795/1‑88 (MQ=255)
                                                                                   aTGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGTGTTAATCAGCAATTCAAATAGCGCCACGCGtcct                                                                                        <  2:527795/88‑1 (MQ=255)
                                                                                   aTGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAAt                                                                                                                    <  2:405013/60‑1 (MQ=255)
                                                                                   aTGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAAt                                                                                                                    >  1:405013/1‑60 (MQ=255)
                                                                                     gCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAAt                                                                                                                    <  2:1238856/58‑1 (MQ=255)
                                                                                     gCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAAt                                                                                                                    >  1:1238856/1‑58 (MQ=255)
                                                                                     gCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcg                                                                                              >  1:492994/1‑80 (MQ=255)
                                                                                     gCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCAcg                                                                                              <  2:492994/80‑1 (MQ=255)
                                                                                     gCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGc                                   <  2:1268055/139‑1 (MQ=255)
                                                                                       ggTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCTCCACCGCGGGTGCGTTAATCAGc                                                                                                                >  1:857085/1‑60 (MQ=255)
                                                                                       ggTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGc                                                                                                                <  2:857085/60‑1 (MQ=255)
                                                                                       ggTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAAt                                                                                                             >  2:1315960/1‑63 (MQ=255)
                                                                                       ggTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAAt                                                                                                             <  1:1315960/63‑1 (MQ=255)
                                                                                         tAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGt                               >  1:1103438/1‑139 (MQ=255)
                                                                                               gtggGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTc                                                                       >  2:512916/1‑93 (MQ=255)
                                                                                               gtggGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTc                                                                       <  1:512916/93‑1 (MQ=255)
                                                                                                  ggTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGc                      <  1:1105134/139‑1 (MQ=255)
                                                                                                      ttCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACgg                  <  2:1103438/139‑1 (MQ=255)
                                                                                                       tCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                     >  1:438197/1‑103 (MQ=255)
                                                                                                       tCCCTTCGCGAATCAAATTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCAc                                                     <  2:438197/103‑1 (MQ=255)
                                                                                                                      aaTTCCCCACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTgcgc  <  2:1282103/139‑1 (MQ=255)
                                                                                                                             cACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCttt                                                                >  2:626161/1‑70 (MQ=255)
                                                                                                                             cACCGCGGGTGCGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCttt                                                                <  1:626161/70‑1 (MQ=255)
                                                                                                                                        |                                                                                                                        
CGTCACAAACGCCCTTCCCCCACCCGGTGCTGATAACTCTGCTGAAACGTTATCATCCCCACCTGCTGCCCGGTTTGAATAACATGCGGTAACTGGTGGGTTTTCCCTTCGCGAATCAAATTCCCCACCGCGGGTGTGTTAATCAGCAATTCAAATAGCGCCACGCGTCCTTCCTGTTTATCCACTTCCAGCTTTTGTGACAGCACTGCCCGTAAACTACCTGCCAGTTGATTACGCACGGGGTCTTTTTCCTGCGC  >  NC_000913/3094098‑3094354

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: ‑