Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 695,215 C→T 10.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
*AP009048695,2150CT10.9% 133.5 / 6.5 46A55T (GCT→ACT) miaBisopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base T (2/3);  total (23/23)
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.79e-01

ATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGAT  >  AP009048/695145‑695283
                                                                      |                                                                    
aTAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGc                                                                      >  1:3519191/1‑71 (MQ=255)
aTAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGc                                                                      >  1:3196733/1‑71 (MQ=255)
 tAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGcc                                                                     <  1:2005041/71‑1 (MQ=255)
          cTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTcgcg                                                             >  1:1484552/1‑70 (MQ=255)
             gATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGAtt                                                         >  1:2586261/1‑71 (MQ=255)
             gATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGAtt                                                         <  1:780166/71‑1 (MQ=255)
              aTTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGAtt                                                         >  1:3204107/1‑70 (MQ=255)
               ttcttcTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGAtt                                                         >  1:53732/1‑69 (MQ=255)
                tcttctCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGc                                                                      >  1:2445414/1‑55 (MQ=255)
                    ctctTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGc                                                     <  1:3801353/68‑1 (MQ=255)
                    ctctTTTAACAGTTTCAAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAgg                                                  >  1:2918102/1‑71 (MQ=255)
                       ttttAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTcgcg                                                             <  1:418571/57‑1 (MQ=255)
                         ttAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGa                                                           >  1:4404984/1‑57 (MQ=255)
                         ttAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGa                                                           >  1:766588/1‑57 (MQ=255)
                          tAACAGTTTCCAGCGACCCAACTGATGGAATACTTTTTCCTGAGTCTtc                                                                  <  1:2949740/49‑1 (MQ=39)
                            aCAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcag                                           >  1:4360332/1‑70 (MQ=255)
                             cAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGgtg                                                <  1:1586059/64‑1 (MQ=255)
                              aGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTcagcag                                        <  1:824381/71‑1 (MQ=255)
                              aGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagaag                                        <  1:2962685/71‑1 (MQ=255)
                               gTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagcag                                        >  1:3852402/1‑70 (MQ=255)
                                tttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTtctc                                                                <  1:1796967/45‑1 (MQ=255)
                                tttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTtctc                                                                <  1:81443/45‑1 (MQ=255)
                                 ttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTcg                                                               <  1:916213/45‑1 (MQ=255)
                                      gCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAg                                                      >  1:1545592/1‑49 (MQ=255)
                                      gCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAg                                                   <  1:4003229/52‑1 (MQ=255)
                                           ccAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGct                           >  1:2736116/1‑71 (MQ=255)
                                             aaCTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGgtg                                                >  1:1016580/1‑48 (MQ=255)
                                              aCTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttat                        >  1:3351497/1‑69 (MQ=255)
                                               cTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGa                                                       >  1:2270465/1‑39 (MQ=255)
                                               cTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagcag                                        >  1:3205972/1‑54 (MQ=255)
                                               cTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttctt                       <  1:1297944/71‑1 (MQ=255)
                                                       aGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACg               >  1:4553470/1‑71 (MQ=255)
                                                        gACTTTTTCCTGAGTCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGt              >  1:1882285/1‑71 (MQ=255)
                                                        gACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGt              >  1:3415491/1‑71 (MQ=255)
                                                        gACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGt              <  1:437418/71‑1 (MQ=255)
                                                        gACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGt              <  1:4099318/71‑1 (MQ=255)
                                                         aCTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcag                                           <  1:134050/41‑1 (MQ=255)
                                                         aCTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTc             <  1:4427814/71‑1 (MQ=255)
                                                         aCTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTc             >  1:3330158/1‑71 (MQ=255)
                                                          cTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTc             <  1:3762728/70‑1 (MQ=255)
                                                           tttttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGc                            >  1:3123114/1‑54 (MQ=255)
                                                            ttttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGt          >  1:3729296/1‑71 (MQ=255)
                                                             tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTc                                             <  1:2506045/35‑1 (MQ=255)
                                                             tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTcc                     >  1:3648395/1‑59 (MQ=255)
                                                             tttCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCg            >  1:3090435/1‑68 (MQ=255)
                                                               tCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCa        <  1:154096/70‑1 (MQ=255)
                                                               tCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCa        <  1:145508/70‑1 (MQ=255)
                                                                  tGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATcc                              <  1:3521544/45‑1 (MQ=255)
                                                                    aGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGa   <  1:3909127/70‑1 (MQ=255)
                                                                    aGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGAt  <  1:1280768/71‑1 (MQ=255)
                                                                      |                                                                    
ATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGAT  >  AP009048/695145‑695283

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: