Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,849,315 A→C Y540S (TAC→TCC)  flhA → formate hydrogenlyase transcriptional activator FlhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,849,3150AC93.8% 45.2 / ‑3.2 16Y540S (TAC→TCC) flhAformate hydrogenlyase transcriptional activator FlhA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (8/7);  total (9/7)
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.28e-01

GACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGA  >  NZ_CP009273/2849236‑2849386
                                                                               |                                                                       
gacgTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAAc                                                               >  2:44559/1‑90 (MQ=255)
 acgTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACg                                                              >  2:227534/1‑90 (MQ=255)
          tAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGa                                                     >  1:260600/1‑90 (MQ=255)
          tAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGa                                                     >  1:301517/1‑90 (MQ=255)
                 cgcgACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCt                                              >  2:261272/1‑90 (MQ=255)
                  gcgACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTg                                             >  1:216782/1‑90 (MQ=255)
                       tAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTgccgcc                                        <  2:296751/90‑1 (MQ=255)
                                   gAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCg                            <  2:216782/90‑1 (MQ=255)
                                       aaaaTGGCCCCCCCCCGTGAGTCCCGTAGCGATCTCTTTCCCCGCCTGAACGTTTTCCCGATTCCCCTGCCGCCACTACGCGAGCGTCCg                        <  1:531243/90‑1 (MQ=255)
                                          aTGGTCGCCGCCCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGaa                     <  2:66114/90‑1 (MQ=255)
                                             gTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGat                  >  1:173163/1‑90 (MQ=255)
                                               cgccgcCCGTGAGTTCCGTAGCGATCTTTTTTCCCGCCTGAACGTATTCCCGATTCCCCTGCCGCCACTACGCGAGCGTCCGGAAGatat                <  1:242820/90‑1 (MQ=37)
                                                    cccGTGATTCCCGTAGCGATCTTTTTTCCCGCCTGAACGTATTCCCGATTCCCCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCgc           <  1:61056/89‑1 (MQ=37)
                                                             tCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGa  >  2:329124/1‑90 (MQ=255)
                                                                     gATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCgctgct       >  1:195176/1‑77 (MQ=255)
                                                                     gATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCgctgct       <  2:195176/77‑1 (MQ=255)
                                                                               |                                                                       
GACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGA  >  NZ_CP009273/2849236‑2849386

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