Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,623,127 (T)5→4 intergenic (+287/‑24) yfgG → / → yfgH YfgG family protein/glycine zipper 2TM domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,623,1230T.92.3% 84.7 / 3.0 26intergenic (+283/‑28)yfgG/yfgHYfgG family protein/glycine zipper 2TM domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (13/11);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

TGCGTGAATTTTAATAAATTTAATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTC  >  NZ_CP009273/2623042‑2623205
                                                                                 |                                                                                  
tGCGTGAATTTTAATAAATTTAATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTaa                                                                           <  1:142411/90‑1 (MQ=255)
                          aTCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGaa                                                 >  1:294495/1‑90 (MQ=255)
                               tttATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTaa                                                                           >  1:622800/1‑59 (MQ=255)
                               tttATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTaa                                                                           <  2:622800/59‑1 (MQ=255)
                                 tataCGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTaaa                                          <  1:167386/90‑1 (MQ=255)
                                  ataCGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTaaaa                                         >  2:883975/1‑90 (MQ=255)
                                       cAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATg                                     >  2:615270/1‑90 (MQ=255)
                                       cAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGt                                    >  1:396746/1‑90 (MQ=255)
                                        gaTACATTTAGTGTCCTCTTTTGGTGATGTTAAATGTCTTATTTTCGTAATGGGTATAACAAGGAATATGGAGAAAATTTAAAAAGTGtc                                   >  1:733103/2‑90 (MQ=255)
                                         aTACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGtct                                  >  1:415470/1‑90 (MQ=255)
                                           aCATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGtcttc                                <  1:130747/90‑1 (MQ=255)
                                            cATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGtcttc                                <  1:615270/90‑1 (MQ=255)
                                                 aCTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATaa                                                                  >  1:759454/1‑50 (MQ=255)
                                                 aCTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATaa                                                                  <  2:759454/50‑1 (MQ=255)
                                                 aCTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCtg                          >  1:670708/1‑90 (MQ=255)
                                                   tttCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCtgtg                        >  2:229897/1‑90 (MQ=255)
                                                    ttCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTaa                                                                           >  1:915653/1‑38 (MQ=38)
                                                    ttCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTaa                                                                           <  2:915653/38‑1 (MQ=38)
                                                     tCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAAt                                                                          >  1:696878/1‑38 (MQ=38)
                                                     tCCTCTTTTGATGATCTTAAATGTCTTATTTTCGTAAt                                                                          <  2:696878/38‑1 (MQ=38)
                                                      ccTCTTTTGATGATCTTAAATGTCTTATTTTCGTAAtg                                                                         >  2:380195/1‑38 (MQ=38)
                                                      ccTCTTTTGATGATCTTAAATGTCTTATTTTCGTAAtg                                                                         <  1:380195/38‑1 (MQ=38)
                                                                   tCTTAAATGTCTTATTTTCGTAATGTGTATAACaa                                                               >  2:43134/1‑35 (MQ=255)
                                                                   tCTTAAATGTCTTATTTTCGTAATGTGTATAACaa                                                               <  1:43134/35‑1 (MQ=255)
                                                                     ttAAATGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTc      <  1:883975/90‑1 (MQ=255)
                                                                         aTGTCTTATTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCActc  <  2:586528/90‑1 (MQ=255)
                                                                                 |                                                                                  
TGCGTGAATTTTAATAAATTTAATCTATCCCTTTATACGCAATACATTTACTTTCCTCTTTTGATGATCTTAAATGTCTTATTTTTCGTAATGTGTATAACAAGGAATAGTGATGAAATTTAAAAAATGTCTTCTGCCTGTGGCAATGTTAGCGTCATTCACTC  >  NZ_CP009273/2623042‑2623205

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