Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,785,687 G→C G75G (GGC→GGG kbl ← glycine C‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,785,6870GC96.5% 184.0 / ‑1.1 58G75G (GGC→GGGkblglycine C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (34/22);  total (35/23)
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.80e-01

ATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAA  >  NZ_CP009273/3785600‑3785772
                                                                                       |                                                                                     
aTCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATccc                                                                                     <  1:237216/90‑1 (MQ=255)
    tCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGaaa                                                                                 <  2:271528/90‑1 (MQ=255)
      cATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAcc                                                                               <  2:312031/90‑1 (MQ=255)
            cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                         >  1:104790/1‑90 (MQ=255)
            cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                         >  2:351774/1‑90 (MQ=255)
             aGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgag                                                                          <  2:1514676/88‑1 (MQ=255)
             aGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgag                                                                          >  1:1514676/1‑88 (MQ=255)
              ggAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAAt                                                                       >  1:610263/1‑90 (MQ=255)
               gAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATc                                                                      <  1:887667/90‑1 (MQ=255)
                   gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  1:503734/1‑90 (MQ=255)
                   gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  2:1413590/1‑90 (MQ=255)
                   gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  2:1559861/1‑90 (MQ=255)
                   gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  1:1259543/1‑90 (MQ=255)
                    ccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTc                                                                 >  2:582333/1‑90 (MQ=255)
                    ccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTc                                                                 >  2:104111/1‑90 (MQ=255)
                      gccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTccc                                                               >  1:300620/1‑90 (MQ=255)
                              gtGTTAAAGCTCTTTATGGCGGCCCTGGGTCCCGCAAAAAAACCGCCCCGAACCCACCCCGAACCCGTGAGAACCCATTCCCGCCTTTgc                                                       <  1:940904/89‑1 (MQ=255)
                               ttttCAAAGTCTTTATGGCTGTCCTGAGTGGCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTgcc                                                      <  2:1519358/88‑1 (MQ=255)
                                   cAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGccgccg                                                  <  1:1413590/90‑1 (MQ=255)
                                       ctctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAat                                              <  2:22435/90‑1 (MQ=255)
                                        tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  2:542739/90‑1 (MQ=255)
                                        tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  2:522062/90‑1 (MQ=255)
                                        tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  1:331784/90‑1 (MQ=255)
                                         ctTTATGGCTGTCCTGAGGGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcc                                            >  2:940904/1‑89 (MQ=255)
                                           ttATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcaga                                          <  1:578263/90‑1 (MQ=255)
                                             aTGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcagatc                                        <  1:900241/90‑1 (MQ=255)
                                             aTGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcagatc                                        <  1:72816/90‑1 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:693726/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGGGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:771640/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:1443894/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:1344437/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:8975/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:183336/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1225138/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:125159/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1457007/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1470409/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:201369/1‑90 (MQ=255)
                                                   gTCCTGAGCGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:100873/1‑90 (MQ=255)
                                                     ccTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatgat                                >  1:975922/1‑90 (MQ=255)
                                                          gTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGcc                           <  2:300620/90‑1 (MQ=255)
                                                          gTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGcc                           >  1:466212/1‑90 (MQ=255)
                                                          gTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGcc                           >  2:641283/1‑90 (MQ=255)
                                                             ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  2:1147856/1‑90 (MQ=255)
                                                             ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:883425/1‑90 (MQ=255)
                                                             ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:395842/1‑90 (MQ=255)
                                                               gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  1:784740/1‑90 (MQ=255)
                                                               gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  1:486104/90‑1 (MQ=255)
                                                                  aaTAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCgag                   >  2:182418/1‑90 (MQ=255)
                                                                    tAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  1:404673/1‑85 (MQ=255)
                                                                    tAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  2:404673/85‑1 (MQ=255)
                                                                    tAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAt                 <  2:975922/90‑1 (MQ=255)
                                                                      aaccGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAg               <  2:1099546/87‑1 (MQ=255)
                                                                       aaCGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgt              >  2:386060/1‑90 (MQ=255)
                                                                           cACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgtt          <  2:784740/90‑1 (MQ=255)
                                                                           cACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgtt          >  1:1369669/1‑90 (MQ=255)
                                                                           cACCGAAGCCATCCCGAAACCGTGAAATTCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgtt          <  1:1487604/90‑1 (MQ=255)
                                                                                  gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  2:435266/90‑1 (MQ=255)
                                                                                   ccATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaaa  <  2:1259543/90‑1 (MQ=255)
                                                                                       |                                                                                     
ATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAA  >  NZ_CP009273/3785600‑3785772

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