Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 459,256 T→C D123D (GAT→GAC lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913459,2560TC98.2% 194.4 / ‑4.2 56D123D (GAT→GAClonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (28/27);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGT  >  NC_000913/459163‑459350
                                                                                             |                                                                                              
cGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCgg                                                                                          >  2:982278/1‑100 (MQ=255)
 gaggagTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCggg                                                                                         <  1:211484/100‑1 (MQ=255)
    gggTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAAc                                                                                      >  2:1324347/1‑100 (MQ=255)
     ggTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACa                                                                                     <  2:840128/100‑1 (MQ=255)
      gTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACGCTTTTCTGCGAAGGCGGAGTAGCTGGAGGCGCCGACCATTGACGAGCGGGAACAg                                                                                    <  2:1375396/100‑1 (MQ=255)
      gTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAg                                                                                    >  2:603207/1‑100 (MQ=255)
      gTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAg                                                                                    <  2:531680/100‑1 (MQ=255)
         aCAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGCGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGaa                                                                                 >  2:1388637/1‑100 (MQ=255)
          cAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAg                                                                                <  1:233817/100‑1 (MQ=255)
             cgcgcgcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGc                                                                             >  2:1126641/1‑100 (MQ=255)
             cgcgcgcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGc                                                                             >  2:1507596/1‑100 (MQ=255)
                 cgcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGt                                                                         >  2:313293/1‑100 (MQ=255)
                  gcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTg                                                                        >  1:737722/1‑100 (MQ=255)
                    gTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCg                                                                      <  1:798026/100‑1 (MQ=255)
                    gTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCg                                                                      <  1:882331/100‑1 (MQ=255)
                      aTTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTa                                                                    >  2:787157/1‑100 (MQ=255)
                       tttCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTAc                                                                   >  1:1536095/1‑100 (MQ=255)
                       tttCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTAc                                                                   <  2:1031606/100‑1 (MQ=255)
                       tttCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTAc                                                                   >  2:1288685/1‑100 (MQ=255)
                         tCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTg                                                                 <  1:817920/100‑1 (MQ=255)
                         tCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTg                                                                 <  1:982284/100‑1 (MQ=255)
                          cTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGc                                                                >  2:1328733/1‑100 (MQ=255)
                            gcgcTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCaa                                                              >  1:828863/1‑100 (MQ=255)
                              gcTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATc                                                            <  1:624599/100‑1 (MQ=255)
                               ctctctGACAATGGCGAACACTTTTCTGCGATGGCGGAGTATCTGTAGTCGCCGACCATTGGCGAGCGGGAACAGGATGTGCTGGTGCGTGCTGCAATCa                                                           >  2:1568043/1‑100 (MQ=255)
                                 ctctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGc                                                         >  1:731649/1‑100 (MQ=255)
                                  tctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGcc                                                        <  1:313315/100‑1 (MQ=255)
                                   ctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCa                                                       >  1:1470365/1‑100 (MQ=255)
                                    tGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAg                                                      <  2:49941/100‑1 (MQ=255)
                                           ggCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAg                                               <  1:270437/100‑1 (MQ=255)
                                            gCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAgg                                              >  1:712322/1‑100 (MQ=255)
                                             cGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGc                                             >  2:143218/1‑100 (MQ=255)
                                                acacTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTAc                                          <  1:105954/100‑1 (MQ=255)
                                                 cacTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACa                                         <  2:985422/100‑1 (MQ=255)
                                                 cacTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACa                                         <  2:1295279/100‑1 (MQ=255)
                                                   cTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATc                                       >  2:798176/1‑100 (MQ=255)
                                                   cTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATc                                       >  1:1292331/1‑100 (MQ=255)
                                                    ttttCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCa                                      >  1:885571/1‑100 (MQ=255)
                                                     tttCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAAACAGCCAGTTCGAAGGCTACATCaa                                     <  2:633244/100‑1 (MQ=255)
                                                        cTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCt                                  <  1:1290666/100‑1 (MQ=255)
                                                          gCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGa                                <  2:1473754/100‑1 (MQ=255)
                                                                       aTCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATccc                   <  1:187883/100‑1 (MQ=255)
                                                                       aTCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATccc                   <  1:275453/100‑1 (MQ=255)
                                                                        tCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTGGAAGGGTACATCAAGCTGAACAAAAGAAGccc                   >  2:1389356/1‑99 (MQ=255)
                                                                          tGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCccacc                >  2:170577/1‑100 (MQ=255)
                                                                           ggAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCccacca               <  1:769689/100‑1 (MQ=255)
                                                                            gAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAg              <  1:170577/100‑1 (MQ=255)
                                                                               tCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAg           <  1:313293/100‑1 (MQ=255)
                                                                                cGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGt          >  2:1089100/1‑100 (MQ=255)
                                                                                  ccGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGc        <  1:24031/100‑1 (MQ=255)
                                                                                  ccGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGc        >  1:664381/1‑100 (MQ=255)
                                                                                    gACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTg      <  1:1424995/100‑1 (MQ=255)
                                                                                    gACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTg      <  1:6882/100‑1 (MQ=255)
                                                                                      ccATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGAc    >  1:1194273/1‑100 (MQ=255)
                                                                                      caaTTGACGAGCGAGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGAc    >  1:1066425/3‑100 (MQ=255)
                                                                                        aTTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGt  >  1:284063/1‑100 (MQ=255)
                                                                                             |                                                                                              
CGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGT  >  NC_000913/459163‑459350

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: