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.89.3% 82.8 / 9.9 28coding (56/294 nt)BW25113_RS06140YcgL domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base . (16/9);  total (18/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGG  >  NZ_CP009273/1223136‑1223312
                                                                                         |                                                                                       
tATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGa                                                                                         >  2:1404599/1‑90 (MQ=255)
 aTGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGaa                                                                                        >  2:156288/1‑90 (MQ=255)
  tGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGaaa                                                                                       >  1:700737/1‑90 (MQ=255)
            tGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGaaaaa                                                                                     >  2:938252/1‑82 (MQ=255)
            tGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGaaaaa                                                                                     <  1:938252/82‑1 (MQ=255)
                  cTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTcg                                                                      >  2:311330/1‑90 (MQ=255)
                     aaaGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCgt                                                                   <  2:700737/90‑1 (MQ=255)
                     aaaGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCgt                                                                   <  2:432411/90‑1 (MQ=255)
                          gTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTcc                                                              >  2:567389/1‑90 (MQ=255)
                                  aTGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAAc                                                      >  2:78863/1‑90 (MQ=255)
                                  aTGTCTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAAc                                                      >  2:1294428/1‑90 (MQ=255)
                                    gTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGa                                                        >  2:691835/1‑86 (MQ=255)
                                    gTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGa                                                        <  1:691835/86‑1 (MQ=255)
                                    gTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACtg                                                    >  2:1038501/1‑90 (MQ=255)
                                       ttGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACtgatg                                                 >  1:1014231/1‑90 (MQ=255)
                                       ttGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACtgatg                                                 >  1:133522/1‑90 (MQ=255)
                                         gtgtgATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGaa                                               >  1:295801/1‑90 (MQ=255)
                                                ttATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGtt                                                                >  2:1039210/1‑66 (MQ=255)
                                                ttATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGtt                                                                <  1:1039210/66‑1 (MQ=255)
                                                  aTCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTgg                                      <  1:567389/90‑1 (MQ=255)
                                                    cGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTc                                    <  1:80997/90‑1 (MQ=255)
                                                    cGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTc                                    <  1:787370/90‑1 (MQ=255)
                                                                   gACCAGACCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCga                     >  2:1459562/1‑90 (MQ=255)
                                                                         aCCTATTTATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTc               >  2:65516/1‑90 (MQ=255)
                                                                            tATTTATATGTCGAAAAAAAGACGATTTTTCGCgt                                                                   >  2:1132160/1‑35 (MQ=38)
                                                                            tATTTATATGTCGAAAAAAAGACGATTTTTCGCgt                                                                   <  1:1132160/35‑1 (MQ=38)
                                                                               ttATATGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGc         <  1:1292220/90‑1 (MQ=255)
                                                                                  tatGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGcc                                <  1:792337/65‑1 (MQ=255)
                                                                                  tatGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGcc                                >  2:792337/1‑65 (MQ=255)
                                                                                  tatGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTgg      >  2:345015/1‑90 (MQ=255)
                                                                                  tatGTCGAAAAAAAGACGATTTTTCGCGTGTTCATGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTgg      >  1:575708/1‑90 (MQ=255)
                                                                                    tGTCGAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGAt    >  1:171059/1‑90 (MQ=255)
                                                                                       cGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATgg  >  1:503519/1‑90 (MQ=255)
                                                                                         |                                                                                       
TATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGG  >  NZ_CP009273/1223136‑1223312

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