Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 459,642 A→G E252G (GAA→GGA)  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,6420AG95.2% 197.7 / ‑1.2 63E252G (GAA→GGA) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (26/34);  total (27/36)
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.94e-01

AGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGA  >  NC_000913/459546‑459737
                                                                                                |                                                                                               
aGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGc                                                                                              <  1:138176/100‑1 (MQ=255)
 gATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGcc                                                                                             >  2:229213/1‑100 (MQ=255)
  aTGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGccc                                                                                            >  2:1126173/1‑100 (MQ=255)
   tGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGTTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCt                                                                                           >  1:1550323/1‑100 (MQ=255)
   tGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCt                                                                                           <  1:716605/100‑1 (MQ=255)
     gagaAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGa                                                                                         >  1:1014767/1‑100 (MQ=255)
      agaAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGaa                                                                                        >  2:1275863/1‑100 (MQ=255)
         aaTCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAgcg                                                                                     <  2:358694/100‑1 (MQ=255)
          aTCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAgcgc                                                                                    >  2:207631/1‑100 (MQ=255)
           tCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCa                                                                                   >  2:7113/1‑100 (MQ=255)
                gCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATc                                                                              <  1:732894/100‑1 (MQ=255)
                  gTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGa                                                                            <  2:1530800/100‑1 (MQ=255)
                    gAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGAcg                                                                          <  2:263140/100‑1 (MQ=255)
                     aGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGAcgc                                                                         >  1:1283650/1‑100 (MQ=255)
                      gTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGAcgcg                                                                        >  2:1286965/1‑100 (MQ=255)
                        aCTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACgcggc                                                                      >  2:770643/1‑100 (MQ=255)
                           aTCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACGCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGaa                                                                   >  1:1115421/1‑100 (MQ=255)
                            tCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAg                                                                  <  1:229213/100‑1 (MQ=255)
                             ntGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAATCGACGCGGCGAAGAt                                                                >  2:538474/2‑100 (MQ=255)
                             cTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGa                                                                 >  1:1530319/1‑100 (MQ=255)
                             cTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGa                                                                 >  1:510654/1‑100 (MQ=255)
                               gAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATg                                                               >  1:163561/1‑100 (MQ=255)
                               gAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATg                                                               >  1:214988/1‑100 (MQ=255)
                                 aCGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                             <  1:718853/100‑1 (MQ=255)
                                    aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                          <  1:537309/100‑1 (MQ=255)
                                    aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGa                                                           >  2:384387/1‑99 (MQ=255)
                                      cAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAag                                                        <  2:1014761/100‑1 (MQ=255)
                                       aaaTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGCGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAaga                                                       <  2:864579/100‑1 (MQ=255)
                                       aaaTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAaga                                                       <  1:1399849/100‑1 (MQ=255)
                                       aaaTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAaga                                                       >  1:75204/1‑100 (MQ=255)
                                         aTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAgg                                                     <  1:35073/100‑1 (MQ=255)
                                             aaGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCaaa                                                 <  2:214988/100‑1 (MQ=255)
                                              aGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCaaaa                                                >  1:107856/1‑100 (MQ=255)
                                                cTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAaga                                              <  1:1442983/100‑1 (MQ=255)
                                                 tATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAagag                                             <  2:1019685/100‑1 (MQ=255)
                                                     cAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAg                                         >  2:224348/1‑100 (MQ=255)
                                                         aagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGa                                     <  2:272757/100‑1 (MQ=255)
                                                          agaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGaa                                    <  1:1454983/100‑1 (MQ=255)
                                                              cTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCag                                <  1:1286970/100‑1 (MQ=255)
                                                                cGGTGAAATGGACGACGCGCCGGACGAAACCGAAGCCCTGAAGCGCAAAATCCTCGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCagag                              <  1:1401273/100‑1 (MQ=255)
                                                                cGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCagag                              >  2:1190449/1‑100 (MQ=255)
                                                                   tGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTg                           <  1:7113/100‑1 (MQ=255)
                                                                    gAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGACAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGc                          >  2:737323/1‑100 (MQ=255)
                                                                    gAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTgt                          <  1:211144/100‑2 (MQ=255)
                                                                    gAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGc                          <  2:70266/100‑1 (MQ=255)
                                                                     aaaTGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCa                         <  1:1487171/100‑1 (MQ=255)
                                                                          gacgacGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGc                    <  1:715309/100‑1 (MQ=255)
                                                                          gacgacGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGc                    <  1:1313616/100‑1 (MQ=255)
                                                                           acgacgCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCt                   <  1:1026040/100‑1 (MQ=255)
                                                                              acgCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGATGCAGAGTTGCAGAAGCTGaa                <  2:965096/100‑1 (MQ=255)
                                                                                 cgcCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGaa                <  2:997377/97‑1 (MQ=255)
                                                                                    cGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat          <  1:784429/100‑1 (MQ=255)
                                                                                     ggACGAAAACGGGGCCCTGAAGCGCAAAATCGACGCGGGGAAGATGCCGAAAGTGGCAAAAGAGAAAGCGGAGGCAGTGTTGCAGAAGCTGAAAatgatg         <  2:1254007/100‑1 (MQ=255)
                                                                                     ggACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgatg         >  1:457104/1‑100 (MQ=255)
                                                                                       aCGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtc       <  1:207631/100‑1 (MQ=255)
                                                                                        ngAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct      >  2:481963/2‑100 (MQ=255)
                                                                                        cGAAAACGGAGCCCTGAAGCGCAAACTCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct      >  1:1198409/1‑100 (MQ=255)
                                                                                        cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct      >  1:1414367/1‑100 (MQ=255)
                                                                                          aaaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc    <  2:1530310/100‑1 (MQ=255)
                                                                                          aaaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAAAGATGtct      >  2:1242307/1‑98 (MQ=255)
                                                                                           aaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCg   <  2:152828/100‑1 (MQ=255)
                                                                                           aaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCg   <  2:1283645/100‑1 (MQ=255)
                                                                                            aaCGGAGCCCTGAAGCGCAAAATCGAAGCGGCGAAGAGGGCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGa  <  1:1523518/100‑1 (MQ=255)
                                                                                                |                                                                                               
AGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGA  >  NC_000913/459546‑459737

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