Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,500 Δ1 bp 16.2% intergenic (‑176/+171) yjhU ← / ← yjhF putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,520,5000A.16.2% 113.6 / 17.5 37intergenic (‑176/+171)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base A (18/13);  new base . (3/3);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

ATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGG  >  NC_000913/4520359‑4520647
                                                                                                                                             |                                                                                                                                                    
aTCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTc                                                                                                                                              <  4:604134/149‑1 (MQ=255)
     aaCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCt                                                                                                                                         <  3:246492/149‑1 (MQ=255)
      aCCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                                                                        <  4:528976/149‑1 (MQ=255)
          tGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACa                                                                                                                                    >  2:562963/1‑149 (MQ=255)
                         ttGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTa                                                                                                                                 <  4:611561/137‑1 (MQ=255)
                         ttGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTa                                                                                                                                 >  3:611561/1‑137 (MQ=255)
                             tGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGaa                                                                                                                       >  4:319798/1‑143 (MQ=255)
                             tGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGaa                                                                                                                       <  3:319798/143‑1 (MQ=255)
                              gCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAccc                                                                                                                >  3:269144/1‑149 (MQ=255)
                                 gATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                                                                      >  4:97644/1‑124 (MQ=255)
                                 gATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTATAAAGCCGGAGACATGCTCCGGCTTTa                                                                                                                                      <  3:97644/124‑1 (MQ=255)
                                 gATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCAATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                                                                        >  2:455579/1‑122 (MQ=39)
                                 gATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCAATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                                                                        <  1:455579/122‑1 (MQ=39)
                                   ttGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCaa                                                                                                           >  2:70206/1‑149 (MQ=255)
                                      ccAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                                                                        <  1:562963/149‑1 (MQ=255)
                                                     tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCaa                                                                                         <  4:269144/149‑1 (MQ=255)
                                                               ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                                                                          >  1:377437/1‑122 (MQ=255)
                                                               ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACACATCAAg                                                                                                          <  2:377437/122‑1 (MQ=255)
                                                                             aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                                                                 >  1:379562/1‑149 (MQ=255)
                                                                                      aCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCa                                                        >  3:26258/1‑149 (MQ=255)
                                                                                       ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc                                                       >  1:221137/1‑149 (MQ=255)
                                                                                        cAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGc                                                                     <  3:102557/134‑1 (MQ=255)
                                                                                        cAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGc                                                                     >  4:102557/1‑134 (MQ=255)
                                                                                                  tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACgg                                            >  2:469959/1‑149 (MQ=255)
                                                                                                          aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGtc                                    >  4:447385/1‑149 (MQ=255)
                                                                                                                        ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTAt                      >  4:402722/1‑149 (MQ=255)
                                                                                                                         cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATg                     >  2:357958/1‑149 (MQ=255)
                                                                                                                          ccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATgg                    >  1:661753/1‑149 (MQ=255)
                                                                                                                              aCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGa                <  2:221137/149‑1 (MQ=255)
                                                                                                                                tGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACa              <  1:469959/149‑1 (MQ=255)
                                                                                                                                  tAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGt            <  3:402722/149‑1 (MQ=255)
                                                                                                                                     aGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa         >  4:444217/1‑149 (MQ=255)
                                                                                                                                     aGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaa         >  1:488065/1‑149 (MQ=255)
                                                                                                                                      gCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGTGACAGTGaac        >  3:320886/1‑149 (MQ=255)
                                                                                                                                        cGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGaacaa      <  1:357958/149‑1 (MQ=255)
                                                                                                                                         ggAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGcc                                       >  3:151547/1‑115 (MQ=255)
                                                                                                                                         ggAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGcc                                       <  4:151547/115‑1 (MQ=255)
                                                                                                                                            gaCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAgg  <  2:379562/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
ATCTTAACCCTGTTTACCCGGAGAGTTGATGCAGATTGCCAGTGGCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGACGGTTGCCGTCACTCACTAAATTATGGGGGACAGTGAACAACAGG  >  NC_000913/4520359‑4520647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: