Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 426,040 A→T T344S (ACG→TCG)  queA → S‑adenosylmethionine:tRNA ribosyltransferase‑isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913426,0400AT97.2% 118.0 / ‑3.5 36T344S (ACG→TCG) queAS‑adenosylmethionine:tRNA ribosyltransferase‑isomerase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (16/19);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

TCGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAG  >  NC_000913/425944‑426131
                                                                                                |                                                                                           
tCGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGt                                                                                          >  1:1152304/1‑100 (MQ=255)
 cGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTa                                                                                         >  1:397434/1‑100 (MQ=255)
        tGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCg                                                                                  <  2:1190821/100‑1 (MQ=255)
            ggTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAgg                                                                              <  2:952756/100‑1 (MQ=255)
              ttATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCa                                                                            >  2:93998/1‑100 (MQ=255)
                  cAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTa                                                                        <  2:656679/100‑1 (MQ=255)
                  cAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTa                                                                        >  2:930502/1‑100 (MQ=255)
                  cAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTa                                                                        >  2:929019/1‑100 (MQ=255)
                  cAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGCAATGCGATGTTTATCTCGTACAATCCGCAGGCAATTa                                                                        >  1:551509/1‑100 (MQ=255)
                   aaCACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTaa                                                                       >  2:723143/1‑100 (MQ=255)
                      acacCATGAGCGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGa                                                                    <  2:744463/100‑1 (MQ=255)
                          cATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTTCGTTGATGCGATTTTTATCTTGTACAATCCGCAGTAAATTAATGa                                                                    >  2:682272/1‑96 (MQ=255)
                           aTGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGcgcg                                                               <  1:668864/100‑1 (MQ=255)
                           aTGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGcgcg                                                               <  1:778871/100‑1 (MQ=255)
                                 gCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCgggg                                                         <  1:1000267/100‑1 (MQ=255)
                                    tataAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGt                                                      >  2:330945/1‑100 (MQ=255)
                                           cagcagTAGAAGAGAAATATCGATTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCcg                                               <  2:240588/100‑1 (MQ=255)
                                               ggTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCgcggc                                           <  1:753802/100‑1 (MQ=255)
                                               ggTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCgcggc                                           >  2:1085071/1‑100 (MQ=255)
                                                    aagaGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTgg                                      >  2:75749/1‑100 (MQ=255)
                                                     agagAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGt                                     >  2:1286910/1‑100 (MQ=255)
                                                     agagAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGt                                     >  1:586182/1‑100 (MQ=255)
                                                       agaAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGttt                                   <  2:420056/100‑1 (MQ=255)
                                                        gaAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTa                                  >  2:112340/1‑100 (MQ=255)
                                                          aaTATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTaaa                                <  1:723148/100‑1 (MQ=255)
                                                               cGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGtt                           <  1:404650/100‑1 (MQ=255)
                                                                gCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTg                          <  2:579412/100‑1 (MQ=255)
                                                                gCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTg                          <  1:94000/100‑1 (MQ=255)
                                                                gCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTg                          >  2:276969/1‑100 (MQ=255)
                                                                 cTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTgg                         <  1:776919/100‑1 (MQ=255)
                                                                 cTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTgg                         >  1:323202/1‑100 (MQ=255)
                                                                   tttttAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGAc                       <  1:118647/100‑1 (MQ=255)
                                                                           tACTGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTc               <  1:534589/100‑1 (MQ=255)
                                                                           tACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTc               <  1:1260388/100‑1 (MQ=255)
                                                                                    gCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAg       <  1:220100/99‑1 (MQ=255)
                                                                                        tGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGgag  >  2:168652/1‑100 (MQ=255)
                                                                                                |                                                                                           
TCGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAG  >  NC_000913/425944‑426131

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