Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 411,489 T→A 9.3% I331I (ATA→ATT)  araJ ← L‑arabinose‑inducible putative transporter, MFS family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913411,4890TA 9.3% 142.3 / 5.9 43I331I (ATA→ATT) araJL‑arabinose‑inducible putative transporter, MFS family
Reads supporting (aligned to +/- strand):  ref base T (16/23);  new base A (2/2);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GCTTATAGCGACCATACAGCAGCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCAAAAATAAGCGACGTTGTTTTCATGCCGCCGCAGAAAAA  >  NC_000913/411344‑411635
                                                                                                                                                 |                                                                                                                                                  
gcTTATAGCGACCATACAGCAGCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTAtt                                                                                                                                                 <  4:25357/149‑1 (MQ=255)
    ataGCGACCATACAGCAGCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGcc                                                                                                                                             <  3:400966/149‑1 (MQ=255)
      aGCGACCATACAGCAGCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCa                                                                                                                                           >  3:15433/1‑149 (MQ=255)
      aGCGACCATACAGCAGCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCa                                                                                                                                           >  1:314066/1‑149 (MQ=255)
      aGCGACCATACAGCAGCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCa                                                                                                                                           >  1:75681/1‑149 (MQ=255)
      aGCGACCATACAGCAGCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCAATTTGCCCa                                                                                                                                           >  1:136879/1‑149 (MQ=255)
      aGCGACCATACAGCAGCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCAATTTGCCCa                                                                                                                                           >  3:263611/1‑149 (MQ=255)
       gCGACCATACAGCAGCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCAc                                                                                                                                          >  4:308720/1‑149 (MQ=255)
                cagcagCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCAc                                                                                                                                 <  2:469017/149‑1 (MQ=255)
                 agcagcAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCAcc                                                                                                                                <  1:133542/149‑1 (MQ=255)
                    agcaACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCtaa                                                                                                                             <  1:24261/149‑1 (MQ=255)
                     gcaACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCtaat                                                                                                                            >  1:165387/1‑149 (MQ=255)
                                aTCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGcc                                                                                                                 <  3:530773/149‑1 (MQ=255)
                                   gcagcaAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCttt                                                                                                              <  2:485750/149‑1 (MQ=255)
                                   gcagcaAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCttt                                                                                                              <  4:286469/149‑1 (MQ=255)
                                    cagcaAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTg                                                                                                             <  4:212351/149‑1 (MQ=255)
                                      gcaAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGc                                                                                                           >  2:458764/1‑149 (MQ=255)
                                                             cagCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTa                                                                                    <  2:289816/149‑1 (MQ=255)
                                                             cagCGCCACGTAATTATATGCCAGCCCCAGCGTCAGAATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCGCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTa                                                                                    <  3:502257/149‑1 (MQ=255)
                                                                  ccACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGt                                                                               <  2:22152/149‑1 (MQ=255)
                                                                         aTTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGaaa                                                                        <  3:81734/149‑1 (MQ=255)
                                                                            atatGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAggg                                                                     >  2:328443/1‑149 (MQ=255)
                                                                            atatGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAggg                                                                     >  2:67291/1‑149 (MQ=255)
                                                                                gCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCaaa                                                                 <  3:140557/149‑1 (MQ=255)
                                                                                gCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCaaa                                                                 >  3:357194/1‑149 (MQ=255)
                                                                                gCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTCGTAGCGGTGCTGAAAGGGCaaa                                                                 <  2:131401/149‑1 (MQ=255)
                                                                                  cAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAaata                                                               >  2:209235/1‑149 (MQ=255)
                                                                                  aaGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTAACGAGGTTAGACGCTATTTGCCCACCTGCGGCACCTAAAAACTGTCCGCCTTTGGCGTTTTGTAGAAACGATATTTGTAGCGGTGCTGAAAGGGCAaata                                                               >  4:349473/2‑149 (MQ=255)
                                                                                     ccccAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATc                                                            <  4:405197/149‑1 (MQ=255)
                                                                                                      tACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCAATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATaaaa                                           <  4:263611/149‑1 (MQ=255)
                                                                                                      tACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCAATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATaaaa                                           <  2:136879/149‑1 (MQ=255)
                                                                                                         cTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATc                                                            <  2:464471/129‑1 (MQ=255)
                                                                                                         cTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATc                                                            >  1:464471/1‑129 (MQ=255)
                                                                                                          tCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCaa                                       >  4:440194/1‑149 (MQ=255)
                                                                                                           ccGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCaaa                                      <  4:204468/149‑1 (MQ=255)
                                                                                                                     cgcCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCAAAAATAAGCGAc                            >  3:181380/1‑149 (MQ=255)
                                                                                                                           cGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCAAAAATAAGCGACgttgtt                      <  4:5851/149‑1 (MQ=255)
                                                                                                                              cgcTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCAAAAATAAGCGACGTTGTTTTc                   <  1:28383/149‑1 (MQ=255)
                                                                                                                               gcTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCAAAAATAAGCGACGTTGTTTTCa                  <  2:314066/149‑1 (MQ=255)
                                                                                                                                   ccGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCAAAAATAAGCGACGTTGTTTTCATgcc              >  4:267599/1‑149 (MQ=255)
                                                                                                                                            aaCGCTATTTGCCCACCTGCGGCGCCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCAAAAATAAGCGACGTTGTTTTCATGCCGCCGCAGaa     >  4:36664/1‑149 (MQ=255)
                                                                                                                                               gCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCAAAAATAAGCGACGTTGTTTTCATGCCGCCGCAGaaaaa  <  4:148358/149‑1 (MQ=255)
                                                                                                                                               gCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCAAAAATAAGCGACGTTGTTTTCATGCCGCCGCAGaaaaa  <  3:447864/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCTTATAGCGACCATACAGCAGCAACGACGACATCGCAGCAAACGAAAGCAGGGCGGCAGGCAGCGCCACGTAATTATATGCCAGCCCCAGCGTCAGCATCATACCTCCGCAATATGCGCCGACGGCGCTACCGAGGTTAAACGCTATTTGCCCACCTGCGGCACCTAATAACTCTCCGCCTTTGGCGTTTTGTAGTAACAATATTTGTAGCGGTGCTGAAAGGGCAAATAATCCCGCGCAACAAATAAAAGCAAAAATAAGCGACGTTGTTTTCATGCCGCCGCAGAAAAA  >  NC_000913/411344‑411635

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: ‑