Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,760,878 G→C 9.2% Q628H (CAG→CAC)  fpvA_3 → Ferripyoverdine receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,760,8780GC 9.2% 148.3 / 6.5 52Q628H (CAG→CAC) fpvA_3Ferripyoverdine receptor
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base C (3/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

ATATCTTCATGCCGCAATCGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCACCCACGGTATCGACCTGGAGTTGATGGGCGCACTG  >  KT2440/4760745‑4761001
                                                                                                                                     |                                                                                                                           
atatCTTCATGCCGCAATGGGAACTGGCGCGGGACCGCTCCATCATCCGCCCAATCGACGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGa                                                                                                                            <  7:148378/135‑1 (MQ=255)
   tCTTCATGCCGCAATCGGAACTGGCGCGTGACCGCGCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGAGCACTCAACGCCAGCGCGGCGATCTTCCAGATCg                                                                                                                         >  1:318835/1‑135 (MQ=255)
    cTTCATGCCGCAATCGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCACATCGa                                                                                                                        >  7:37896/1‑135 (MQ=255)
    cTTCATGCCGCAATCGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCACATCGa                                                                                                                        >  7:31858/1‑135 (MQ=255)
         tGCCGCAATCGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAg                                                                                                                    >  8:278653/1‑134 (MQ=255)
         tGCCGCAATCGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAg                                                                                                                    >  7:29838/1‑134 (MQ=255)
           ccGCAATCGGAACTGGCGCGTGCCCGCTCCATAATCCGCCCAATCGAAGTCAAGAACTACGAGATAGGCCTCAAGGGCGAGTACTTCGACGGCGCACACAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGa                                                                                                                  >  5:18990/1‑134 (MQ=255)
           ccGCAATCGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCAACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGaa                                                                                                                 >  4:366121/1‑135 (MQ=255)
              cAATCGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGCTCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAAcc                                                                                                             <  3:151086/136‑1 (MQ=255)
                 tCGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACcgcg                                                                                                          >  4:171499/1‑136 (MQ=255)
                 tCGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACcgcg                                                                                                          >  4:106326/1‑136 (MQ=255)
                  cGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACcgcg                                                                                                          >  2:70752/1‑135 (MQ=255)
                      aCTGGCGCGTGACCGCTCCATCAGCCCCCCAATCCAAGCCAAGAACAACGAGACCGGCACCAAGGGCGAGTACTTCCACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGcgccgc                                                                                                      <  5:115894/135‑1 (MQ=255)
                        tGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCCGCGCGGCGCTCTTCCAGATCGACCAGGAATCCCGCGCCGCAg                                                                                                    >  3:221081/1‑135 (MQ=255)
                                  ccGCTCCATCATCCCCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAAc                                                                                         <  3:219594/136‑1 (MQ=255)
                                       ccATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAg                                                                                      >  7:256861/1‑134 (MQ=255)
                                        catcatCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCACATCGCCCAGGAAAACCGCGCCGCAGAAGCGTCTAACCAg                                                                                      >  1:149154/1‑133 (MQ=255)
                                                     aTCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACa                                                                      >  4:103541/1‑136 (MQ=255)
                                                         aaGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATc                                                                    <  8:29838/134‑1 (MQ=255)
                                                         aaGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCAc                                                                  <  4:312117/136‑1 (MQ=255)
                                                               aGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGAACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCt                                                             <  1:263477/135‑1 (MQ=255)
                                                                  aCTACGATAGCGGCCGCACGGGCGAGTCCTGCGACGGCGCAGTCAACTCCAGCGCGGCTATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGa                                                         <  3:37627/136‑1 (MQ=255)
                                                                  aCTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGa                                                         <  3:174178/136‑1 (MQ=255)
                                                                   cTACGAGAGCGGCAGCATGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGa                                                         <  1:134021/135‑1 (MQ=255)
                                                                      cGAGATCGGCATCAAGGGCGAGTACTTCGAAGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCATGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCTACATCACCTGCTACGAAg                                                       <  4:60680/134‑1 (MQ=255)
                                                                       gagaTCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAg                                                       <  7:223210/133‑1 (MQ=255)
                                                                       gagaTCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAg                                                       <  7:278653/133‑1 (MQ=255)
                                                                         gaTCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCg                                                  <  3:106326/136‑1 (MQ=255)
                                                                         gaTCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCg                                                  <  3:171499/136‑1 (MQ=255)
                                                                          agcGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGAGCGACAAGGAAACCCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCgg                                                 <  3:44115/134‑1 (MQ=255)
                                                                            cGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGAGCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCggg                                                <  6:217660/135‑1 (MQ=255)
                                                                               cATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGaa                                                        <  8:235363/124‑1 (MQ=255)
                                                                               cATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGaa                                                        >  7:235363/1‑124 (MQ=255)
                                                                                aTCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAg                                            <  5:77385/135‑1 (MQ=255)
                                                                                aTCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAg                                            <  8:339325/135‑1 (MQ=255)
                                                                                 tCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAg                                            >  8:119034/1‑134 (MQ=255)
                                                                                 tCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGCTTGCGTCGCCATCACCTGCAACGAAGCCTCGGGCAAg                                            >  8:307553/1‑134 (MQ=255)
                                                                                 tCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTACTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAg                                            >  8:93241/1‑134 (MQ=255)
                                                                                       gCGATTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCAc                                    <  3:215446/136‑1 (MQ=255)
                                                                                       gCGAGTACTTCAACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCAc                                    <  3:366121/136‑1 (MQ=255)
                                                                                        cGAGTACTTCGACGGCGCACTCACCGCCAGCGCGGCGAGCTTCCAGATCGACCAGGAAACCCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCAc                                    <  6:184451/135‑1 (MQ=255)
                                                                                               ttCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCACATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCACCCACgg                              <  8:31858/134‑1 (MQ=255)
                                                                                               ttCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCACATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCACCCACgg                              <  8:37896/134‑1 (MQ=255)
                                                                                                   tcGGCTCACTCACCGCCAGAGCGGCTTTCTTCAAGATCGAACAGGAAAACCGCGCCTCAGAATCTTCTACAAAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCTCCCACGGTATCg                         <  6:18990/134‑1 (MQ=255)
                                                                                                   aCGGCGCACTCAACGCCAGCGCGGCGATCTTCCCGATCGAGAAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAATGTTGCGTCGACTTCTCCTGCTACGAAGCCTCGGGCAAGGTACGCACCCACGGTATCg                         <  8:121005/135‑1 (MQ=255)
                                                                                                    cGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCACCCACGGTATCGAc                       <  3:363612/136‑1 (MQ=255)
                                                                                                          aCTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCACCCACGGTATCGACCTGGa                  >  7:360587/1‑135 (MQ=255)
                                                                                                          aCTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTATAACCAGGAAGCTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAACGTACGCACCCACGGTATCGACCTGGa                  >  8:304069/1‑135 (MQ=255)
                                                                                                                 gCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCACCCACGGTATCGACCTGGAGTTGATg           >  8:357472/1‑135 (MQ=255)
                                                                                                                    aGCGCGGCGATCTTCCATATCGACCATGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGAACCCACGGTATCGACCGGGAGTTGATGGgc        >  1:116392/1‑135 (MQ=255)
                                                                                                                         ggcgATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGCAGCTTCTAACCAGGAAGGTTGCGTCGCCATCACCTGCTACGAAGCCTCGGGCAAGGTCCGCACCCACGGTATCGACCTGGAGTTGATGGGCGCACt   >  1:249486/1‑135 (MQ=255)
                                                                                                                         ggcgATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCACCCACGGTATCGACCTGGAGTTGATGGGCGCACTg  >  3:179498/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                           
ATATCTTCATGCCGCAATCGGAACTGGCGCGTGACCGCTCCATCATCCGCCCAATCGAAGGCAAGAACTACGAGATCGGCATCAAGGGCGAGTACTTCGACGGCGCACTCAACGCCAGCGCGGCGATCTTCCAGATCGACCAGGAAAACCGCGCCGCAGAAGCTTCTAACCAGGAAGGTTGCGTCGACATCACCTGCTACGAAGCCTCGGGCAAGGTACGCACCCACGGTATCGACCTGGAGTTGATGGGCGCACTG  >  KT2440/4760745‑4761001

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