Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 459,715 A→G 16.2% Y535C (TAC→TGC)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048459,7150AG16.2% 164.4 / 23.9 68Y535C (TAC→TGC) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (20/37);  new base G (4/7);  total (24/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACGTAATGCACTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCAAAGCGGTT  >  AP009048/459645‑459779
                                                                      |                                                                
gAACGTAATGCACTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTtacta                                                                  >  1:3378057/1‑71 (MQ=255)
   cGTAATGCACTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTacac                                                               <  1:2151903/71‑1 (MQ=255)
     tAATGCACTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTGCAccc                                                             <  1:3543018/71‑1 (MQ=255)
       aTGCACTGAAAAAAGGTTAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGt                                                           <  1:2901884/71‑1 (MQ=255)
       aTGCACTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGt                                                           <  1:3850853/71‑1 (MQ=255)
        tGCACTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGt                                                           <  1:162622/70‑1 (MQ=255)
          cACTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGt                                                           >  1:184060/1‑68 (MQ=255)
          cACTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAg                                                        >  1:416768/1‑71 (MQ=255)
            cTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAggc                                                      >  1:2232850/1‑71 (MQ=255)
             tGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAggcg                                                     >  1:2995417/1‑71 (MQ=255)
              gAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTGCACCCGTGAggcgg                                                    >  1:1181965/1‑71 (MQ=255)
               aaaaaaGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGa                                                         >  1:1321866/1‑65 (MQ=255)
               aaaaaaGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAggcggg                                                   >  1:1257522/1‑71 (MQ=255)
               aaaaaaGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAggcgg                                                    >  1:230276/1‑70 (MQ=255)
               aaaaaaGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAggc                                                      <  1:1550027/68‑1 (MQ=255)
                 aaaaGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAg                                                        >  1:625953/1‑64 (MQ=255)
                 aaaaGGTGAGCTGACCGTCGACGATAGCGCCATTATAGGCATTATTCGTTACTGCACCCGTGAggcgggcg                                                 >  1:3302761/1‑71 (MQ=255)
                  aaaGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGgcgt                                                <  1:2930484/71‑1 (MQ=255)
                  aaaGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGgcgt                                                <  1:2563264/71‑1 (MQ=255)
                   aaGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGgcgtg                                               <  1:2696006/71‑1 (MQ=255)
                      gTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGgcgtgcg                                             <  1:4091520/70‑1 (MQ=255)
                       tGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGgcgt                                                <  1:512834/66‑1 (MQ=255)
                       tGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGcgtgcgtg                                           <  1:3558923/71‑1 (MQ=255)
                       tGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGcgtgcgtg                                           <  1:519567/71‑1 (MQ=255)
                       tGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGcgtgcgtg                                           <  1:951240/71‑1 (MQ=255)
                        gAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTGCACCCGTGAGGCGGgcgtgcg                                             >  1:2589852/1‑68 (MQ=255)
                         aGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTgg                                          <  1:915444/70‑1 (MQ=255)
                         aGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTgg                                          <  1:871607/70‑1 (MQ=255)
                            tGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTGCACCCGTGAGGCGGGCGTGCGTGGTc                                        <  1:3867926/69‑1 (MQ=255)
                             gACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGa                                                         >  1:2555711/1‑51 (MQ=255)
                              aCCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTg                                      <  1:1903336/69‑1 (MQ=255)
                              aCCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTg                                      <  1:231596/69‑1 (MQ=255)
                                cGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGc                                  >  1:2169902/1‑71 (MQ=255)
                                  tCGACGATAGCGCCATTATCGGCATTATTCGTTACTGCACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGt                                <  1:3040441/71‑1 (MQ=255)
                                  tCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGgcgt                                                <  1:3665111/55‑1 (MQ=255)
                                  tCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGgcgt                                                <  1:3659833/55‑1 (MQ=255)
                                  tCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGt                                >  1:3181184/1‑71 (MQ=255)
                                  tCGACGATAGAGCCATTATCGGCATTATTCGTTACTACACCCGTGAggcgg                                                    <  1:4066487/51‑1 (MQ=255)
                                   cgacgaTAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGt                                >  1:4164242/1‑70 (MQ=255)
                                    gacgaTAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGt                                <  1:4086539/69‑1 (MQ=255)
                                    gacgaTAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTg                               <  1:2021657/70‑1 (MQ=255)
                                    gacgaTAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTg                               >  1:1573961/1‑70 (MQ=255)
                                     acgaTAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGa                                                         <  1:2224603/43‑1 (MQ=255)
                                        aTAGCGCCATTATCGGCATTATTCGTTACTGCACCCGTGAggc                                                      <  1:3599815/43‑1 (MQ=255)
                                        aTAGCGCCATTATCGGCATTATTCGTTACTGCACCCGTGAGGCGGGCGTGCGTgg                                          <  1:559475/55‑1 (MQ=255)
                                         tAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAggcggg                                                   >  1:2024956/1‑45 (MQ=255)
                                         tAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAAtct                         <  1:2559445/71‑1 (MQ=255)
                                         tAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAAtc                          >  1:3449807/1‑70 (MQ=255)
                                             gcCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCaa                     <  1:3228037/71‑1 (MQ=255)
                                             gcCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCaa                     <  1:3155616/71‑1 (MQ=255)
                                             gcCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCaa                     <  1:4108864/71‑1 (MQ=255)
                                              ccATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCaaa                    <  1:182611/71‑1 (MQ=255)
                                                 ttATCGGCATTATTCGTTACTACACCCGTGAGGCGGgcgtgcg                                             <  1:1805428/43‑1 (MQ=255)
                                                 ttATCGGCATTATTCGTTACTACACCCGTGAGGCGGgcgtgcg                                             <  1:840665/43‑1 (MQ=255)
                                                 ttATCGGCATTATTCGTTACTACACCAGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACtg                 >  1:1943612/1‑71 (MQ=255)
                                                  tATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACtg                 >  1:1196987/1‑70 (MQ=255)
                                                   aTCGGCATTATTCGTTACTACACCCGTGAGGCGGGcgtgcgtg                                           >  1:1027983/1‑43 (MQ=255)
                                                   aTCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACtgt                <  1:662984/70‑1 (MQ=255)
                                                   aTCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACtgt                <  1:637412/70‑1 (MQ=255)
                                                    tCGGCATTATTCGTTACTGCACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACtgtg               <  1:3326083/70‑1 (MQ=255)
                                                          ttattCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCa          >  1:2204248/1‑69 (MQ=255)
                                                            attCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCAAAg       <  1:1752294/70‑1 (MQ=255)
                                                            attCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCAAAg       <  1:435235/70‑1 (MQ=255)
                                                            attCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCAAAg       <  1:1730280/70‑1 (MQ=255)
                                                             ttCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCAAAGCg     <  1:1745701/71‑1 (MQ=255)
                                                              tCGTTACTGCACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCaaa        >  1:1155994/1‑67 (MQ=255)
                                                               cGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCaaa                    >  1:1937867/1‑54 (MQ=255)
                                                                gTTACTGCACCCGTGATGCGGGCGTGCGTGGTCTGGAGCGTg                               <  1:2834008/42‑1 (MQ=38)
                                                                gTTACTACACCTGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCAAAGCGGtt  <  1:464731/71‑1 (MQ=255)
                                                                gTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTg                               <  1:2847158/42‑1 (MQ=255)
                                                                  tactacACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCaaa                    >  1:4393497/1‑51 (MQ=255)
                                                                  tactacACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCAAAGCGGtt  >  1:2621926/1‑69 (MQ=255)
                                                                   actacACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGaaa                            >  1:1889349/1‑42 (MQ=255)
                                                                   actacACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCAAAGCGGtt  >  1:2843452/1‑68 (MQ=255)
                                                                     tacACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGAGTCGCaaa        <  1:476736/60‑1 (MQ=255)
                                                                      |                                                                
GAACGTAATGCACTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGTGCGTGGTCTGGAGCGTGAAATCTCCAAACTGTGTCGCAAAGCGGTT  >  AP009048/459645‑459779

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