Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,550,546 C→T 7.5% A263V (GCC→GTC)  PP_4881 → iron ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,550,5460CT 7.5% 153.9 / 2.2 54A263V (GCC→GTC) PP_4881iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (27/22);  new base T (2/2);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

TGGCTGCAGTACTGGTTAACAACTACTACTGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAG  >  NC_002947/5550408‑5550688
                                                                                                                                          |                                                                                                                                              
tGGCTGCAGTACTGGTTAACAACTACTACTGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAacac                                                                                                                                      >  2:341435/1‑149 (MQ=255)
         tACTGGTTAACAACTACTACTGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCa                                                                                                                                              >  4:271547/1‑132 (MQ=255)
         tACTGGTTAACAACTACTACTGGTACACACTTGAACACGAACGCGACAAGCTGAACACCAAAAAAAACAACCTGGCCGAAGCCAAAGAAGGCAACCTGGAAAAAAAAACCGGCGCAGGCGTGAACAAAGACCGCAAACCCACAAaagaa                                                                                                                             >  8:221119/1‑149 (MQ=255)
           cTGGTTAACAACTACTACTGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGc                                                                                                                           <  2:286142/149‑1 (MQ=255)
             ggTTAACAACTACTACTGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGcc                                                                                                                          >  6:87977/1‑148 (MQ=255)
              tttAAAAACTAATAATTGTACGCACTTGAACGCGAATGGTTCAATTTGGTAACCAATCTCTAATACCTGGTCTTCGGCGTTTCTGGCAACCTTGTCACTATTTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGccc                                                                                                                         <  7:221119/147‑1 (MQ=255)
                     actactactGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTg                                                                                                                 <  1:234951/149‑1 (MQ=255)
                         ctactGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCt                                                                                                               <  6:153232/147‑1 (MQ=255)
                         ctactGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTc                                                                                                              >  4:292738/1‑148 (MQ=255)
                            ctGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTTGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACt                                                                                                          <  3:292738/149‑1 (MQ=255)
                             tGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGTCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTg                                                                                                         >  4:261145/1‑149 (MQ=255)
                             tGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTg                                                                                                         >  6:14465/1‑149 (MQ=255)
                             tGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCAGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTg                                                                                                         >  5:313022/1‑149 (MQ=255)
                                   gCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGTGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCaaa                                                                                                                                         <  1:91592/111‑1 (MQ=255)
                                   gCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCaaa                                                                                                                                         >  2:91592/1‑111 (MQ=255)
                                   gCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGc                                                                                                   >  1:329554/1‑149 (MQ=255)
                                     aCTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGc                                                                                                   >  6:304424/1‑147 (MQ=255)
                                     aCTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGc                                                                                                   >  1:18078/1‑147 (MQ=255)
                                     aCTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCa                                                                                                 <  5:304424/149‑1 (MQ=255)
                                       ttGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGTCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGc                                                                                               >  6:235822/1‑149 (MQ=255)
                                         gaacgcgaacgcGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCg                                                                                              <  3:157612/148‑1 (MQ=255)
                                           acgcgaacgcgGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCg                                                                                              >  4:218123/1‑146 (MQ=255)
                                              cgaacgcgGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGTCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGaagaag                                                                                        <  5:235822/149‑1 (MQ=255)
                                               gaacgcgGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGc                                                                                                                <  2:294740/124‑1 (MQ=255)
                                               gaacgcgGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGc                                                                                                                >  1:294740/1‑124 (MQ=255)
                                                aacgcgGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGc                                                                                      <  5:87977/149‑1 (MQ=255)
                                                    cgGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAgc                                                                                  <  1:341435/149‑1 (MQ=255)
                                                        aaGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGt                                                                              <  3:218123/149‑1 (MQ=255)
                                                          gCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAgc                                                                                  <  4:217176/143‑1 (MQ=255)
                                                          gCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAgc                                                                                  >  3:217176/1‑143 (MQ=255)
                                                            tGGACACCAAGCTATACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGAt                                                                           >  4:256096/1‑148 (MQ=255)
                                                                cacCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCAccc                                                                      <  2:329554/149‑1 (MQ=255)
                                                                   cAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGt                                                                              >  8:304447/1‑138 (MQ=255)
                                                                   cAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGt                                                                              <  7:304447/138‑1 (MQ=255)
                                                                   cAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGa                                                                   >  3:289725/1‑149 (MQ=255)
                                                                          tactacCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGc                                                                                      >  7:318081/1‑123 (MQ=255)
                                                                          tactacCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGc                                                                                      <  8:318081/123‑1 (MQ=255)
                                                                          tactacCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCAccgc                                                            >  2:17420/1‑149 (MQ=255)
                                                                               ccTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGa                                                         <  4:289725/147‑1 (MQ=255)
                                                                                 tGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTAc                                                     >  3:68909/1‑149 (MQ=255)
                                                                                    ccGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCa                                                  <  3:256096/149‑1 (MQ=255)
                                                                                            gATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGTCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACaa                                          <  3:261145/149‑1 (MQ=255)
                                                                                            gATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGcac                                            >  1:22448/1‑147 (MQ=255)
                                                                                            gATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACaa                                          <  5:14465/149‑1 (MQ=255)
                                                                                                  ggCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCAt                                    >  7:112708/1‑149 (MQ=255)
                                                                                                   gCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATg                                   <  6:313022/149‑1 (MQ=255)
                                                                                                       ccTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCCCTGCACAAGGGCATGGttt                               >  2:23309/1‑149 (MQ=255)
                                                                                                        cTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTc                              <  2:18078/149‑1 (MQ=255)
                                                                                                          ggTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCt                             >  4:310538/1‑148 (MQ=255)
                                                                                                              aCTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACcg                        >  1:135960/1‑149 (MQ=255)
                                                                                                                tATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACcgcg                      >  1:19339/1‑149 (MQ=255)
                                                                                                                              gCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTc        >  8:324939/1‑149 (MQ=255)
                                                                                                                                gTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTcga      <  3:310538/149‑1 (MQ=255)
                                                                                                                                    tCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAg  >  5:202457/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                              
TGGCTGCAGTACTGGTTAACAACTACTACTGGTACGCACTTGAACGCGAACGCGGCAAGCTGGACACCAAGCTCTACTACCTGGCCGACGGCGATGCTGGCAACCTGGTCACTATCTCCGGCGCAGGCGTGGTCAAGGCCAGCAAACACCCGAAAGAAGCCCAGGCCCTGCTCAACTGGATGGCCAGCGAAGAAGGCCAGCGCGTGATCACCCAGACCACCGCCGAGTACCCACTGCACAAGGGCATGGTTTCTGACCGCGGCCTGAAACCGTTCGACGAG  >  NC_002947/5550408‑5550688

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