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,2560TC92.7% 213.0 / 5.3 68D123D (GAT→GAClonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (25/38);  total (27/41)
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.63e-01

CTGGTCGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGTCG  >  NC_000913/459158‑459352
                                                                                                  |                                                                                                
cTGGTCGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATtgatg                                                                                                 >  1:442601/1‑100 (MQ=255)
  ggTCGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATTGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAg                                                                                               >  1:1078713/1‑100 (MQ=255)
  ggTCGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGTGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAg                                                                                               >  1:296582/1‑100 (MQ=255)
   gTCGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGc                                                                                              <  2:820988/100‑1 (MQ=255)
     cGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGGCGCCGACCATTGACGAGCgg                                                                                            >  1:1708824/1‑100 (MQ=255)
      gAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCggg                                                                                           <  1:287319/100‑1 (MQ=255)
      gAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCggg                                                                                           <  1:1254650/100‑1 (MQ=255)
       aGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGa                                                                                          <  2:886822/100‑1 (MQ=255)
       aGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGa                                                                                          <  2:8696/100‑1 (MQ=255)
           gTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAg                                                                                      >  1:131506/1‑100 (MQ=255)
            ttACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAgg                                                                                     <  1:1295791/100‑1 (MQ=255)
             tACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGTGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAgg                                                                                     >  1:1362604/1‑99 (MQ=255)
             tACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGa                                                                                    <  1:422536/100‑1 (MQ=255)
               cAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAg                                                                                  <  2:1512427/100‑1 (MQ=255)
                aGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGt                                                                                 <  1:1500647/100‑1 (MQ=255)
                  cgcgcgcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGc                                                                               >  1:451280/1‑100 (MQ=255)
                  cgcgcgcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGc                                                                               <  1:1260665/100‑1 (MQ=255)
                  cgcgcgcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGc                                                                               <  1:1006319/100‑1 (MQ=255)
                   gcgcgcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCt                                                                              <  2:1523150/100‑1 (MQ=255)
                      cgcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGt                                                                           >  1:1397211/1‑100 (MQ=255)
                       gcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTg                                                                          <  2:1650159/100‑1 (MQ=255)
                       gcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTg                                                                          <  2:1723075/100‑1 (MQ=255)
                       gcgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTg                                                                          <  2:920353/100‑1 (MQ=255)
                       gcgTATTTCTGCGCTCTCTGACAACGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTg                                                                          <  2:748276/100‑1 (MQ=255)
                        cgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGc                                                                         >  1:1739786/1‑100 (MQ=255)
                        cgTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGc                                                                         <  1:336454/100‑1 (MQ=255)
                         gTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCg                                                                        <  2:637096/100‑1 (MQ=255)
                           aTTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTa                                                                      >  1:712877/1‑100 (MQ=255)
                           aTTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTa                                                                      <  1:650558/100‑1 (MQ=255)
                             ttCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACt                                                                    >  1:707741/1‑100 (MQ=255)
                              tCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTg                                                                   >  1:434431/1‑100 (MQ=255)
                               cTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGc                                                                  <  1:1781972/100‑1 (MQ=255)
                               cTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGc                                                                  <  2:210223/100‑1 (MQ=255)
                                 gcgcTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCaa                                                                <  2:451282/100‑1 (MQ=255)
                                 gcgcTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCaa                                                                <  2:558169/100‑1 (MQ=255)
                                  cgcTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAAt                                                               >  2:1402039/1‑100 (MQ=255)
                                      ctctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGc                                                           <  2:1036406/100‑1 (MQ=255)
                                       tctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGcc                                                          <  1:908957/100‑1 (MQ=255)
                                       tctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGcc                                                          <  1:1419301/100‑1 (MQ=255)
                                       tctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGcc                                                          <  2:1328295/100‑1 (MQ=255)
                                        ctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCa                                                         >  2:563856/1‑100 (MQ=255)
                                        ctGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCa                                                         >  1:1000376/1‑100 (MQ=255)
                                            cAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTc                                                     >  1:1019032/1‑100 (MQ=255)
                                            cAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGAAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTc                                                     >  1:1040892/1‑100 (MQ=255)
                                                 gCGAACACTTTTCTGTGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAgg                                                <  2:296584/100‑1 (MQ=255)
                                                 gCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAgg                                                <  2:138504/100‑1 (MQ=255)
                                                   gAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCt                                              >  2:1360199/1‑100 (MQ=255)
                                                    aaCACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTa                                             <  1:1169864/100‑1 (MQ=255)
                                                        cTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATc                                         >  1:724192/1‑100 (MQ=255)
                                                             cTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCt                                    <  2:1820543/100‑1 (MQ=255)
                                                              tGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTg                                   >  2:263757/1‑100 (MQ=255)
                                                              tGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTg                                   >  1:96163/1‑100 (MQ=255)
                                                                      cGGAGTATCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACaaaaa                           <  1:988128/100‑1 (MQ=255)
                                                                            aTCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATccc                     <  2:274866/100‑1 (MQ=255)
                                                                            aTCTGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATccc                     <  2:724193/100‑1 (MQ=255)
                                                                               tGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCccacc                  <  2:1862380/100‑1 (MQ=255)
                                                                               tGGAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCccacc                  >  1:1021527/1‑100 (MQ=255)
                                                                                ggAGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCccacca                 >  1:1095290/1‑100 (MQ=255)
                                                                                  aGTCGCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGa               <  2:1040893/100‑1 (MQ=255)
                                                                                    tCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAg             <  2:592034/100‑1 (MQ=255)
                                                                                      gCCGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTg           >  1:75052/1‑100 (MQ=255)
                                                                                       ccGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGc          <  2:1276396/100‑1 (MQ=255)
                                                                                        cGACCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCt         >  2:1452043/1‑100 (MQ=255)
                                                                                          aCCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGa       <  2:1000377/100‑1 (MQ=255)
                                                                                          aCCATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGa       <  1:1402031/100‑1 (MQ=255)
                                                                                            cATTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACg     <  2:131507/100‑1 (MQ=255)
                                                                                             aTTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGt    >  2:355425/1‑100 (MQ=255)
                                                                                             aTTGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGt    >  2:823672/1‑100 (MQ=255)
                                                                                              ttGACGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGt    >  2:508501/1‑99 (MQ=255)
                                                                                               tgacgaGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGTCg  >  2:846030/1‑100 (MQ=255)
                                                                                                  |                                                                                                
CTGGTCGAGGGGTTACAGCGCGCGCGTATTTCTGCGCTCTCTGACAATGGCGAACACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGAAGTGCTGACGTCG  >  NC_000913/459158‑459352

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