Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,434,529 A→C F186V (TTT→GTT)  fabB ← beta‑ketoacyl‑ACP synthase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,434,5290AC79.3% 70.7 / 13.1 29F186V (TTT→GTT) fabBbeta‑ketoacyl‑ACP synthase I
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base C (9/14);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGAT  >  NZ_CP009273/2434440‑2434610
                                                                                         |                                                                                 
cGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCaaa                                                                                   <  1:135722/90‑1 (MQ=255)
     tgtCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGa                                                                              <  1:173934/90‑1 (MQ=255)
      gtCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGAt                                                                             >  2:201804/1‑90 (MQ=255)
       tCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATg                                                                            <  1:42546/90‑1 (MQ=255)
                        cAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTg                                                           >  2:652076/1‑90 (MQ=255)
                          gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGa                                                         >  2:98278/1‑90 (MQ=255)
                          gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGa                                                         >  2:202767/1‑90 (MQ=255)
                          gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGa                                                         <  1:812895/90‑1 (MQ=255)
                          gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGa                                                         >  1:464313/1‑90 (MQ=255)
                          gcgcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGa                                                         >  2:189259/1‑90 (MQ=255)
                            gcACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATc                                                       <  1:201804/90‑1 (MQ=255)
                                    tGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTAc                                               <  2:97763/90‑1 (MQ=255)
                                    tGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTAc                                               >  2:557735/1‑90 (MQ=255)
                                      cGTCGAATTCGCAAGCCATTTCCCAGCCCAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTg                                             <  2:532415/90‑1 (MQ=255)
                                      cGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTg                                             <  1:189259/90‑1 (MQ=255)
                                       gTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGc                                            >  2:729202/1‑90 (MQ=255)
                                          gAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                         <  1:557735/90‑1 (MQ=255)
                                          gAATTCGCAAGCCATTTCCCAGCACAGCTCTTCCCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGtt                                         <  2:614075/90‑1 (MQ=255)
                                                  aaGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATAc                                 <  2:774997/90‑1 (MQ=255)
                                                    gCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAg                               <  2:400290/90‑1 (MQ=255)
                                                         ttCCCAGCACAGCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGc                          <  1:331881/90‑1 (MQ=255)
                                                         ttCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGc                          <  1:729202/90‑1 (MQ=255)
                                                             cAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGaa                      >  2:37974/1‑90 (MQ=255)
                                                                    gCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGctct                                                 <  1:814811/56‑1 (MQ=255)
                                                                    gCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGctct                                                 >  2:814811/1‑56 (MQ=255)
                                                                    gCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTAc                                       <  2:64897/66‑1 (MQ=255)
                                                                    gCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTAc                                       >  1:64897/1‑66 (MQ=255)
                                                                    gCTCTTCGCCGCCGCCAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGcac               >  1:330449/1‑90 (MQ=255)
                                                                            ccgccgccAGCAACCACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGTACACGCGGAg       <  1:635531/90‑1 (MQ=255)
                                                                                 gccAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGAt  <  1:37974/90‑1 (MQ=255)
                                                                                         |                                                                                 
CGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTACCGATACAGTGTGCGGAAGTCGCACACGCGGAGCTGAT  >  NZ_CP009273/2434440‑2434610

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