Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,147,479 G→A *114* (TAG→TAA
T280I (ACT→ATT) 
frwD →
yijO ←
putative enzyme IIB component of PTS
AraC family putative transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,147,4790GA96.4% 82.3 / ‑5.8 28*114* (TAG→TAA
T280I (ACT→ATT) 
frwD
yijO
putative enzyme IIB component of PTS
AraC family putative transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/13);  minor base T (1/0);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGC  >  NC_000913/4147346‑4147597
                                                                                                                                     |                                                                                                                      
gTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGt                                                                                                                   <  2:574047/139‑1 (MQ=255)
      aGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCt                                                                                                             >  1:273777/1‑139 (MQ=255)
           tttGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAg                                                                                                        <  2:273777/139‑1 (MQ=255)
                aCATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGc                                                                                                   <  2:514203/139‑1 (MQ=255)
                              tgCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCgg                                                                                     <  1:461789/139‑1 (MQ=255)
                                             aCGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGc                                                                      >  2:346152/1‑139 (MQ=255)
                                                  ttCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGa                                                                    >  2:313659/1‑136 (MQ=255)
                                                  ttCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGa                                                                    <  1:313659/136‑1 (MQ=255)
                                                   tCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGt                                                                >  2:53661/1‑139 (MQ=255)
                                                             gCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGGGAATTGGTTTTTTCGGTCAGTTGCCTGTGGTATTCCCGGCGATATTCGGACGGCGAGCGTCGGGTGtttt                                                      >  1:563937/1‑139 (MQ=255)
                                                               cGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTgg                                                                                                    >  1:648040/1‑91 (MQ=255)
                                                               cGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTgg                                                                                                    <  2:648040/91‑1 (MQ=255)
                                                                 ccacGGGTTATGGACGCGGTTCGCAAACTTCTTTTCTCCCCCCAACAAACCCAACCTTTTCTTGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACg                                                  <  2:563937/136‑1 (MQ=255)
                                                                                gcgGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGaa                                   >  1:47939/1‑139 (MQ=255)
                                                                                    tGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGtt                                                                                                      <  2:7150/68‑1 (MQ=255)
                                                                                    tGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGtt                                                                                                      >  1:7150/1‑68 (MQ=255)
                                                                                          aaGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTg                         <  2:500/139‑1 (MQ=255)
                                                                                            gTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTc                       >  1:367113/1‑139 (MQ=255)
                                                                                              gCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGc                                       >  2:657656/1‑121 (MQ=255)
                                                                                              gCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGc                                       <  1:657656/121‑1 (MQ=255)
                                                                                              gCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAACTAGTTGCTGTCGa                     >  2:70897/1‑139 (MQ=255)
                                                                                                tttCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACa                   >  2:68393/1‑139 (MQ=255)
                                                                                                       gccgcAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGc            >  2:724286/1‑139 (MQ=255)
                                                                                                         cgcAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTc                                                                                                          <  1:204862/43‑1 (MQ=255)
                                                                                                         cgcAGCAAACCCATCTTATTCTGGAGTAATCGGTTTTTCTGTc                                                                                                          >  2:204862/1‑43 (MQ=255)
                                                                                                                 aCCCATCTTATTCTGGAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTgcgc  <  1:53661/139‑1 (MQ=255)
                                                                                                                                gAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCg             >  1:694963/1‑113 (MQ=255)
                                                                                                                                gAGTAATCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCg             <  2:694963/113‑1 (MQ=255)
                                                                                                                                     |                                                                                                                      
GTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCCGCAGCAAACCCATCTTATTCTGGAGTAGTCGGTTTTTCTGTCAGTTGGCTGTGGTACTGCCGGCGATATTCCGACGGCGAGCGTTCGGTGTTTTTACGAAACAGACGGCAGAAGTAGTTGCTGTCGACAAAGCCGCAGGCGTGCGC  >  NC_000913/4147346‑4147597

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