Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 418,094 C→T 18.9% 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
*minE418,0940CT18.9% 125.4 / 22.0 58A55T (GCT→ACT) miaBisopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand):  ref base C (27/20);  new base T (7/4);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

TAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCC  >  minE/418025‑418166
                                                                     |                                                                        
tAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGcc                                                                         >  1:1380386/1‑71 (MQ=255)
     aGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTtctc                                                                    <  1:2288230/71‑1 (MQ=255)
       gTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTcgc                                                                  <  1:1038452/71‑1 (MQ=255)
            gATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGAtt                                                             >  1:315252/1‑71 (MQ=255)
                 ttctCTTTTAACAGTTTCCAGCGACCCAACTGATTGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                          >  1:3280437/1‑69 (MQ=255)
                 ttctCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGa                                                           >  1:2194322/1‑68 (MQ=255)
                   ctctTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAgg                                                      >  1:772600/1‑71 (MQ=255)
                    tctTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGgt                                                     >  1:1140102/1‑71 (MQ=255)
                      ttttAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGgtgt                                                   >  1:2532081/1‑71 (MQ=255)
                        ttAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGgtg                                                    <  1:153746/68‑1 (MQ=255)
                        ttAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGgtg                                                    <  1:1175868/68‑1 (MQ=255)
                           aCAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTcagc                                              <  1:3513740/71‑1 (MQ=255)
                           aCAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTcag                                               >  1:1852091/1‑70 (MQ=255)
                           aCAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcag                                               <  1:2120492/70‑1 (MQ=255)
                              gTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTcagcag                                            >  1:3932344/1‑70 (MQ=255)
                              gTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagcag                                            >  1:2305042/1‑70 (MQ=255)
                               tttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCtt                                                                       >  1:2561200/1‑42 (MQ=255)
                                ttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                          >  1:3247712/1‑54 (MQ=255)
                                 tCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTcagcagcag                                         >  1:2512855/1‑70 (MQ=255)
                                 tCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagcagcag                                         >  1:1615283/1‑70 (MQ=255)
                                     gCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACAt                                    <  1:3861403/71‑1 (MQ=255)
                                      cGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATc                                   >  1:983126/1‑71 (MQ=255)
                                          ccAACTGATGGAGGACTTTTTCCTGAGCATTCTCGCGGATTGAGCatg                                                      >  1:1615649/1‑46 (MQ=255)
                                             aCTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttct                            >  1:584972/1‑71 (MQ=255)
                                             aCTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttct                            <  1:1277958/71‑1 (MQ=255)
                                              cTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                          <  1:569519/40‑1 (MQ=255)
                                              cTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                          <  1:240913/40‑1 (MQ=255)
                                              cTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagcagc                                           >  1:3808039/1‑55 (MQ=255)
                                              cTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttctt                           <  1:12708/71‑1 (MQ=255)
                                                 aTGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCg                        >  1:3630372/1‑71 (MQ=255)
                                                 aTGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCg                        >  1:3001757/1‑71 (MQ=255)
                                                 aTGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttctt                           <  1:2238280/68‑1 (MQ=255)
                                                 aTGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCttcttc                          <  1:2941882/69‑1 (MQ=255)
                                                 aTGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCg                        >  1:1073148/1‑71 (MQ=255)
                                                  tGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGAGTTcagcagcag                                         >  1:1439379/1‑53 (MQ=255)
                                                    gAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGcc                      <  1:2577980/70‑1 (MQ=255)
                                                    gAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGcc                      <  1:855586/70‑1 (MQ=255)
                                                      aGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACg                   >  1:1791439/1‑71 (MQ=255)
                                                      aGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttct                            <  1:2089315/62‑1 (MQ=255)
                                                      aGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACg                   <  1:2376508/71‑1 (MQ=255)
                                                       gACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGtt                                                  <  1:1760460/39‑1 (MQ=255)
                                                       gACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGt                  <  1:719346/71‑1 (MQ=255)
                                                       gACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGt                  >  1:2622032/1‑71 (MQ=255)
                                                        aCTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTc                 >  1:758344/1‑71 (MQ=255)
                                                          tttttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCgg               >  1:1991568/1‑71 (MQ=255)
                                                            tttCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGt              >  1:3730590/1‑70 (MQ=255)
                                                            tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGt              <  1:3833203/70‑1 (MQ=255)
                                                            tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGt              >  1:3911314/1‑70 (MQ=255)
                                                            tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGt              >  1:3152349/1‑70 (MQ=255)
                                                             ttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGt              >  1:721199/1‑69 (MQ=255)
                                                               ccTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACg                   >  1:727309/1‑62 (MQ=255)
                                                                cTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGtt         <  1:19890/71‑1 (MQ=255)
                                                                cTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGtt         <  1:1217395/71‑1 (MQ=255)
                                                                cTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGt          <  1:3667299/70‑1 (MQ=255)
                                                                 tGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttc                             <  1:3809149/50‑1 (MQ=255)
                                                                   aGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttctt                           >  1:2510071/1‑50 (MQ=255)
                                                                   aGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTg        >  1:343976/1‑69 (MQ=255)
                                                                   aGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGa       >  1:2330112/1‑70 (MQ=255)
                                                                   aGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGa       >  1:1961010/1‑70 (MQ=255)
                                                                     ccTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTC‑‑CCACGTCGGTCAGTTGATAGcc  <  1:657038/71‑1 (MQ=255)
                                                                     |                                                                        
TAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCC  >  minE/418025‑418166

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: