Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,421,635 C→G A56P (GCA→CCA)  ytfE ← iron‑sulfur cluster repair protein YtfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,421,6350CG93.0% 133.6 / 0.3 43A56P (GCA→CCA) ytfEiron‑sulfur cluster repair protein YtfE
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base G (12/28);  total (13/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

TGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGCCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCC  >  NZ_CP009273/4421550‑4421720
                                                                                     |                                                                                     
tGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGt                                                                                   >  1:452171/1‑90 (MQ=255)
 ggTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGtt                                                                                  <  2:873726/90‑1 (MQ=255)
    aGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAg                                                                               <  1:481223/90‑1 (MQ=255)
     gcgcACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGc                                                                              >  2:57947/1‑90 (MQ=255)
            atgatATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATg                                                                       >  2:356871/1‑90 (MQ=255)
             tgatATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTATTCAGCGAGCTTTGGCAGTTCAGCTTCAATGa                                                                      <  1:462177/90‑1 (MQ=255)
              gatATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGAc                                                                     <  1:619522/90‑1 (MQ=255)
                   ggTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCaa                                                                <  1:356871/90‑1 (MQ=255)
                   ggTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCaa                                                                <  1:38074/90‑1 (MQ=255)
                   ggTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGCCAGTTCAGCTTCAATGACCTCaa                                                                <  1:439910/90‑1 (MQ=255)
                    gTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAAc                                                               >  1:195419/1‑90 (MQ=255)
                      cGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATg                                                                       >  2:427009/1‑80 (MQ=255)
                      cGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATg                                                                       <  1:427009/80‑1 (MQ=255)
                        atgatTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATcc                                                           <  2:147718/90‑1 (MQ=255)
                        atgatTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATcc                                                           <  2:240008/90‑1 (MQ=255)
                         tgatTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCa                                                          <  2:452171/90‑1 (MQ=255)
                            tttCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGtt                                                       <  1:403857/90‑1 (MQ=255)
                             ttCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTc                                                      >  2:50623/1‑90 (MQ=255)
                             ttCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTc                                                      >  1:884685/1‑90 (MQ=255)
                                         ggCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTgccgc                                          <  2:464216/90‑1 (MQ=255)
                                           cgcTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGccgccg                                        <  1:728259/90‑1 (MQ=255)
                                             cTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCcgcg                                      >  2:354279/1‑90 (MQ=255)
                                                  ccAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCgcgcgcgc                                 <  1:50623/90‑1 (MQ=255)
                                                  ccAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTAAACATCCAGTTCTTTACGTGCCGCCgcgcgcgc                                 <  1:629507/90‑1 (MQ=255)
                                                    aGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCa                               >  2:294302/1‑90 (MQ=255)
                                                               aTCGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTa                    >  2:173497/1‑90 (MQ=255)
                                                                 cGGTTGTTCAGCGAGCTTTGGCAGTTCAGCTTCAATGCCCTCAACACCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTAcc                  <  2:160734/90‑1 (MQ=255)
                                                                    ttgttCAGCGAGCTTTGGCAGTTCAGCTTCAATGCCCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGcc               <  2:42241/90‑1 (MQ=255)
                                                                    ttgttCAGCGAGCTTTGGCAGTCCAGCTTCAATGCCCCCACCACCCAGTTCTTTACGTGCCCCCGCGCGCGCCAGCGTCTGCTTACCGcc               <  2:474378/90‑1 (MQ=255)
                                                                    ttgttCAGCGAGCTTCGCCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGcc               >  2:188201/1‑90 (MQ=255)
                                                                         cAGCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCATTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCAcagc          <  2:195419/90‑1 (MQ=255)
                                                                           gCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCAcagcag        <  1:16090/90‑1 (MQ=255)
                                                                           gCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCAcagcag        >  2:425373/1‑90 (MQ=255)
                                                                           gCGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGt       >  2:772132/1‑90 (MQ=255)
                                                                           gCGAGCTTTGGAAGTTCAGCTTCAATGCCCTCAACTTCCAGTTCTTTACGTGCCGCAGCGCGCGCCAGCGTCTGCTTACCGCCAcagcag        <  1:330143/90‑1 (MQ=255)
                                                                            cGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGt       <  1:808839/90‑1 (MQ=255)
                                                                            cGAGCTTTGGCAGTTCAGCTTCAATGACCTCAACACCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGt       <  2:705549/90‑1 (MQ=255)
                                                                             gAGCTTTGGCAGTTCAGCTTAAATGCCCTCAACACCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTa      <  2:425173/90‑1 (MQ=255)
                                                                              aGCTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTaa     <  1:564612/90‑1 (MQ=255)
                                                                                cTTTGGCAGTTCGGCTTCAATGACCTCAACATCCAGTTCTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATcc  <  1:772132/90‑1 (MQ=255)
                                                                                cTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATc   <  1:201820/90‑1 (MQ=255)
                                                                                cTTTGGCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATc   <  1:294302/90‑1 (MQ=255)
                                                                                cTTTGCCAGTTCAGCTCCACTGCCCTCACCACCCAGTTTTTTCCGTCCCGCCGCCCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATc   <  2:626772/90‑1 (MQ=255)
                                                                                     |                                                                                     
TGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGCCAGTTCAGCTTCAATGACCTCAACATCCAGTTCTTTACGTGCCGCCGCGCGCGCCAGCGTCTGCTTACCGCCACAGCAGTAATCC  >  NZ_CP009273/4421550‑4421720

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