Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,223,232 (A)8→7 coding (63/294 nt) BW25113_RS06140 → YcgL domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,223,2250A.97.9% 198.3 / ‑1.3 47coding (56/294 nt)BW25113_RS06140YcgL domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (24/22);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

TATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGG  >  NZ_CP009273/1223136‑1223312
                                                                                         |                                                                                       
tATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGa                                                                                         >  2:467931/1‑90 (MQ=255)
  tGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGaaa                                                                                       <  1:1597849/90‑1 (MQ=255)
      aaaaCCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGaaaaaaa                                                                                   >  1:1519603/1‑90 (MQ=255)
            tGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGAt                                                                            >  1:1672545/1‑90 (MQ=255)
            tGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGAt                                                                            >  2:2410315/1‑90 (MQ=255)
            tGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGAt                                                                            >  2:1514560/1‑90 (MQ=255)
            tGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGAt                                                                            >  1:2304768/1‑90 (MQ=255)
             ggCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGAtt                                                                           >  2:1053435/1‑90 (MQ=255)
             ggCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGAtt                                                                           >  1:2130689/1‑90 (MQ=255)
                     aaaGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCgt                                                                   <  1:2240421/90‑1 (MQ=255)
                     aaaGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCgt                                                                   <  2:1996179/90‑1 (MQ=255)
                     aaaGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCgt                                                                   <  2:1742853/90‑1 (MQ=255)
                            aaaaGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTg                                                            >  1:2503292/1‑90 (MQ=255)
                                  aTGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAAc                                                      >  2:881439/1‑90 (MQ=255)
                                  aTGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAAc                                                      <  2:1377194/90‑1 (MQ=255)
                                  aTGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAAc                                                      <  2:2503292/90‑1 (MQ=255)
                                       ttGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACtgatg                                                 >  1:2145730/1‑90 (MQ=255)
                                         gtgtgATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGaa                                               >  2:1098717/1‑90 (MQ=255)
                                         gtgtgATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGaa                                               >  1:2296161/1‑90 (MQ=255)
                                         gtgtgATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGaa                                               >  1:953730/1‑90 (MQ=255)
                                         gtgtgATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGaa                                               >  1:693131/1‑90 (MQ=255)
                                          tgtgATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGaaa                                              <  1:2310941/90‑1 (MQ=255)
                                          tgtgATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGaaa                                              <  1:2402751/90‑1 (MQ=255)
                                              aTTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAGAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGtt                                          >  2:643391/1‑90 (MQ=255)
                                                    cGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTc                                    <  2:2296161/90‑1 (MQ=255)
                                                    cGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTc                                    <  1:525361/90‑1 (MQ=255)
                                                    cGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTc                                    <  2:2145730/90‑1 (MQ=255)
                                                    cGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTc                                    <  2:2304768/90‑1 (MQ=255)
                                                    cGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTc                                    <  1:467931/90‑1 (MQ=255)
                                                    cGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTc                                    <  1:2271367/90‑1 (MQ=255)
                                                       agcagcAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTcg                                                                      >  2:1785582/1‑53 (MQ=255)
                                                       agcagcAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTcg                                                                      <  1:1785582/53‑1 (MQ=255)
                                                        gcagcaAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGcc                                <  2:2009812/90‑1 (MQ=255)
                                                        gcagcaAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGcc                                <  2:953730/90‑1 (MQ=255)
                                                         cagcaAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCt                               <  2:977922/90‑1 (MQ=255)
                                                           gcaAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCa                             <  2:198607/90‑1 (MQ=255)
                                                               gCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTa                         <  1:2217190/90‑1 (MQ=255)
                                                                   gACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTa                         >  2:2485058/1‑86 (MQ=255)
                                                                   gACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTa                         <  1:2485058/86‑1 (MQ=255)
                                                                        gACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGAtt                >  2:2221238/1‑90 (MQ=255)
                                                                        gACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGAtt                >  1:2199490/1‑90 (MQ=255)
                                                                         aCCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTc               <  1:2642701/90‑1 (MQ=255)
                                                                                  tatGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTgg      >  1:1197039/1‑90 (MQ=255)
                                                                                  tatGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTgg      >  2:2315666/1‑90 (MQ=255)
                                                                                  tatGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTgg      >  1:2077872/1‑90 (MQ=255)
                                                                                  tatGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTgg      >  2:342919/1‑90 (MQ=255)
                                                                                    tGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGAt    >  1:749149/1‑90 (MQ=255)
                                                                                    tGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGAt    <  2:693131/90‑1 (MQ=255)
                                                                                    tGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGAt    >  1:2336256/1‑90 (MQ=255)
                                                                                       cGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATgg  >  2:2652592/1‑90 (MQ=255)
                                                                                         |                                                                                       
TATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGG  >  NZ_CP009273/1223136‑1223312

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