Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,227,331:1 +T 10.3% coding (1776/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,3311.T10.3% 209.9 / 19.7 58coding (1776/2391 nt)PP_4606ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base . (29/23);  new base T (4/2);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

GTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTC  >  NC_002947/5227191‑5227477
                                                                                                                                             |                                                                                                                                                   
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                                                                                  >  2:1336547/1‑85 (MQ=255)
                                                                                                                          gCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTg                                                                                  <  1:1336547/85‑1 (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                                              <  1:570748/119‑1 (MQ=255)
                                                                                                                            cctcctCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTg                                              >  2:570748/1‑119 (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)
                                                                                                                                         cACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTc  >  2:118527/1‑150 (MQ=255)
                                                                                                                                          aCCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGc                                                             >  2:1149701/1‑90 (MQ=255)
                                                                                                                                          aCCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGc                                                             <  1:1149701/90‑1 (MQ=255)
                                                                                                                                           ccGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGAc                                                                                        >  1:920660/1‑62 (MQ=255)
                                                                                                                                           ccGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGAc                                                                                        <  2:920660/62‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GTCTCGCTTGTGACGGCGTCAAAGTTCGAGACCAGTACATTGCGCTTGGTGATATCGAACAATGCCAGGGTGCCTGTGATGCGGCCTGGCATGTCGAGCTTGGCGCCCAGCTCCCACGCCTTGCCCTCCTCAGGTGCCACCGAGGAGTCGAGCACGACATTGCCGGTCAGCGGCGCGATGGTCGAGTTGGGCTTGAACGACTCGCTGTAGCTGCCATAGAACGAGAGCTGGTCGTCGACTTTGTAAACGATACCGGCGTGAGGTACCCAGGCCTGGCCACTGATGTC  >  NC_002947/5227191‑5227477

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