Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 454,444 A→G E34G (GAA→GGA)  lon → endopeptidase La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273454,4440AG100.0% 131.9 / NA 37E34G (GAA→GGA) lonendopeptidase La
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/16);  total (21/16)

GTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGAAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGCTTCAACGGATGA  >  NZ_CP009273/454357‑454528
                                                                                       |                                                                                    
gTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGgaa                                                                                    <  1:539974/90‑1 (MQ=255)
 ttCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCTGGCGGGgaaa                                                                                   >  1:176470/1‑90 (MQ=255)
         gCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATcc                                                                           <  2:141001/90‑1 (MQ=255)
          cATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCg                                                                          >  1:1009697/1‑90 (MQ=255)
           aTTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGt                                                                         <  2:567307/90‑1 (MQ=255)
                    cccGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGaa                                                                <  1:252436/90‑1 (MQ=255)
                      cGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAgc                                                              <  2:510991/90‑1 (MQ=255)
                      cGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAgc                                                              <  1:450219/90‑1 (MQ=255)
                            gCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGaa                                                                <  1:555398/82‑1 (MQ=255)
                            gCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTATGGaa                                                                >  2:555398/1‑82 (MQ=255)
                              cGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATgg                                                      >  2:18834/1‑90 (MQ=255)
                                cTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGAc                                                    <  2:282461/90‑1 (MQ=255)
                                 tGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGAcc                                                   <  1:202616/90‑1 (MQ=255)
                                 tGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGAcc                                                   >  1:785811/1‑90 (MQ=255)
                                 tGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGAcc                                                   >  2:888573/1‑90 (MQ=255)
                                  gcgcGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCTTTGTCTGGAAGCGGCGATGGACCa                                                  >  2:740947/1‑90 (MQ=255)
                                                 ttATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGgaa                                                                                    <  1:707788/41‑1 (MQ=255)
                                                 ttATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGgaa                                                                                    >  2:707788/1‑41 (MQ=255)
                                                 ttATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAgc                                                              >  2:999840/1‑63 (MQ=255)
                                                 ttATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAgc                                                              <  1:999840/63‑1 (MQ=255)
                                                 ttATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGAt                                                        >  1:989487/1‑69 (MQ=255)
                                                 ttATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGAt                                                        <  2:989487/69‑1 (MQ=255)
                                                 ttATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATaaa                                           <  1:617982/82‑1 (MQ=255)
                                                 ttATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATaaa                                           >  2:617982/1‑82 (MQ=255)
                                                 ttATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTAt                                   >  1:640050/1‑90 (MQ=255)
                                                   aTCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGt                                                                         >  2:910248/1‑50 (MQ=255)
                                                   aTCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGt                                                                         <  1:910248/50‑1 (MQ=255)
                                                                 aTCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCagaaa                   >  1:506519/1‑90 (MQ=255)
                                                                  tCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTcg                           <  1:157400/81‑1 (MQ=255)
                                                                  tCCCCTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTcg                           >  2:157400/1‑81 (MQ=255)
                                                                      cTTATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGc              >  2:305891/1‑90 (MQ=255)
                                                                       ttATTTGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGCt             <  1:948788/90‑1 (MQ=255)
                                                                            tGTCGGTCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGCTTCAAc        >  1:984842/1‑90 (MQ=255)
                                                                            tGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGCTTCAAc        >  2:356041/1‑90 (MQ=255)
                                                                            tGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGCTTCAAc        >  1:404481/1‑90 (MQ=255)
                                                                            tGTCGGGCGGGGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGAGCAGAAAGAAGCTTCAAc        <  2:640050/90‑1 (MQ=255)
                                                                              tCGGGCGGGGAAAAACTATCCGTTGTCTGGAAGCGGCGTTGGACCATGATAAAAAAATTATGCTGGCCGCGCAGAAAGAAGCTTCAACgg      >  2:833537/1‑90 (MQ=255)
                                                                               cgggcgggGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGCTTCAACGGa     >  2:859888/1‑90 (MQ=255)
                                                                                gggcgggGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGCTTCAGCGGAt    >  2:1067166/1‑90 (MQ=255)
                                                                                gggcgggGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGCTTCAACGGAt    <  2:619032/90‑1 (MQ=255)
                                                                                  gcgggGAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGCTTCAACGGATGa  >  1:347372/1‑90 (MQ=255)
                                                                                       |                                                                                    
GTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTGGTTTATCCGCACATGGTCATCCCCTTATTTGTCGGGCGGGAAAAATCTATCCGTTGTCTGGAAGCGGCGATGGACCATGATAAAAAAATTATGCTGGTCGCGCAGAAAGAAGCTTCAACGGATGA  >  NZ_CP009273/454357‑454528

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