Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,781,358 +G coding (144/474 nt) trmL → tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,3541.G95.7% 79.2 / ‑2.9 23coding (140/474 nt)trmLtRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (14/8);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCTATC  >  NC_000913/3781215‑3781488
                                                                                                                                             |                                                                                                                                      
aTGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCcgcgc                                                                                                                                         <  2:337559/140‑1 (MQ=255)
                tttACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCa                                                                                                                         >  1:68144/1‑140 (MQ=255)
                                                      atcCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATcgcg                                                                                   >  2:122397/1‑140 (MQ=255)
                                                      atcCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATcgcg                                                                                   <  2:380470/140‑1 (MQ=255)
                                                                     aaTACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAg                                                                    >  2:194901/1‑140 (MQ=255)
                                                                                tCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAg                                                         >  1:58437/1‑140 (MQ=255)
                                                                                            catcGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGcc                                             <  1:228541/140‑1 (MQ=255)
                                                                                                 aaCCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGcc                                             <  2:295443/135‑1 (MQ=255)
                                                                                                 aaCCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGcc                                             >  1:295443/1‑135 (MQ=255)
                                                                                                  cccGATGGGTTTTGCCTGGGACGATAAGCGCCTGCGCCGCGGGGGGGCGGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTcacc                                       <  2:335872/139‑1 (MQ=255)
                                                                                                     gATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGATTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACg                                    <  1:194901/140‑1 (MQ=255)
                                                                                                      aTGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGa                                   >  1:65317/1‑140 (MQ=255)
                                                                                                               gCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAgg                               >  2:363787/1‑135 (MQ=255)
                                                                                                               gCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAgg                               <  1:363787/135‑1 (MQ=255)
                                                                                                               gCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTacac                          <  2:58437/140‑1 (MQ=255)
                                                                                                                    ggACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAgg                               >  2:126369/1‑130 (MQ=255)
                                                                                                                    ggACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAgg                               <  1:308243/130‑1 (MQ=255)
                                                                                                                    ggACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAgg                               <  1:126369/130‑1 (MQ=255)
                                                                                                                    ggACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAgg                               >  2:308243/1‑130 (MQ=255)
                                                                                                                    ggACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCt                     >  2:228215/1‑140 (MQ=255)
                                                                                                                    ggACGATAAGCGCCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCt                     >  1:258678/1‑140 (MQ=255)
                                                                                                                               gcCTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGt          <  2:258678/140‑1 (MQ=255)
                                                                                                                               cccTGCGCCGCGCGGGGGCTGGACTATCACGATTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGt          <  2:68144/139‑1 (MQ=255)
                                                                                                                                ccTGCGCCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTa         >  2:26192/1‑140 (MQ=255)
                                                                                                                                   gcgcCGCGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGc      <  1:122397/140‑1 (MQ=255)
                                                                                                                                    cgccgcGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCt     >  1:378852/1‑140 (MQ=255)
                                                                                                                                    cgccgcGCGGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCt     >  1:101553/1‑140 (MQ=255)
                                                                                                                                       cgcgcgGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCTATc  <  2:101553/140‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                      
ATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCTATC  >  NC_000913/3781215‑3781488

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: