Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,315,139 T→A 27.0% D244V (GAT→GTT)  infB ← fused protein chain initiation factor 2, IF2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,315,1390TA27.0% 100.7 / 37.9 59D244V (GAT→GTT) infBfused protein chain initiation factor 2, IF2
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base A (5/11);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTGCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAG  >  W3110S.gb/3315071‑3315205
                                                                    |                                                                  
tgttgCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACg                                                                  >  1:2007345/1‑71 (MQ=255)
 gttgCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGc                                                                 >  1:479489/1‑71 (MQ=255)
 gttgCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACACGGCCACGGCCGCCTTCGACTTCACGATCGCt                                                                >  1:2740660/1‑71 (MQ=255)
    gCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCttt                                                              >  1:3393934/1‑71 (MQ=255)
    gCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCtt                                                               >  1:4701588/1‑70 (MQ=255)
         tcttcGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGtct                                                         >  1:3252607/1‑71 (MQ=255)
             cGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGtcttct                                                      <  1:1778460/70‑1 (MQ=255)
                     ctgctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTc                                                             <  1:3512559/55‑1 (MQ=255)
                     ctgctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTc                                                             <  1:3506797/55‑1 (MQ=255)
                        ctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCa                                          >  1:3088797/1‑71 (MQ=255)
                        ctTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCa                                          >  1:421836/1‑71 (MQ=255)
                           tCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATg                                        <  1:2872160/70‑1 (MQ=255)
                            cgcgTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGGGCATGtt                                      >  1:3869779/1‑71 (MQ=255)
                            cgcgTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGt                                                           >  1:2189097/1‑50 (MQ=255)
                             gcgTTACGACCTCGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTg                                     <  1:3107840/71‑1 (MQ=255)
                             gcgTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGtt                                      >  1:4519627/1‑70 (MQ=255)
                              cgTTACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCa                                          <  1:436932/65‑1 (MQ=255)
                               gTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                   <  1:1956235/71‑1 (MQ=255)
                               gTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                   <  1:3720215/71‑1 (MQ=255)
                                ttACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCtt                                                               <  1:2133196/42‑1 (MQ=255)
                                 tACGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGtctt                                                        >  1:118328/1‑48 (MQ=255)
                                  aCGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTATTCTGCCTGGCGAGCATGTTGAGAag                                >  1:4143146/1‑71 (MQ=255)
                                  aCGACCACGGCCACGGCCGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGa                                             >  1:1650567/1‑58 (MQ=255)
                                   cGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCttt                                                              >  1:4561664/1‑40 (MQ=255)
                                   cGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCttt                                                              >  1:590660/1‑40 (MQ=255)
                                   cGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAagt                               <  1:1292450/71‑1 (MQ=255)
                                      ccacggccacggCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAagtagt                            <  1:4155192/71‑1 (MQ=255)
                                       cacggccacggcCGCCTTCGACTTCACGAACGCTTTCGtc                                                          <  1:347611/40‑1 (MQ=255)
                                         cggccacggccGCCTTCGACTTCACGATCGCTTTc                                                             <  1:3630628/35‑1 (MQ=255)
                                          ggccacggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGAc                         >  1:3245778/1‑70 (MQ=255)
                                          ggccacggccGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACg                        >  1:4186334/1‑71 (MQ=255)
                                           gccacggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                   <  1:997337/59‑1 (MQ=255)
                                           gccacggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                   <  1:2631034/59‑1 (MQ=255)
                                             cacggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGCCGTGa                     >  1:4306674/1‑71 (MQ=255)
                                               cggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATa                   <  1:4423616/71‑1 (MQ=255)
                                                ggccGCCTTCGACTTCACGAACGCTTTCGTCTTCTGCCTGGCGAGCa                                          <  1:3396286/47‑1 (MQ=255)
                                                 gccgccTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAagta                              >  1:4775665/1‑58 (MQ=255)
                                                       ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgaga                                  >  1:4655079/1‑48 (MQ=255)
                                                       ttCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGAc                         >  1:3251654/1‑57 (MQ=255)
                                                        tCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTAGAGAAGTAGTGAc                         >  1:2382029/1‑56 (MQ=255)
                                                         cGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCa                                          >  1:1859400/1‑38 (MQ=255)
                                                         cGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGaa         >  1:1165167/1‑71 (MQ=255)
                                                         cGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGa          <  1:676893/70‑1 (MQ=255)
                                                          gACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGa                     >  1:1252423/1‑58 (MQ=255)
                                                          gACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAAt        <  1:781176/71‑1 (MQ=255)
                                                            cTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCt      <  1:1781305/71‑1 (MQ=255)
                                                             ttCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTgaga                                  <  1:2657846/42‑1 (MQ=255)
                                                             ttCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGa                     <  1:738129/55‑1 (MQ=255)
                                                             ttCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGa                     <  1:1105665/55‑1 (MQ=255)
                                                              tCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCg               <  1:1436806/60‑1 (MQ=255)
                                                               cACGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCtt     <  1:2404010/69‑1 (MQ=255)
                                                                aCGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGTTAATCGc              >  1:1775776/1‑59 (MQ=255)
                                                                aCGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAg  <  1:82048/71‑1 (MQ=255)
                                                                 cGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                   <  1:4354334/37‑1 (MQ=255)
                                                                 cGAACGCTTTCGTCTTCTGCCTGGCGAGCATGTTgag                                   <  1:855521/37‑1 (MQ=255)
                                                                  gATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATa                   <  1:4667004/52‑1 (MQ=255)
                                                                   aTCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAagtagt                            <  1:2590553/42‑1 (MQ=255)
                                                                   aTCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAagtagt                            <  1:2644539/42‑1 (MQ=255)
                                                                   aTCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAg  <  1:1509650/68‑1 (MQ=255)
                                                                    |                                                                  
TGTTGCCTTTCTTCGGACGCGCTGCTTTCGCGTTACGACCACGGCCACGGCCGCCTTCGACTTCACGATCGCTTTCGTCTTCTGCCTGGCGAGCATGTTGAGAAGTAGTGACGTGATAATCGCTGGAATCTTCAG  >  W3110S.gb/3315071‑3315205

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: