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,3150AC96.4% 85.2 / ‑3.6 28Y540S (TAC→TCC) flhAformate hydrogenlyase transcriptional activator FlhA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (14/13);  total (15/13)
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.76e-01

TTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACC  >  NZ_CP009273/2849228‑2849397
                                                                                       |                                                                                  
ttCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTccc                                                                                  <  2:1516179/90‑1 (MQ=255)
ttCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTccc                                                                                  <  1:1504461/90‑1 (MQ=255)
        gacgTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAAc                                                                          >  2:845295/1‑90 (MQ=255)
         acgTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACg                                                                         <  2:1580210/90‑1 (MQ=255)
         acgTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACg                                                                         >  2:1181634/1‑90 (MQ=255)
             gCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTAtt                                                                     >  2:46293/1‑90 (MQ=255)
             gCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTAtc                                                                     >  1:1442250/1‑89 (MQ=255)
              cGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTc                                                                    <  2:1386582/90‑1 (MQ=255)
                        ccgcGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCAcc                                                          >  1:1500879/1‑90 (MQ=255)
                        ccgcGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCAcc                                                          >  1:883095/1‑90 (MQ=255)
                          gcgACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTg                                                        >  2:398193/1‑90 (MQ=255)
                           cgACCACCCGCGTTCTGAAAAAAAGGGCCGCCGCCCGTGAGTTCCGTAGCGTTCTTTATTCCCGCCTGAACGTATTCCCGATTCACCTgc                                                       <  1:1323218/90‑1 (MQ=255)
                            gACTACCCGCGATCTAAAAAAAATGGTCCCCGCCCGTGAGTTCCGTAGCGATCTCTTTTCCCGCCTGAACGTATTCCCGATTCCCCTgcc                                                      <  1:845295/90‑1 (MQ=255)
                             ccTACCCGCGATCTGAAAAAAATGGTCGCCGCCCGTGAGTTCCGTAGCGATCTCTTTTCCCGCCTGAACGTATTCCCGATTCACCTgccg                                                     <  1:21155/89‑1 (MQ=255)
                                      gATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTAcgc                                            >  2:1521642/1‑90 (MQ=255)
                                      gATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTAcgc                                            >  1:424003/1‑90 (MQ=255)
                                       aTCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTccc                                                                                  >  2:895734/1‑51 (MQ=255)
                                       aTCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTccc                                                                                  <  1:895734/51‑1 (MQ=255)
                                          tGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGc                                        <  1:398193/90‑1 (MQ=255)
                                           aaaaaaaaTGGTCCCCGCCCGTGATTCCCGTAGCGATCTCTTTCCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCg                                       <  1:1022746/89‑1 (MQ=255)
                                               aaaaTGGTCGCCGCCCGTGATTCCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCg                                   <  1:27492/90‑1 (MQ=255)
                                                  aTGGTCGCCGCCCGTGATTCCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGaa                                <  1:292037/90‑1 (MQ=255)
                                                     gTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGat                             >  2:916034/1‑90 (MQ=255)
                                                     gTCGCCGACCGTGAGTCCCGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGat                             <  2:349542/90‑1 (MQ=255)
                                                                 gAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGctgctg                 >  1:1054185/1‑90 (MQ=255)
                                                                      ccGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCCCCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGaa            <  2:1374572/90‑1 (MQ=255)
                                                                      ccGTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGaa            <  1:1521642/90‑1 (MQ=255)
                                                                        gTAGCGATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAg          >  1:617965/1‑90 (MQ=255)
                                                                             gATCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCttt     >  1:711900/1‑90 (MQ=255)
                                                                              aTCTCTATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTa    >  1:165879/1‑90 (MQ=255)
                                                                                ctctATTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTAcc  >  2:1095832/1‑90 (MQ=255)
                                                                                    aTTCCCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAg                                         <  2:620212/47‑1 (MQ=255)
                                                                                       |                                                                                  
TTCAGACGGACGTGCGTCTAATCGCCGCGACTAACCGCGATCTGAAAAAAATGGTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCACTACGCGAGCGTCCGGAAGATATTCCGCTGCTGGCGAAAGCCTTTACC  >  NZ_CP009273/2849228‑2849397

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