Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,535,339 T→C F525S (TTC→TCC)  feoB → Fe(2+) transporter permease subunit FeoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,535,3390TC79.3% 48.7 / 12.8 29F525S (TTC→TCC) feoBFe(2+) transporter permease subunit FeoB
Reads supporting (aligned to +/- strand):  ref base T (5/1);  new base C (16/7);  total (21/8)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

ACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCG  >  NZ_CP009273/3535250‑3535420
                                                                                         |                                                                                 
acacGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCAtttt                                                                                   >  1:187247/1‑90 (MQ=255)
   cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCt                                                                                >  2:633333/1‑90 (MQ=255)
     ttAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGa                                                                              >  1:549176/1‑90 (MQ=255)
               gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCaa                                                                    >  1:252226/1‑90 (MQ=255)
                aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAAc                                                                   <  2:227575/90‑1 (MQ=255)
                  tATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAg                                                                 <  2:829820/90‑1 (MQ=255)
                         aCCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAACGCTTTCAACAGCTTCTCg                                                          >  2:615659/1‑90 (MQ=255)
                             ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                      >  2:182026/1‑90 (MQ=255)
                             ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                      >  2:544849/1‑90 (MQ=255)
                                  cGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGGGATCATCATCGCCAGCATTTCCCTGAGCGCTTTCAACACCTTCTCGCTGAGCggg                                                 >  1:228174/1‑90 (MQ=255)
                                            gCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgtcg                                       >  2:105646/1‑90 (MQ=255)
                                            gCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgtcg                                       >  2:283157/1‑90 (MQ=255)
                                             cTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGa                                      >  1:660711/1‑90 (MQ=255)
                                              ttCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGAt                                     <  2:807464/90‑1 (MQ=255)
                                               tCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                                    >  2:741621/1‑90 (MQ=255)
                                                      gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCaa                             >  2:496444/1‑90 (MQ=255)
                                                      gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCaa                             >  2:313528/1‑90 (MQ=255)
                                                                     gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc              >  1:167365/1‑90 (MQ=255)
                                                                     gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc              >  1:645880/1‑90 (MQ=255)
                                                                     gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc              >  1:347024/1‑90 (MQ=255)
                                                                     gATCATCATCGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc              >  1:90353/1‑90 (MQ=255)
                                                                      atcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg             <  1:182026/90‑1 (MQ=255)
                                                                       tcatcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt            >  2:776757/1‑90 (MQ=255)
                                                                       tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt            >  2:27059/1‑90 (MQ=255)
                                                                       tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg             >  2:658376/1‑89 (MQ=255)
                                                                       tcatcatcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg             <  1:658376/89‑1 (MQ=255)
                                                                            atcGTCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc       <  1:105646/90‑1 (MQ=255)
                                                                                tCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGcc   <  1:496444/90‑1 (MQ=255)
                                                                                tCAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGcc   <  2:420214/90‑1 (MQ=255)
                                                                                 cAGCATTTCCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCg  >  1:95383/1‑90 (MQ=255)
                                                                                         |                                                                                 
ACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCG  >  NZ_CP009273/3535250‑3535420

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