Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 886,189 A→G G7G (GGT→GGC grxA ← glutaredoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273886,1890AG100.0% 274.9 / NA 74G7G (GGT→GGCgrxAglutaredoxin 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (42/32);  total (42/32)

ACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGT  >  NZ_CP009273/886103‑886274
                                                                                      |                                                                                     
aCATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCCTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCa                                                                                    <  1:2934396/90‑1 (MQ=255)
   tactgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCaaaa                                                                                 >  2:468587/1‑90 (MQ=255)
   tactgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCaaaa                                                                                 >  1:1764578/1‑90 (MQ=255)
     ctgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTATGGGCAACCCGAACGGCCa                                                                                    <  1:2639797/85‑1 (MQ=255)
     ctgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTATGGGCAACCCGAACGGCCa                                                                                    >  2:2639797/1‑85 (MQ=255)
     ctgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAAt                                                                               >  2:140922/1‑90 (MQ=255)
     ctgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAAt                                                                               >  1:2448670/1‑90 (MQ=255)
     ctgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAAt                                                                               >  2:2856381/1‑90 (MQ=255)
      tgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATa                                                                              >  2:243589/1‑90 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCa                                                                                    >  2:1280678/1‑83 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCa                                                                                    <  1:1280678/83‑1 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  2:1105167/1‑90 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  2:2687307/1‑90 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  2:2642728/1‑90 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  2:1338316/1‑90 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  1:431413/1‑90 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             <  1:2980610/90‑1 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             >  1:478395/1‑90 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                             <  1:1338316/90‑1 (MQ=255)
        atactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAAc                                                                            >  1:2713198/1‑90 (MQ=255)
        atactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAAc                                                                            >  2:2686182/1‑90 (MQ=255)
            tgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGtt                                                                        <  2:2557569/90‑1 (MQ=255)
             gaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGttt                                                                       >  2:2492642/1‑90 (MQ=255)
               aaaTCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGc                                                                     >  2:2038487/1‑90 (MQ=255)
                     tcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGttt                                                                       <  2:1547312/82‑1 (MQ=255)
                     tcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGttt                                                                       >  1:1547312/1‑82 (MQ=255)
                     tcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCat                                                                   >  2:1565613/1‑86 (MQ=255)
                     tcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCat                                                                   <  1:1565613/86‑1 (MQ=255)
                         gTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                           <  1:2686182/90‑1 (MQ=255)
                         gTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                           <  1:2492642/90‑1 (MQ=255)
                         gTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                           <  1:621365/90‑1 (MQ=255)
                         gTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                           <  1:207750/90‑1 (MQ=255)
                           tCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTtctctc                                                         <  1:2856381/90‑1 (MQ=255)
                                 cTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATa                                                   <  1:468587/90‑1 (MQ=255)
                                      tttCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATCATTTCTCTCCTCATAGAttt                                              <  1:2687307/90‑1 (MQ=255)
                                        tctcAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGtt                                                                        >  2:2989115/1‑62 (MQ=255)
                                        tctcAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGtt                                                                        <  1:2989115/62‑1 (MQ=255)
                                          taaGCCAGATCTTTTGCCGGTACACAGTAAGGGCACCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTATGcc                                         <  1:931026/88‑1 (MQ=255)
                                           cAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                         >  2:2107561/1‑90 (MQ=255)
                                           cAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                         <  1:519483/90‑1 (MQ=255)
                                           cAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                         <  1:340119/90‑1 (MQ=255)
                                                 gATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAAt                                   >  2:2126526/1‑90 (MQ=255)
                                                     ttttGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATc                               <  1:534777/90‑1 (MQ=255)
                                                      tttGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCa                              <  1:732102/90‑1 (MQ=255)
                                                           acgcacACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTa                         >  2:2059378/1‑90 (MQ=255)
                                                                       aGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGc             <  1:140922/90‑1 (MQ=255)
                                                                        gggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCt            >  2:1062779/1‑90 (MQ=255)
                                                                        gggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCt            <  2:674377/90‑1 (MQ=255)
                                                                        gggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGACCACGCTAAAATGTATTCGCt            >  2:902679/1‑90 (MQ=255)
                                                                         ggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTg           <  1:2488927/90‑1 (MQ=255)
                                                                         ggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTg           <  1:2642728/90‑1 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                           <  2:738217/39‑1 (MQ=39)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                           >  1:738217/1‑39 (MQ=39)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                           <  2:564034/39‑1 (MQ=39)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                           >  1:564034/1‑39 (MQ=39)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATa                                                   <  2:98695/47‑1 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATa                                                   >  1:98695/1‑47 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTAtt                >  1:1119580/1‑82 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTAtt                <  2:1119580/82‑1 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:2170909/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:974629/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:1185746/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:1652918/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:2042476/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:630847/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  1:63946/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:1627254/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:2098805/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa        >  2:2273397/1‑90 (MQ=255)
                                                                              cccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGt      <  1:243589/90‑1 (MQ=255)
                                                                              cccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGt      <  1:2209975/90‑1 (MQ=255)
                                                                                  aaCGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGAttt                                              <  1:629713/46‑1 (MQ=255)
                                                                                  aaCGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGAttt                                              >  2:629713/1‑46 (MQ=255)
                                                                                  aaCGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGt  >  2:818521/1‑90 (MQ=255)
                                                                                      |                                                                                     
ACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGT  >  NZ_CP009273/886103‑886274

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