Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,573,250 G→A 9.9% R699C (CGC→TGC)  PP_3155 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,573,2500GA 9.9% 132.9 / 3.6 50R699C (CGC→TGC) PP_3155ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (24/21);  new base A (2/3);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATC  >  NC_002947/3573103‑3573392
                                                                                                                                                   |                                                                                                                                                  
aTGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTgcgc                                                                                                                                                   >  3:173570/1‑149 (MQ=255)
        ctCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑Gcc                                                                                                                                         >  7:131342/1‑149 (MQ=255)
        ctCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcagg                                                                                                                                           >  7:346103/1‑146 (MQ=255)
         tCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCg                                                                                                                                        >  7:177990/1‑149 (MQ=255)
           cAACGCCGGTATAAAGCAGCCAAACGAGGGCATAGAGTGTTACTGGTTGTTATACAGAACAGCGTTCTTACATTACGGGGCCACTGCGCGGCCCGGGCATTTCTTATGTTATAGTAAACGCCGGCCTGCACCGCGCACGGTgctggccg                                                                                                                                        <  2:314721/149‑7 (MQ=255)
           cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCggg                                                                                                                                      >  2:278526/1‑149 (MQ=255)
           cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcaggccg                                                                                                                                        <  8:346103/149‑7 (MQ=255)
           cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGCGCACGGTgcaggc                                                                                                                                          <  3:144529/147‑5 (MQ=255)
             aCGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGt                                                                                                                                    <  8:177990/149‑1 (MQ=255)
                         agcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGt                                                                                                                        <  8:131342/149‑1 (MQ=255)
                             gcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAACGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGaa                                                                                                                    >  3:283852/1‑149 (MQ=255)
                              ccAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAAc                                                                                                                   <  1:36255/149‑1 (MQ=255)
                                aaaCGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAAcgc                                                                                                                 <  1:278526/149‑1 (MQ=255)
                                 aaCGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAAcgcg                                                                                                                >  2:105858/1‑149 (MQ=255)
                                          aTAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCCGGGTATACATAGCTGTTGAACGCGCCCTCGTCa                                                                                                       >  3:54039/1‑149 (MQ=255)
                                            agagTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑a                                                                                                   >  3:241459/1‑149 (MQ=255)
                                               gTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑Acc                                                                                                 >  7:338069/1‑148 (MQ=255)
                                                     cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCttgtt                                                                                          >  3:204751/1‑149 (MQ=255)
                                                            ttCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAaca                                                                                     <  4:204751/147‑1 (MQ=255)
                                                                   gcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCa                                                                            <  4:120381/149‑1 (MQ=255)
                                                                           ttCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGa                                                                    <  6:251808/149‑1 (MQ=255)
                                                                              tGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGAtt                                                                  >  1:20468/1‑148 (MQ=255)
                                                                                aCATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGtat                                                                                                                                  >  3:267389/1‑84 (MQ=255)
                                                                                aCATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGtat                                                                                                                                  <  4:267389/84‑1 (MQ=255)
                                                                                aCATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCaa                                                               >  7:91850/1‑149 (MQ=255)
                                                                                aCATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCaa                                                               >  2:108208/1‑149 (MQ=255)
                                                                                          gacACTGATCGGCACCGGCAATTCAGAAGGTATAGGAAAAGCCGGCCTGCACCGTGCGCGGTGC‑‑GCAGGGGTATACATAGATGTTGAACGCGCCCTCGTCAT‑‑ACCACTTGTTGAACACGTTTTTCACTTCCAGATTAAAAAGCACgtgt                                                     >  2:56529/3‑149 (MQ=255)
                                                                                           ccACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGtt                                                    >  8:162863/1‑149 (MQ=255)
                                                                                            cACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTc                                                   >  1:249912/1‑149 (MQ=255)
                                                                                             aCTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCCGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCg                                                  <  4:54039/149‑1 (MQ=255)
                                                                                                 cgcgGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGtt                                                <  2:260228/147‑1 (MQ=255)
                                                                                                  gcgGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGAc                                             <  2:20468/149‑1 (MQ=255)
                                                                                                    ggCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCt                                           >  5:330835/1‑149 (MQ=255)
                                                                                                      ccccGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTg                                         <  4:241459/149‑1 (MQ=255)
                                                                                                        ccGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGCCCTTg                                         >  2:118215/1‑147 (MQ=255)
                                                                                                         cGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATTCACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGt                                        <  1:105858/149‑1 (MQ=255)
                                                                                                         cGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTg                                         >  8:348139/1‑146 (MQ=255)
                                                                                                              aTTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGttttt                                                                              <  3:285423/104‑1 (MQ=255)
                                                                                                              aTTCAGAAGGTATAGGAAACGCAGGCCAGCACCGTGCGCCGTGC‑‑GCCGGGGTATACATAGCTGTTTAACGCGCCCTCGACAT‑‑ACCCCTTGTTGAAAAAGttttt                                                                              >  4:285423/1‑104 (MQ=255)
                                                                                                               ttCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTc                                <  1:108208/149‑1 (MQ=255)
                                                                                                                   gAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCACTCGTCAT‑‑ACCCCTTGTTGCACACGTTTTTCACGTCCAGCTTCAACCGAACGTATTCGTTGACATTGTAGAAGCTaaa                              >  4:63532/1‑147 (MQ=255)
                                                                                                                    aaGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAAc                                                              <  2:84500/114‑1 (MQ=255)
                                                                                                                    aaGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAAc                                                              >  1:84500/1‑114 (MQ=255)
                                                                                                                      ggTATAGGAAACGCCTGCCTGAACCTTTCGCGGTGC‑‑GACGTGGTATACATAGCTGTTTAACGCGCCCTCGTCAT‑‑ACCCATTGTTGAACACGTTTTTCACGTCAAGATTCAAACGCACGTTTTCGTTGACCTTGTAGAAGCTCAAAAGGt                         <  1:56529/149‑1 (MQ=255)
                                                                                                                      ggTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGt                         <  7:348139/149‑1 (MQ=255)
                                                                                                                          taGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGAc                     <  6:330835/149‑1 (MQ=255)
                                                                                                                            ggAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACACCTTGTTGAACACGTTTTTCACGTCCAGATTGAACAGCACGTGTTCGTTGAACTTGTAGAAGCTCAACAGGTCGACCa                   >  3:275724/1‑149 (MQ=255)
                                                                                                                             gAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGATCACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCAc                  >  7:314914/1‑149 (MQ=255)
                                                                                                                                   ccGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTa            <  8:91850/149‑1 (MQ=255)
                                                                                                                                     ccggcgCACCGCGCACGGTGCAGGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTg              >  5:97322/6‑147 (MQ=255)
                                                                                                                                             ccGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATc  <  8:338069/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGC‑‑GCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCAT‑‑ACCCCTTGTTGAACACGTTTTTCACGTCCAGATTCAACCGCACGTGTTCGTTGACCTTGTAGAAGCTCAACAGGTCGACCACGCTGTAGCGCTCCATC  >  NC_002947/3573103‑3573392

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: ‑