Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 695,215 C→T 31.8% A55T (GCT→ACT)  miaB ← isopentenyl‑adenosine A37 tRNA methylthiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb695,2150CT31.8% 77.8 / 50.4 60A55T (GCT→ACT) miaBisopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base T (12/7);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

AATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGAT  >  W3110S.gb/695147‑695283
                                                                    |                                                                    
aaTCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCt                                                                    >  1:3495205/1‑71 (MQ=255)
    aGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTtctc                                                                <  1:1767932/71‑1 (MQ=255)
        cTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTcgcg                                                             >  1:5046005/1‑70 (MQ=255)
        cTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTACTGAGCCTTCTcgcg                                                             >  1:3074331/1‑70 (MQ=255)
                ttctCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAg                                                      >  1:3346974/1‑69 (MQ=255)
                ttctCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAg                                                      >  1:156679/1‑69 (MQ=255)
                ttctCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                      >  1:555600/1‑69 (MQ=255)
                ttctCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                      >  1:2689640/1‑69 (MQ=255)
                  ctctTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGc                                                     <  1:2773992/68‑1 (MQ=255)
                      tttAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGc                                                     <  1:3060633/64‑1 (MQ=255)
                      tttAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGc                                                     <  1:3641307/64‑1 (MQ=255)
                       ttAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGtt                                              >  1:3576330/1‑70 (MQ=255)
                        tAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGtt                                              >  1:1162057/1‑69 (MQ=255)
                        tAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGtt                                              >  1:748868/1‑69 (MQ=255)
                          aCAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTcagc                                          <  1:1885621/71‑1 (MQ=255)
                          aCAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagc                                          <  1:2470147/71‑1 (MQ=255)
                             gTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCgg                                                            <  1:3751540/50‑1 (MQ=255)
                             gTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCgg                                                            <  1:2411861/50‑1 (MQ=255)
                              tttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTcag                                           <  1:1902261/66‑1 (MQ=255)
                              tttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTcag                                           <  1:4762280/66‑1 (MQ=255)
                              tttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                      <  1:2505378/55‑1 (MQ=255)
                              tttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                      <  1:2547945/55‑1 (MQ=255)
                              tttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                      <  1:1345750/55‑1 (MQ=255)
                                tCCAGCGACCCAACGGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagcagcag                                     >  1:3203795/1‑70 (MQ=255)
                                 ccAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                      <  1:4539874/52‑1 (MQ=255)
                                  cAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAgcagcagca                                   <  1:1746275/70‑1 (MQ=255)
                                      gACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagcagcag                                     <  1:4016702/64‑1 (MQ=255)
                                      gACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATc                               <  1:3094969/70‑1 (MQ=255)
                                      gACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATc                               <  1:2560249/70‑1 (MQ=255)
                                        cccAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGc                            <  1:499728/71‑1 (MQ=255)
                                          cAACTGATGGAAGACTTTTTCCTGAGACTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGctt                          >  1:3681733/1‑71 (MQ=255)
                                           aaCTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTcagcagca                                      >  1:1251540/1‑58 (MQ=255)
                                           aaCTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTAcagcagc                                       >  1:4310394/1‑57 (MQ=255)
                                            aCTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttct                        >  1:304058/1‑71 (MQ=255)
                                            aCTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttct                        >  1:789429/1‑71 (MQ=255)
                                            aCTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttct                        >  1:2258623/1‑71 (MQ=255)
                                            aCTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttct                        >  1:4019436/1‑71 (MQ=255)
                                              tGATGGAAGACTTTTTCCTGAGTCTTCTCGGGGATTGAGCAGGTGTTcagca                                         >  1:2943248/1‑52 (MQ=255)
                                              tGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCttcttc                      >  1:4028872/1‑71 (MQ=255)
                                               gATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTcc                     >  1:1659082/1‑71 (MQ=255)
                                                 tGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCATCAGCACATCCGcttct                        <  1:546750/66‑1 (MQ=255)
                                                  ggAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGc                   >  1:975773/1‑70 (MQ=255)
                                                   gAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCGCATCCGCTTCTTCCGcc                  >  1:519721/1‑70 (MQ=255)
                                                   gAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGcc                  >  1:1136143/1‑70 (MQ=255)
                                                    aaGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCAc                >  1:2142581/1‑71 (MQ=255)
                                                     aGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACg               >  1:1179448/1‑71 (MQ=255)
                                                       aCTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGc                   <  1:3036030/65‑1 (MQ=255)
                                                        cTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCg            <  1:1538660/71‑1 (MQ=255)
                                                         tttttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTc                                             <  1:5000465/37‑1 (MQ=255)
                                                           tttCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGctt                          >  1:3964653/1‑54 (MQ=255)
                                                           tttCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGctt                          >  1:4504632/1‑54 (MQ=255)
                                                           tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagcagca                                      >  1:3612965/1‑42 (MQ=255)
                                                           tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATc                               >  1:4751169/1‑49 (MQ=255)
                                                           tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttct                        >  1:2165875/1‑56 (MQ=255)
                                                           tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGt          >  1:2884032/1‑70 (MQ=255)
                                                           tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGt          >  1:4314517/1‑70 (MQ=255)
                                                            ttCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGt          >  1:3928938/1‑69 (MQ=255)
                                                                  aGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTg    <  1:2243115/69‑1 (MQ=255)
                                                                  aGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTg    <  1:1296770/69‑1 (MQ=255)
                                                                  aGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGAt  <  1:3395628/71‑1 (MQ=255)
                                                                    |                                                                    
AATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGAT  >  W3110S.gb/695147‑695283

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