Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,077,887 G→A K63K (AAG→AAA efeO → iron uptake system protein EfeO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,077,8870GA94.1% 179.9 / ‑0.6 51K63K (AAG→AAAefeOiron uptake system protein EfeO
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (20/28);  total (21/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 = 1.00e+00

CTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTC  >  NZ_CP009273/1077757‑1078011
                                                                                                                                  |                                                                                                                            
ctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTg                                                                                                                        <  1:1139707/137‑1 (MQ=255)
ctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTgg                                                                                                                       <  1:79406/138‑1 (MQ=255)
ctGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTgg                                                                                                                       <  1:1227733/138‑1 (MQ=255)
   cTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTACCGACGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGATGGAGTg                                                                                                                   >  1:1150876/1‑139 (MQ=255)
   cTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGa                                                                                                                      <  1:696141/136‑1 (MQ=255)
   cTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTg                                                                                                                   >  1:970085/1‑139 (MQ=255)
        aTGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGgag                                                                                                               <  2:374206/138‑1 (MQ=255)
        aTGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGgag                                                                                                               <  1:638390/138‑1 (MQ=255)
               aCGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTc                                                                                                         <  1:505295/137‑1 (MQ=255)
               aCGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCaa                                                                                                       <  1:1142521/139‑1 (MQ=255)
                    gcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCg                                                                                                  <  1:461660/139‑1 (MQ=255)
                        aTGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAt                                                                                              <  2:575252/139‑1 (MQ=255)
                        aTGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAt                                                                                              <  1:1065322/139‑1 (MQ=255)
                           tgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGgag                                                                                                               >  1:1150840/1‑119 (MQ=255)
                           tgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGgag                                                                                                               <  2:1150840/119‑1 (MQ=255)
                             ccGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAtggtgg                                                                                         >  2:108205/1‑139 (MQ=255)
                                  ggTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTg                                                                                                >  2:57303/1‑127 (MQ=255)
                                  ggTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTg                                                                                                <  1:57303/127‑1 (MQ=255)
                                  ggTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAtggtggtgg                                                                                      >  1:83905/1‑137 (MQ=255)
                                      aaaGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAagag                                                                                 <  2:699266/138‑1 (MQ=255)
                                      aaaGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAagag                                                                                 <  2:304558/138‑1 (MQ=255)
                                             ccGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACTCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGaaaa                                                                         >  2:855519/1‑139 (MQ=255)
                                             ccGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCg                                                                                                  <  2:769367/114‑1 (MQ=255)
                                             ccGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCg                                                                                                  >  1:769367/1‑114 (MQ=255)
                                                                       aTGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaaaa                                               >  1:192923/1‑139 (MQ=255)
                                                                        tGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAAt                                              >  2:683081/1‑139 (MQ=255)
                                                                        tGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAAt                                              >  1:549771/1‑139 (MQ=255)
                                                                                      aaCGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAatat                                                                      >  2:474141/1‑101 (MQ=255)
                                                                                      aaCGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAatat                                                                      <  1:474141/101‑1 (MQ=255)
                                                                                      aaCGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTa                                 >  1:1702/1‑138 (MQ=255)
                                                                                         gCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGc                             <  1:1289933/139‑1 (MQ=255)
                                                                                         gCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGc                             <  2:192923/139‑1 (MQ=255)
                                                                                                  acacAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGgag                                                                                                               >  1:710134/1‑48 (MQ=255)
                                                                                                  acacAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGgag                                                                                                               <  2:710134/48‑1 (MQ=255)
                                                                                                  acacAGTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAAt                    <  2:549771/139‑1 (MQ=255)
                                                                                                       gTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGat               <  1:855519/139‑1 (MQ=255)
                                                                                                       gTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGat               <  2:1115625/139‑1 (MQ=255)
                                                                                                       gTTCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                >  1:58592/1‑138 (MQ=255)
                                                                                                        ttCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAt                                                                                              >  1:536124/1‑59 (MQ=255)
                                                                                                        ttCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAt                                                                                              <  2:536124/59‑1 (MQ=255)
                                                                                                        ttCATTATTCAGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGatat             <  1:701213/139‑1 (MQ=255)
                                                                                                                  aGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGc                                                      >  1:1162271/1‑89 (MQ=255)
                                                                                                                  aGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGc                                                      <  2:1162271/89‑1 (MQ=255)
                                                                                                                  aGAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCg     <  2:349315/138‑1 (MQ=255)
                                                                                                                   gAACCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGt   <  2:1251151/139‑1 (MQ=255)
                                                                                                                    aaCCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGcccc                                                                >  2:633989/1‑77 (MQ=255)
                                                                                                                    aaCCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGcccc                                                                <  1:633989/77‑1 (MQ=255)
                                                                                                                    aaCCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGc                                        >  1:1189046/1‑101 (MQ=255)
                                                                                                                    aaCCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGc                                        <  2:1189046/101‑1 (MQ=255)
                                                                                                                    aaCCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTc  >  2:1140584/1‑139 (MQ=255)
                                                                                                                    aaCCACAGCCAGAAAGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTc  >  1:66785/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                            
CTGCTTTTATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTC  >  NZ_CP009273/1077757‑1078011

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