Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 883,744 T→G 22.3% A79A (GCT→GCG)  U712_04295 → putative siderophore transport system permease protein yfiZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881883,7440TG22.3% 57.1 / 6.4 56A79A (GCT→GCG) U712_04295putative siderophore transport system permease protein yfiZ
Reads supporting (aligned to +/- strand):  ref base T (14/28);  new base G (11/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACGTCATGGAGACAGGTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAGGCATTAGTATGGAGCT  >  CP006881/883608‑883877
                                                                                                                                        |                                                                                                                                     
gTACGTCATGGAGACAGGTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGctct                                                                                                                                     >  5:84029/1‑139 (MQ=255)
   cGTCATGGAGACAGGTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCAt                                                                                                                                  >  8:325310/1‑139 (MQ=255)
        tGGAGACAGGTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAgg                                                                                                                             <  7:325310/139‑1 (MQ=255)
        tGGAGACAGGTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCGCTCATGCAgg                                                                                                                             >  8:327316/1‑139 (MQ=255)
          gagaCAGGTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGgcg                                                                                                                           >  1:54164/1‑139 (MQ=255)
            gaCAGGTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTAATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCt                                                                                                                         <  3:311984/139‑1 (MQ=255)
             aCAGGTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCGTGCCGCCGCGGGCGCGCTCCTGCCCgcgcgc                                                                                                                        >  1:138969/1‑137 (MQ=255)
               aGGTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCATGCAGCCGCGGGCGCTCACCTGCAGGCGCtga                                                                                                                       >  8:16109/1‑136 (MQ=255)
               aGGTCTATCAGGCGTTTACTTCCGGCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCCGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCCGCGCCTCGCTTGCAGCCGCGGGCGCAGGCCTGCAGGCGCTcac                                                                                                                      >  2:26026/1‑139 (MQ=255)
                ggTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTcaca                                                                                                                     <  1:26026/139‑1 (MQ=255)
                      tCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCCGCGCCTCGCATGCCCCCCCGGGGGCGCGCCTGCAGGCCGGCACcaaaaa                                                                                                                >  6:135912/1‑138 (MQ=255)
                           cGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCGTGCCGCCGCGGGCGGGCGCCAGCCGGCGCGGACCAAAAATCCGca                                                                                                          >  4:355605/1‑138 (MQ=255)
                              ttACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCGCTCATGCCGGCGCGCACAAAAAATCAGCGTGc                                                                                                       >  2:109724/1‑139 (MQ=255)
                                aCTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTTCAGCCGCGGGCGCTCTCATGCAGGGGCGCACAAAAAATACCCTTGCGt                                                                                                     >  4:164389/1‑139 (MQ=255)
                                    ccttcAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCCGCCGCGGGCGCGCGCAGGGCGGCGCGCCCACAAAATCCGCGTGCGGCAcc                                                                                                 >  4:148588/1‑139 (MQ=255)
                                    ccttcAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCAcc                                                                                                 >  6:141329/1‑139 (MQ=255)
                                              aaCAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCCGCCGCGGGCGCGGGCCAGCCGGCGCTCACAAAAAAAACGCTTGCGTCCCCGGGGCttttg                                                                                       >  2:65314/1‑138 (MQ=255)
                                              aaCAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCCTGCCGGCGCTCCCAAAAAATCCGCTTGCGGCACCGGGGATTTTc                                                                                       >  7:107927/1‑139 (MQ=255)
                                              aaCAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTc                                                                                       <  3:325901/139‑1 (MQ=255)
                                                cAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTTCAGCCGCGGGCGCTCTACTTCCTGCGCTCCCAAAAAAACCGCGGGCGGCCCCGGGGAGTTTAgg                                                                                     >  4:148044/1‑139 (MQ=255)
                                                cAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCCTGCAGGCGCGCACAAAAAATACGGGTGCGTCACCGGGGATTTTCgg                                                                                     >  6:70922/1‑139 (MQ=255)
                                                cAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCgg                                                                                     >  5:83043/1‑139 (MQ=255)
                                                    tGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGAtt                                                                                 <  2:54164/139‑1 (MQ=255)
                                                                  tcaAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAgg                                                                   >  5:199952/1‑139 (MQ=255)
                                                                  tcaAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAgg                                                                   <  5:135912/139‑1 (MQ=255)
                                                                     aGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCACGCTTGCAGCCCCGGGCGCGCTCAAGGACGCGCGACCAAAAAATACGCGTTCCGCACCCGGGGTTTTTGGGTTTAACGCCGGGGGAGGgtt                                                                >  2:144719/1‑139 (MQ=255)
                                                                           tCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCGGTCATGACGGCGGTCACCAAAAATACGCATGCGTAACCGGGGATTTTCGGGGTTAACGCCTGGGGAGGgtttttttt                                                          >  1:82661/1‑137 (MQ=255)
                                                                           tCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCGCTCATGCAGGCGCGCACAAAAAATCCGCTTGCGTCCCCGGGGATTTTCGGGATTTACGCCGGGGCAGGttttttttt                                                          >  3:13958/1‑137 (MQ=255)
                                                                                  tCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCAAGGAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGCTTAACGCCGGTGCAGGGGTTTTTAtttgttcc                                                   >  7:264132/1‑133 (MQ=255)
                                                                                    cgcgcgCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGc                                                                            <  6:181859/112‑1 (MQ=255)
                                                                                    cgcgcgCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGc                                                                            >  5:181859/1‑112 (MQ=255)
                                                                                              ggTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGttttt                                        <  6:84029/138‑1 (MQ=255)
                                                                                               gTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTc                                      <  4:165297/139‑1 (MQ=255)
                                                                                                   cGACAGTCGTCGGCGCCTCGCTTGCCGCCCCGGGCGCGGTCAAGCCCGCGCTGACAAAAAATACCCATGCGTCAACGGGGATTTTATGGATTTACGCCGGTGCCGGTTTTTTTATTGTTTGCGGTTAGGTTTTTCttga                                  >  4:139823/1‑136 (MQ=255)
                                                                                                    gACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCat                                 <  6:83043/139‑1 (MQ=255)
                                                                                                       aGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATAtt                              <  5:284165/139‑1 (MQ=255)
                                                                                                         tcgtcgGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGc                                                                            <  4:382034/91‑1 (MQ=255)
                                                                                                         tcgtcgGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGc                                                                            >  3:382034/1‑91 (MQ=255)
                                                                                                         tcgtcgGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATAtt                              >  3:135974/1‑137 (MQ=255)
                                                                                                           gtcgGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACGGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGt                          <  5:247139/139‑1 (MQ=255)
                                                                                                           gtcgGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGt                          <  1:109724/139‑1 (MQ=255)
                                                                                                           gtcgGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGt                          <  4:135974/139‑1 (MQ=255)
                                                                                                           gtcgGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGt                          <  3:148204/139‑1 (MQ=255)
                                                                                                           gtcgGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGt                          <  3:148588/139‑1 (MQ=255)
                                                                                                           gtcgGCGCCTCGCTTGCAGCCGCGGGCGCGCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGt                          <  2:364477/139‑1 (MQ=255)
                                                                                                           ggcgGCGCCTCGCGTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGt                          <  4:13958/137‑1 (MQ=255)
                                                                                                                cgcCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCt                     <  6:199952/139‑1 (MQ=255)
                                                                                                                  ccTCGCGTGCAGCCGCGGGCGCGCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCa                   <  8:198564/139‑1 (MQ=255)
                                                                                                                      gCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAgg                 <  4:355672/137‑1 (MQ=255)
                                                                                                                      gCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAgg                 <  1:65314/137‑1 (MQ=255)
                                                                                                                      gCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAgg                 <  8:195165/137‑1 (MQ=255)
                                                                                                                      gCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAGGCa               <  4:386300/139‑1 (MQ=255)
                                                                                                                       ctggCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGGTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAGGCAt              <  8:224907/136‑1 (MQ=255)
                                                                                                                       cTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAGGCAt              <  7:327316/139‑1 (MQ=255)
                                                                                                                       cTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAGGCAt              <  5:31221/139‑1 (MQ=255)
                                                                                                                                   ggCGCTCTCCTGCAGGAGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTGTTTTTATTGTGGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAGGCATTAGTATGGAGCt  <  3:240102/139‑1 (MQ=255)
                                                                                                                                   ggCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAGGCATTAGTATGGAGCt  <  7:167448/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
GTACGTCATGGAGACAGGTCTATCAGGCGTTTACTTCCTTCAATGGAACAAATGAGCATGTCATCATCAAGGATGTCAGGCTTCCGCGCGCTTTGGTTGCGACAGTCGTCGGCGCCTCGCTTGCAGCCGCGGGCGCTCTCATGCAGGCGCTCACAAAAAATCCGCTTGCGTCACCGGGGATTTTCGGGATTAACGCCGGTGCAGGTTTTTTTATTGTTGCCGGATCGTTTTTTCTGCATATTCAGTCACCTCAGGCATTAGTATGGAGCT  >  CP006881/883608‑883877

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