Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,534,897 A→T 15.2% intergenic (+28/+68) nhoA → / ← yddE arylamine N‑acetyltransferase/PF02567 family protein YddE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,534,8970AT15.2% 111.7 / 13.0 46intergenic (+28/+68)nhoA/yddEarylamine N‑acetyltransferase/PF02567 family protein YddE
Reads supporting (aligned to +/- strand):  ref base A (25/14);  new base T (3/4);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

TGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGA  >  NC_000913/1534752‑1535046
                                                                                                                                                 |                                                                                                                                                     
tGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACg                                                                                                                                                   >  1:454522/1‑150 (MQ=255)
        tGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                           >  2:58865/1‑150 (MQ=255)
          cAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                         >  2:1141620/1‑150 (MQ=255)
           aaGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                        <  1:1152447/150‑1 (MQ=255)
             gAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCAt                                                                                                                                      <  2:375806/150‑1 (MQ=255)
              aaCAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATAAGGCGTTTACGCCGCATCCGGCaac                                                                                                                                     <  1:697466/150‑3 (MQ=255)
               aCAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGACGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATAAGGCGTTTACGCCGCATCCGGCaacc                                                                                                                                    <  2:216148/150‑4 (MQ=255)
                 aGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATAAGGCGTTTACGCCGCATCCGGCaacctg                                                                                                                                  <  1:1332340/150‑6 (MQ=38)
                                 ggATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCTGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGAc                                                                                                                  >  2:606821/1‑150 (MQ=255)
                                 ggATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGAc                                                                                                                  >  2:342150/1‑150 (MQ=255)
                                   atgatgCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGAcat                                                                                                                >  2:1090711/1‑148 (MQ=255)
                                       tgCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACa                                                                                                            <  2:144514/150‑1 (MQ=255)
                                        gCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAg                                                                                                           >  1:830386/1‑150 (MQ=255)
                                          gAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATACGGCATACATTAGCCCGGTTGTAGACCTGACAGGc                                                                                                         >  2:1149492/1‑150 (MQ=255)
                                                     ttACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCa                                                                                              >  1:1409497/1‑150 (MQ=255)
                                                      tACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCTGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAg                                                                                             >  1:1008144/1‑150 (MQ=255)
                                                        ccgtttaTGATTTAGCGCTGGTTATGTCGGCGTTTGATACTCAACAGGAGGCGGGAAAATAATTTATGTCAGGGTGCCGGATGCTTCGTAAACGCCTTATACGGCATACATTAGCACGGTTTTAGAACTTACATGCTAAGCGCATCAGAc                                                                                           <  2:1037749/144‑1 (MQ=255)
                                                        ccGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGAc                                                                                           >  1:502524/1‑150 (MQ=255)
                                                          gTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACaa                                                                                         >  1:806513/1‑150 (MQ=255)
                                                            ggATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACg                                                                                       >  2:1180197/1‑150 (MQ=255)
                                                             gATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACgg                                                                                      >  2:244801/1‑150 (MQ=255)
                                                               tGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGAGAACGGtt                                                                                    <  2:1311462/150‑1 (MQ=255)
                                                                     aGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCagat                                                                              <  2:358252/150‑2 (MQ=255)
                                                                         cTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCa                                                                          >  1:491887/1‑150 (MQ=255)
                                                                           ggTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAAt                                                                        <  1:1255627/150‑1 (MQ=255)
                                                                                 ggcggcGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCa                                                                  >  1:457067/1‑150 (MQ=255)
                                                                                  gcggcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGa                                                                                            >  1:851393/1‑123 (MQ=255)
                                                                                  gcggcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGa                                                                                            <  2:851393/123‑1 (MQ=255)
                                                                                     gcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCg                                                              <  2:454522/150‑1 (MQ=255)
                                                                                            aTACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGaa                                                       >  1:634080/1‑150 (MQ=255)
                                                                                               cGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCa                                                           <  1:251575/143‑1 (MQ=255)
                                                                                               cGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCa                                                           >  2:251575/1‑143 (MQ=255)
                                                                                                 cACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaa                                                  >  2:292574/1‑150 (MQ=255)
                                                                                                           cGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCtg                                        <  1:68114/150‑1 (MQ=255)
                                                                                                                     atTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAAt                              >  1:904705/1‑150 (MQ=255)
                                                                                                                       ttATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATgg                            <  1:1141620/150‑1 (MQ=255)
                                                                                                                             cAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACtt                      >  2:1381548/1‑150 (MQ=255)
                                                                                                                                   gCCGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCg                >  1:558087/1‑150 (MQ=255)
                                                                                                                                    ccGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCgg               >  1:691094/1‑150 (MQ=255)
                                                                                                                                      ggataaggcgTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCggtt             >  1:585658/7‑150 (MQ=255)
                                                                                                                                      ggATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCggtt             <  2:52042/150‑1 (MQ=255)
                                                                                                                                      ggATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCggtt             >  2:819303/1‑150 (MQ=255)
                                                                                                                                           cggcgTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTa        >  2:699083/1‑150 (MQ=255)
                                                                                                                                             gcgTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATc      <  2:1008144/150‑1 (MQ=255)
                                                                                                                                              cgTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCa     >  1:1382952/1‑150 (MQ=255)
                                                                                                                                               gTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCAc    <  1:1180197/150‑1 (MQ=255)
                                                                                                                                                 aaaCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGa  >  2:172400/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGA  >  NC_000913/1534752‑1535046

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