Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,762,710 T→A 10.4% K62I (AAA→ATA)  PP_2416 ← iron ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,762,7100TA10.4% 123.8 / 3.9 46K62I (AAA→ATA) PP_2416iron ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base A (3/2);  total (25/21)
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

CCAAGGCAACCGCGTCGAACACGCTGATCGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGCAGCCCCAGCCCATGGCAG  >  NC_002947/2762570‑2762856
                                                                                                                                            |                                                                                                                                                  
ccAAGGCAACCGCGTCGAACACGCTGATCGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCAACCAGGACCTGCACGCTGCCGCAGGCGGGTATGCGCAGc                                                                                                                                            >  6:177753/1‑149 (MQ=255)
  aaGGCAACCGCGTCGAACACGCTGATCGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGccc                                                                                                                                          >  5:331608/1‑149 (MQ=255)
     gCAACCGCGTCGAACACGCTGATCGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGc                                                                                                                                       <  5:323841/149‑1 (MQ=255)
           gcgTCGAACACGCTGATCGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCaa                                                                                                                                 <  4:276369/149‑1 (MQ=255)
               cGAACACGCTGATCGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCCGCCAGCAACtt                                                                                                                              >  4:45455/1‑148 (MQ=255)
                 aaCACGCTGATCGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGttt                                                                                                                                                   <  2:270055/125‑1 (MQ=255)
                  acacGCTGATCGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAAt                                                                                                                          <  3:221784/149‑1 (MQ=255)
                          aTCGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCCGCCAGCAACTTCAATAGCGAg                                                                                                                    <  3:45455/147‑1 (MQ=255)
                            cGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGAtt                                                                                                                >  7:168528/1‑149 (MQ=255)
                              cATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTa                                                                                                              >  8:176756/1‑149 (MQ=255)
                                tCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTAcc                                                                                                            >  3:114860/1‑149 (MQ=255)
                                         tCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGtt                                                                                                   >  8:193286/1‑149 (MQ=255)
                                          cGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTg                                                                                                  >  7:219290/1‑149 (MQ=255)
                                                      cGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATg                                                                                      >  6:167135/1‑149 (MQ=255)
                                                       gACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGc                                                                                     <  8:219290/149‑1 (MQ=255)
                                                       gACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCACGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGCGGCcgaccaagc                                                                                     >  6:29005/1‑144 (MQ=255)
                                                               ccagcgccTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGtt                                                                             <  7:176756/149‑1 (MQ=255)
                                                                         gggaaCACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGc                                                                     <  5:29005/144‑1 (MQ=255)
                                                                            aaCACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTATGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTc                                                                           <  7:222464/138‑1 (MQ=255)
                                                                            aaCACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTATGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTc                                                                           <  7:222463/138‑1 (MQ=255)
                                                                            aaCACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTATGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTc                                                                           >  8:222463/1‑138 (MQ=255)
                                                                            aaCACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTATGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTc                                                                           >  8:222464/1‑138 (MQ=255)
                                                                                cggcggcggGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAAcgc                                                             >  8:289163/1‑148 (MQ=255)
                                                                                 ggcggcggGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAAcgcg                                                            >  3:88930/1‑148 (MQ=255)
                                                                                      cggGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTc                                                      <  8:168528/149‑1 (MQ=255)
                                                                                          gCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGata                                                  >  5:321872/1‑149 (MQ=255)
                                                                                          gCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGata                                                  <  5:167135/149‑1 (MQ=255)
                                                                                                gtgCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCtt                                            >  1:294399/1‑148 (MQ=255)
                                                                                                 tgCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCg                                               <  4:114860/145‑1 (MQ=255)
                                                                                                 tgCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCt                                             <  7:286104/147‑1 (MQ=255)
                                                                                                              cccAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTcc                              <  7:193286/149‑1 (MQ=255)
                                                                                                               ccAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCg                             >  6:34722/1‑149 (MQ=255)
                                                                                                                  ggAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGc                           >  3:15349/1‑148 (MQ=255)
                                                                                                                      cTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGtt                      >  5:256563/1‑149 (MQ=255)
                                                                                                                      cTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGtt                      >  2:44414/1‑149 (MQ=255)
                                                                                                                             cTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGcagccc               <  2:294399/149‑1 (MQ=255)
                                                                                                                             cTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCCGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGcagccc               <  5:177753/149‑1 (MQ=255)
                                                                                                                              tGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGcagccc               >  7:305688/1‑148 (MQ=255)
                                                                                                                                ccGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTc                                                                                                                             >  1:252194/1‑36 (MQ=255)
                                                                                                                                ccGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTc                                                                                                                             <  2:252194/36‑1 (MQ=255)
                                                                                                                                    ttgCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGCAGCCCCAGCCCa        <  3:248138/147‑1 (MQ=255)
                                                                                                                                      gCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGCAGCCCCAGCCCa        >  4:248138/1‑147 (MQ=255)
                                                                                                                                      gCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGCAGCCCCAGCCCATg      >  6:267686/1‑149 (MQ=255)
                                                                                                                                       cGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGCAGCCCCAGCCCATgg     <  7:289163/149‑1 (MQ=255)
                                                                                                                                          ttttGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGTGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGCAGCCCCAGCCCATGGCAg  <  6:331608/148‑1 (MQ=255)
                                                                                                                                          gTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCa                         >  8:173892/1‑126 (MQ=255)
                                                                                                                                          gTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCa                         <  7:173892/126‑1 (MQ=255)
                                                                                                                                          gTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGCAGCCCCAGCCCATGGCAg  <  1:44414/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CCAAGGCAACCGCGTCGAACACGCTGATCGCATCAAGTGTGTCGGCCTGTTGCTCGACCAGCGCCAGCGCCTGGGCAACACGGCGGCGGGGCATCTGTGCCAGTGGCTCACCCAGGAGCTGCACGCTGCCGCAGGCGGGTTTGCGCAGCCCTGCCAGCAACTTCAATAGCGAGGATTTACCCGAGCCGTTGGGGCCGACGATGCCCAGGGTTTCACCGGCGACAACGCGAAGGTCGATATCGCTGAGCACGGTATTTCCGGCCAGTTGCAGCCCCAGCCCATGGCAG  >  NC_002947/2762570‑2762856

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