Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,726,892 C→A intergenic (+10/+32) wecH → / ← yiaA O‑acetyltransferase WecH/conserved inner membrane protein YiaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,726,8920CA93.9% 106.0 / 1.3 33intergenic (+10/+32)wecH/yiaAO‑acetyltransferase WecH/conserved inner membrane protein YiaA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (17/14);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CCGGGGAATTGAGCTTAAAAATTGGCCAATACTGGATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCACCGCTACTGCACCAAACAAGCTTA  >  NC_000913/3726756‑3727024
                                                                                                                                        |                                                                                                                                      
ccGGGGTATTGAGCTTAAAAATTGGCCAATACTGGATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGt                                                                                                                                     <  2:279136/140‑1 (MQ=255)
ccGGGGAATTGAGCTTAATAATTGGCCAATACTGGATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGt                                                                                                                                     <  2:374240/140‑1 (MQ=255)
         tGAGCTTAAAAATTGGCCAATACTGGATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGa                                                                                                                            <  1:289444/140‑1 (MQ=255)
            gCTTAAAAATTGGCCAATACTGGATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTg                                                                                                                         >  2:466505/1‑140 (MQ=255)
                       ggCCAATACTGGATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTAc                                                                                                              >  2:301926/1‑140 (MQ=255)
                              aCTGGATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGtt                                                                                                       <  2:456154/140‑1 (MQ=255)
                                tGGATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTAt                                                                                                     <  1:301926/140‑1 (MQ=255)
                                  gATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTc                                                                                                   >  2:271006/1‑140 (MQ=255)
                                  gATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTc                                                                                                   >  2:480896/1‑140 (MQ=255)
                                    tattatTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGc                                                                                                 <  2:304020/140‑1 (MQ=255)
                                     attattTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGAAAAAGCCAGGTCACATTGGACCGACCGGGCTTACGTGAGTTATTCGCt                                                                                                <  2:411417/140‑1 (MQ=255)
                                             gATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTAttctt                                                                                        <  1:466505/140‑1 (MQ=255)
                                               tCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTcc                                                                                      <  1:288390/140‑1 (MQ=255)
                                                  tttGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTgg                                                                                   >  1:320185/1‑140 (MQ=255)
                                                  tttGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTgg                                                                                   >  1:433450/1‑140 (MQ=255)
                                                   ttGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTggg                                                                                  >  2:325689/1‑140 (MQ=255)
                                                    tGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGt                                                                                 <  2:364794/140‑1 (MQ=255)
                                                        cgACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAAcc                                                                             <  2:320185/140‑1 (MQ=255)
                                                           cGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCtgt                                                                          >  1:262037/1‑140 (MQ=255)
                                                                                           gTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACg                                          >  2:206059/1‑140 (MQ=255)
                                                                                             aCAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGaa                                        <  1:325689/140‑1 (MQ=255)
                                                                                                     tCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTc                              >  2:152894/1‑140 (MQ=255)
                                                                                                      cGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCt                             <  1:460674/140‑1 (MQ=255)
                                                                                                           gAAACAGATTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCAcc                        >  1:128835/1‑140 (MQ=255)
                                                                                                                  aTTAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCACCGCTACTg                 >  1:12041/1‑140 (MQ=255)
                                                                                                                    tAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCACCGCTACTGCa               <  2:12041/140‑1 (MQ=255)
                                                                                                                    tAGTGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCACCGCTACTGCa               <  2:101738/140‑1 (MQ=255)
                                                                                                                       tGAGTTAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTTTCTGCACCGCTACTGCAc              >  2:286431/1‑140 (MQ=255)
                                                                                                                           ttAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCACCGCTACTGCACCAAACa        >  2:387567/1‑140 (MQ=255)
                                                                                                                           ttAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCACCGCTACTGCACCAAACa        >  1:176978/1‑140 (MQ=255)
                                                                                                                           ttAAGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCACCGCTACTGCACCAAACa        >  2:77781/1‑140 (MQ=255)
                                                                                                                              aGTAAAAGCCAGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCACCGCTACTGCACCAAACAAGc     >  1:250016/1‑140 (MQ=255)
                                                                                                                                 aaaaGCCAGGGCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCACCGCTACTGCACCAAACAAGCTTa  >  1:34224/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                      
CCGGGGAATTGAGCTTAAAAATTGGCCAATACTGGATATTATTTGGATCTTTTGCGCGACGTTGGCAGCGAGTTTGTTACTTTCTATGCTGGTACAACGAATCGACAGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATA‑‑TTTTTCTGCACCGCTACTGCACCAAACAAGCTTA  >  NC_000913/3726756‑3727024

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