Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,880,851 T→A 13.6% intergenic (‑92/+35) leuE ← / ← dmlR leucine efflux protein/transcriptional activator of dmlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,880,8510TA13.6% 127.5 / 12.2 44intergenic (‑92/+35)leuE/dmlRleucine efflux protein/transcriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base A (2/4);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCA  >  NC_000913/1880710‑1880999
                                                                                                                                             |                                                                                                                                                    
gCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTg                                                                                                                                               >  2:185501/1‑149 (MQ=255)
 ccccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGa                                                                                                                                              <  4:287875/149‑1 (MQ=255)
  cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGat                                                                                                                                             <  2:373144/149‑1 (MQ=255)
  cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGat                                                                                                                                             <  1:108623/149‑1 (MQ=255)
   ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGatt                                                                                                                                            <  1:38966/149‑2 (MQ=255)
   ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGa                                                                                                                                              >  2:92132/1‑147 (MQ=255)
   ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGa                                                                                                                                              <  1:92132/147‑1 (MQ=255)
       ccAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCaa                                                                                                                                        <  4:445008/149‑1 (MQ=255)
        cAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAc                                                                                                                                       <  4:342660/149‑1 (MQ=255)
        cAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAc                                                                                                                                       <  3:247954/149‑1 (MQ=255)
        cAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAc                                                                                                                                       <  2:67755/149‑1 (MQ=255)
         aGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAcc                                                                                                                                      >  1:108369/1‑149 (MQ=255)
                   aGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGa                                                                                                                            <  1:281141/149‑1 (MQ=255)
                   aGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGa                                                                                                                            <  2:108369/149‑1 (MQ=255)
                     tAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGACTATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAAt                                                                                                                          <  2:367094/149‑1 (MQ=255)
                     tAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAAt                                                                                                                          <  4:377182/149‑1 (MQ=255)
                           aGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTaa                                                                                                                    <  3:123279/149‑1 (MQ=255)
                             aaCCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATc                                                                                                                  >  1:266296/1‑149 (MQ=255)
                               ccccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAg                                                                                                                <  3:136383/149‑1 (MQ=255)
                                cccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAgg                                                                                                               >  1:116723/1‑149 (MQ=255)
                                                  aTTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAAt                                                                                             >  1:131020/1‑149 (MQ=255)
                                                      aCCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCgg                                                                                         <  2:221608/149‑1 (MQ=255)
                                                         tctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTa                                                                                      >  3:68108/1‑149 (MQ=255)
                                                           tttAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGt                                                                                    >  1:378528/1‑149 (MQ=255)
                                                                ttATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCg                                                                               >  2:199949/1‑149 (MQ=255)
                                                                     tttCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCaaa                                                                          >  4:126677/1‑149 (MQ=255)
                                                                          aTCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACt                                                                     <  4:158929/149‑1 (MQ=255)
                                                                          aTCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACt                                                                     >  1:462430/1‑149 (MQ=255)
                                                                                   gATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTaaaa                                                            >  2:11082/1‑149 (MQ=255)
                                                                                     ttAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAAt                                                          <  2:462430/149‑1 (MQ=255)
                                                                                         aGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTa                                                      <  3:126677/149‑1 (MQ=255)
                                                                                         aGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTa                                                      >  3:125311/1‑149 (MQ=255)
                                                                                          ggTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTAc                                                     <  1:11082/149‑1 (MQ=255)
                                                                                              aaTATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGtt                                                 >  4:274583/1‑149 (MQ=255)
                                                                                                    cTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGc                                           >  3:394508/1‑149 (MQ=255)
                                                                                                     tAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTc                                                        <  3:16559/135‑1 (MQ=255)
                                                                                                     tAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTc                                                        >  4:16559/1‑135 (MQ=255)
                                                                                                      aaCCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCAc                                         <  2:378528/149‑1 (MQ=255)
                                                                                                       aCCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACt                                        >  2:360025/1‑149 (MQ=255)
                                                                                                         cAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACttt                                      >  4:428311/1‑149 (MQ=255)
                                                                                                          aaTCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTc                                     <  1:185501/149‑1 (MQ=255)
                                                                                                            tCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGc                                   >  4:224850/1‑149 (MQ=255)
                                                                                                                          gAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGcc                     <  2:381212/149‑1 (MQ=255)
                                                                                                                                      tCAGGTTAAAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCACTGACGTCGCCAGCCTTGAAACATaaa         >  3:454267/1‑149 (MQ=255)
                                                                                                                                           ttttAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGcc    <  2:116723/149‑1 (MQ=255)
                                                                                                                                             ttAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCa  <  1:360025/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCA  >  NC_000913/1880710‑1880999

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