Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,723,625 T→G 12.2% I190I (ATA→ATC)  ynzC ← FeoB family ferrous iron (Fe2+) uptake protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,723,6250TG12.2% 92.0 / 9.7 41I190I (ATA→ATC) ynzCFeoB family ferrous iron (Fe2+) uptake protein
Reads supporting (aligned to +/- strand):  ref base T (16/20);  new base G (2/3);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

ACAAGAACTACCAAAGCTTATTAAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACGAGCGAATATGTTCCAAGTGATAGCACAC  >  CP000730/2723487‑2723764
                                                                                                                                          |                                                                                                                                           
aCAAGAACTACCAAAGCTTATTAAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTAt                                                                                                                                           <  1:401242/141‑1 (MQ=255)
   aGAACTACCAAAGCTTATTAAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTCTCCa                                                                                                                                        <  2:70090/141‑1 (MQ=255)
   aGAACTACCAAAGCTTATTAAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCa                                                                                                                                        >  2:272110/1‑141 (MQ=255)
    gAACTACCAAAGCTTATTAAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCATACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCAt                                                                                                                                       <  1:266547/141‑1 (MQ=255)
        tACCAAAGCTTATTAAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAATCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCGATCCATAACa                                                                                                                                   >  2:272600/1‑141 (MQ=255)
        tACCAAAGCTTATTAAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAATCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCGATCCATAACa                                                                                                                                   >  2:272576/1‑141 (MQ=255)
              aGCTTATTAAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACg                                                                                                                             >  2:84326/1‑141 (MQ=255)
                 ttattaAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGtt                                                                                                                          >  1:55450/1‑141 (MQ=255)
                    ttaAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCaa                                                                                                                       >  2:201329/1‑141 (MQ=255)
                                 cTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATatcat                                                                                                          <  2:35627/141‑1 (MQ=255)
                                 cTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATatcat                                                                                                          <  1:448407/141‑1 (MQ=255)
                                   gtgtCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATcc                                                                                                        >  2:265903/1‑141 (MQ=255)
                                   gtgtCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATcc                                                                                                        <  1:45163/141‑1 (MQ=255)
                                     gtCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTAACGTCTCAATTGCCCATATCCACCAt                                                                                                      >  1:174323/1‑141 (MQ=255)
                                                ggTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTcc                                                                                           >  1:59581/1‑141 (MQ=255)
                                                   aTGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAATCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCGATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGt                                                                                        <  2:237280/141‑1 (MQ=255)
                                                    tGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGtt                                                                                       <  2:248975/141‑1 (MQ=255)
                                                      gCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTg                                                                                     <  1:341230/141‑1 (MQ=255)
                                                       cATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGa                                                                                    <  2:220490/141‑1 (MQ=255)
                                                           ggTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCaa                                                                                >  1:98945/1‑141 (MQ=255)
                                                                   gctTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCa                                                                        <  2:325001/141‑1 (MQ=255)
                                                                                 agctgcATTACATGCAAGCCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAAt                                                          <  2:359832/140‑1 (MQ=255)
                                                                                          aCATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTAc                                                 <  2:198064/141‑1 (MQ=255)
                                                                                             tCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAAc                                              >  1:243854/1‑141 (MQ=255)
                                                                                             tCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAAc                                              >  1:89818/1‑141 (MQ=255)
                                                                                                 aaCCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGc                                          <  2:35702/141‑1 (MQ=255)
                                                                                                      ttCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGc                                     >  2:309626/1‑141 (MQ=255)
                                                                                                           aTAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATa                                <  2:349262/141‑1 (MQ=255)
                                                                                                           aTAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATa                                <  2:105496/141‑1 (MQ=255)
                                                                                                           aTAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATa                                <  1:4956/141‑1 (MQ=255)
                                                                                                             aGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACg                              >  1:393000/1‑141 (MQ=255)
                                                                                                                 gCACGATATCATTACCTTGTAATCCGATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATTGCAGTTGATAAACTAATTAATATTACAACCGGCAATGCCCATACGAGCg                          <  1:272576/141‑1 (MQ=255)
                                                                                                                 gCACGATATCATTACCTTGTAATCCGATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATTGCAGTTGATAAACTAATTAATATTACAACCGGCAATGCCCATACGAGCg                          <  1:272600/141‑1 (MQ=255)
                                                                                                                  cACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATAATTTTTGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACGAGCGa                         <  1:212991/141‑1 (MQ=255)
                                                                                                                            ttACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACGAGCGAATATGTTCCa               >  1:69151/1‑141 (MQ=255)
                                                                                                                               ccTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACGAGCGAATATGTTCCAAGt            <  2:445730/141‑1 (MQ=255)
                                                                                                                               ccTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACGAGCGAATATGTTCCAAGt            >  1:330384/1‑141 (MQ=255)
                                                                                                                                    tAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACGAGCGAATATGTTCCAAGTGATAg       >  2:424975/1‑141 (MQ=255)
                                                                                                                                    tAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACGAGCGAATATGTTCCAAGTGATAg       >  2:453247/1‑141 (MQ=255)
                                                                                                                                       tCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACGAGCGAATATGTTCCAAGTGATAGcac    <  1:304137/141‑1 (MQ=255)
                                                                                                                                         cTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACGAGCGAATATGTTCCAAGTGATAGCacac  <  2:443905/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
ACAAGAACTACCAAAGCTTATTAAACTCATACACTGTGTCTTCGTGCAGGTATGGCATTGGTGTGCTGCTTGTGAAATAGCTGCTGCATTACATCCAAACCCTTCTAATAGTGGCACGATATCATTACCTTGTAATCCTATCCATAACATTGACGGTTCAATTGCCCATATCATCCATGATTTGAGTCCTGTTTGATCAATTATAGCAGTTGATAAACTAATCAATATTACAACCGGCAATGCCCATACGAGCGAATATGTTCCAAGTGATAGCACAC  >  CP000730/2723487‑2723764

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