Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,227,322 Δ1 bp 10.3% coding (1785/2391 nt) PP_4606 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,227,3220A.10.3% 165.0 / 15.2 58coding (1785/2391 nt)PP_4606ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (27/25);  new base . (4/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

CGCATTGCTGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCAC  >  NC_002947/5227175‑5227470
                                                                                                                                                   |                                                                                                                                                      
cGCATTGCTGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCATGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAgg                                                                                                                                                      <  1:786792/150‑1 (MQ=255)
 gCATTGCTGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGt                                                                                                                                                     <  1:248658/150‑1 (MQ=255)
 gCATTGCTGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGt                                                                                                                                                     <  1:1206662/150‑1 (MQ=255)
                gTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCg                                                                                                                                     >  2:271833/1‑150 (MQ=255)
                         gtgACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACa                                                                                                                            <  2:828168/150‑1 (MQ=255)
                         gtgACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACa                                                                                                                            <  1:1596319/150‑1 (MQ=255)
                            aCGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTg                                                                                                                         <  2:336718/150‑1 (MQ=255)
                                    aaaGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAg                                                                                                                 >  1:1002357/1‑150 (MQ=255)
                                       gTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAgcgg                                                                                                              >  2:1609989/1‑150 (MQ=255)
                                        ttCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAgcggc                                                                                                             >  1:118638/1‑150 (MQ=255)
                                            agaCCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGACAGCGAGGAGTCATGCACGTCAGTGCCGGTCAGCGGcgcgc                                                                                                         >  2:14595/1‑149 (MQ=255)
                                                            cTTGGTGATATCGAACCATGCTAGCGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCtt                                                                                        <  2:427251/150‑1 (MQ=255)
                                                             ttGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTg                                                                                       >  2:565743/1‑150 (MQ=255)
                                                               ggTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGaa                                                                                     <  2:1153124/150‑1 (MQ=255)
                                                                gTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAAc                                                                                    <  1:271833/150‑1 (MQ=255)
                                                                  gaaatCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCGGTGGCACCTgagga                                                                                                                                         >  1:1645854/4‑96 (MQ=255)
                                                                  gATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCGGTGGCACCTgagga                                                                                                                                         <  2:1645854/96‑1 (MQ=255)
                                                                        gAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCt                                                                            >  1:1613589/1‑150 (MQ=255)
                                                                         aaCAATTCCAGGTTGCCTGTGATGCGGCCTGGCATGTCGAGCGTGGCGCCCAGCTCCCACGCCTTACCCTCCGCAGGTGCCACCGAAGAGGCGAGCAAGACATGGCAGGACAGCGGCGCCGTGGGCGAGTTGGGCTTTAACGACTCGCTg                                                                           >  2:153454/1‑150 (MQ=255)
                                                                         aaCAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTg                                                                           >  1:104521/1‑150 (MQ=255)
                                                                                 cAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGcc                                                                   >  2:1155168/1‑150 (MQ=255)
                                                                                 cAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGAACTCGAGCACGAGATTTCCGGTCAGCGGCGCGATGGTAGAGTTGGGTTTGAACGACTCGCTGTAGCTGcc                                                                   >  1:1263277/1‑150 (MQ=255)
                                                                                      tGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGa                                                              <  2:193728/150‑1 (MQ=255)
                                                                                         cTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCCACCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACg                                                           <  1:14595/150‑1 (MQ=255)
                                                                                          tgtgATGCGGCCTGGCATGTCGAGCTTGGCGCTCAGCTCCCACGCCTTGCCCTCCTCGGTGGCACCTgaggag                                                                                                                                        <  2:795884/73‑1 (MQ=37)
                                                                                          tgtgATGCGGCCTGGCATGTCGAGCTTGGCGCTCAGCTCCCACGCCTTGCCCTCCTCGGTGGCACCTgaggag                                                                                                                                        >  1:795884/1‑73 (MQ=37)
                                                                                          tgtgATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGgcgc                                                                                                           <  2:893485/102‑1 (MQ=255)
                                                                                          tgtgATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGgcgc                                                                                                           >  1:893485/1‑102 (MQ=255)
                                                                                           gtgATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACgag                                                         >  1:1314254/1‑150 (MQ=255)
                                                                                            tgATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACgaga                                                        >  2:850319/1‑150 (MQ=255)
                                                                                              aTGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGc                                                      <  2:104521/150‑1 (MQ=255)
                                                                                               tGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCt                                                     <  1:565743/150‑1 (MQ=255)
                                                                                                gCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTg                                                    <  2:392859/150‑1 (MQ=255)
                                                                                                         cATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCAACCGAGGAGTCGAGCACGACATTGCCGGTCCACGGCGTGTTGGTACAGGTGCGGGTTGTCGACGCTCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGAc                                           <  1:850319/150‑1 (MQ=255)
                                                                                                          aTGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACt                                          >  1:695691/1‑150 (MQ=255)
                                                                                                           tGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAAAAAATCGTTGTATCGGCCATAGAACGAGAGCTGGTCGTCGACtt                                         >  2:793438/1‑150 (MQ=255)
                                                                                                            gTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACttt                                        <  1:1155168/150‑1 (MQ=255)
                                                                                                              cGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGt                                      <  2:1002357/150‑1 (MQ=255)
                                                                                                                aGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACGTTGTaa                                    >  1:235945/1‑150 (MQ=255)
                                                                                                                                 ccACGCCTTGCCCTCCTCAGGTGCAACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAgg                   >  1:1212020/1‑150 (MQ=255)
                                                                                                                                  cACGCCTTGCCCTCCTCAGGTGCCACCGGGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGt                  >  2:1596180/1‑150 (MQ=255)
                                                                                                                                   aCGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTa                 <  2:1314254/150‑1 (MQ=255)
                                                                                                                                        ttGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAg            >  2:1532087/1‑150 (MQ=255)
                                                                                                                                        ttGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGGTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGGGAGGGACCCAg            >  2:1474268/1‑150 (MQ=255)
                                                                                                                                         tcgactcctcGGTGGCACCTGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAgg           >  1:488139/5‑150 (MQ=255)
                                                                                                                                          gCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTg                                                                           <  1:1336547/85‑1 (MQ=255)
                                                                                                                                          gCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTg                                                                           >  2:1336547/1‑85 (MQ=255)
                                                                                                                                          gCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGc          >  2:1203607/1‑150 (MQ=255)
                                                                                                                                           gactcctcGGTGGCACCTGAGGAGTCGAGCACGAAATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGcc         >  1:251216/3‑150 (MQ=255)
                                                                                                                                           cccTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAgg           <  2:1291715/148‑1 (MQ=255)
                                                                                                                                           cccTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAgg           >  1:1291715/1‑148 (MQ=255)
                                                                                                                                            cctcctCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTg                                       >  2:570748/1‑119 (MQ=255)
                                                                                                                                            cctcctCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTg                                       <  1:570748/119‑1 (MQ=255)
                                                                                                                                             ctcctcAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGa                                                                                                         <  2:887923/53‑1 (MQ=255)
                                                                                                                                             ctcctcAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGa                                                                                                         >  1:887923/1‑53 (MQ=255)
                                                                                                                                                 tcAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGAt                               >  1:653704/1‑122 (MQ=255)
                                                                                                                                                 tcAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGAt                               <  2:653704/122‑1 (MQ=255)
                                                                                                                                                  cAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTAAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCAc  <  1:202045/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CGCATTGCTGTAGACCGTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCAC  >  NC_002947/5227175‑5227470

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