Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,855,194 A→T 35.5% I90N (ATC→AAC)  ygbA ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,855,1940AT35.5% 42.5 / 48.2 48I90N (ATC→AAC) ygbAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/10);  new base T (4/13);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 5.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

ACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCTTTG  >  W3110S.gb/2855130‑2855260
                                                                |                                                                  
aCGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCg                                                                 <  1:793615/68‑1 (MQ=255)
        tttCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCg                                                                 >  1:3265564/1‑60 (MQ=255)
            cGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCcaca                                                     <  1:3097175/68‑1 (MQ=255)
             ggCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCcaca                                                     <  1:2835584/67‑1 (MQ=255)
              gCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGTTCg                                                                 <  1:862992/54‑1 (MQ=255)
              gCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGTTCg                                                                 <  1:2673838/54‑1 (MQ=255)
              gCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGTTCg                                                                 >  1:298684/1‑54 (MQ=255)
              gCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGTTCg                                                                 >  1:430564/1‑54 (MQ=255)
              gCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGTTCg                                                                 >  1:1778694/1‑54 (MQ=255)
                aGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCAtt                                                <  1:816171/69‑1 (MQ=255)
                aGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCAtt                                                <  1:3261248/69‑1 (MQ=255)
                   tCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATtctc                                             >  1:1616848/1‑69 (MQ=255)
                            ggACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTcggcccggc                                     >  1:621967/1‑68 (MQ=255)
                             gACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGAt                                                              <  1:3480415/42‑1 (MQ=255)
                              aCGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCa                                  >  1:871509/1‑69 (MQ=255)
                                   tGTCATCGATGAGATGACGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCa                             <  1:2637601/69‑1 (MQ=255)
                                   tGTCATCGATGAGATGACGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCa                             <  1:2845198/69‑1 (MQ=255)
                                    gTCATCGATGAGATGACGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCggcccggccc                                   >  1:1832626/1‑62 (MQ=255)
                                      cATCGATGAGATGACGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCAt                            <  1:2895966/67‑1 (MQ=255)
                                      cATCGATGAGATGACGCACGGTTAAGATCGGATGCCGCCACAGCATTCTCGGCCCGGCCCAACGCATaa                          >  1:131253/1‑69 (MQ=255)
                                       aTCGATGAGATGACGCACGGTTAAGATCGGATGACGGCACAGCATTCTCggcccggccc                                   >  1:429238/1‑59 (MQ=255)
                                             gagatgaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt                         >  1:370288/1‑63 (MQ=255)
                                             gagatgaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt                         >  1:371049/1‑63 (MQ=255)
                                             gagatgaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt                         >  1:2736970/1‑63 (MQ=255)
                                             gagatcaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAAt                         >  1:1326228/1‑63 (MQ=255)
                                              agatgaCGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTc                  <  1:2828775/69‑1 (MQ=255)
                                              agatgaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTc                  <  1:2784073/69‑1 (MQ=255)
                                              agatgaCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTc                  <  1:1210727/69‑1 (MQ=255)
                                               gatgaCGCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTc                  <  1:2026483/68‑1 (MQ=255)
                                                   aCGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCAtctc             >  1:1198865/1‑69 (MQ=255)
                                                     gCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCtt           <  1:2618488/69‑1 (MQ=255)
                                                     gCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCtt           <  1:1349684/69‑1 (MQ=255)
                                                     gCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCtt           <  1:790530/69‑1 (MQ=255)
                                                     gCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCtt           <  1:1439811/69‑1 (MQ=255)
                                                     gCACGGTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCtt           <  1:1548255/69‑1 (MQ=255)
                                                        cGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTc                  >  1:2636434/1‑59 (MQ=255)
                                                          gTTAAGTTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGc      <  1:44252/69‑1 (MQ=255)
                                                          gTTAAGATCGGATGACGCCACAGCATTCTCggcccggccc                                   <  1:1847964/40‑1 (MQ=255)
                                                          gTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGc      >  1:1687835/1‑69 (MQ=255)
                                                          gTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGc      >  1:729302/1‑69 (MQ=255)
                                                          gTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGc      >  1:733577/1‑69 (MQ=255)
                                                          gTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGc      >  1:2561935/1‑69 (MQ=255)
                                                              aGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGc                              >  1:602898/1‑41 (MQ=255)
                                                              aGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCtt    >  1:525077/1‑67 (MQ=255)
                                                              aGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCtt    >  1:1367089/1‑67 (MQ=255)
                                                              aGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCtt    >  1:3134851/1‑67 (MQ=255)
                                                              aGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCtt    >  1:2774758/1‑67 (MQ=255)
                                                               gATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCTTTg  <  1:2517998/68‑1 (MQ=255)
                                                                aTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCttt   >  1:2556003/1‑66 (MQ=255)
                                                                aTCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCTTTg  <  1:1212662/67‑1 (MQ=255)
                                                                |                                                                  
ACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCATTCTCGGCCCGGCCCAACGCATAATCTGTTTCATCTCTTCACGCTTTG  >  W3110S.gb/2855130‑2855260

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: