Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,281,677 T→A 6.4% intergenic (+48/+106) yjcE → / ← yjcF putative cation/proton antiporter/pentapeptide repeats protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,281,6770TA 6.4% 216.5 / 3.9 62intergenic (+48/+106)yjcE/yjcFputative cation/proton antiporter/pentapeptide repeats protein
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base A (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATA  >  NC_000913/4281534‑4281822
                                                                                                                                               |                                                                                                                                                 
cTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtttt                                                                                                                                              <  1:1264270/149‑1 (MQ=255)
     cgcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCAcc                                                                                                                                         <  1:220521/149‑1 (MQ=255)
       ccACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCtt                                                                                                                                       >  2:310117/1‑149 (MQ=255)
       ccACACGGGAGATCAACAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCtt                                                                                                                                       <  1:310117/149‑1 (MQ=255)
          cacGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTaa                                                                                                                                    <  1:226937/149‑1 (MQ=255)
                   tCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCttt                                                                                                                                               >  1:1064576/1‑129 (MQ=255)
                   tCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCttt                                                                                                                                               <  2:1064576/129‑1 (MQ=255)
                      gCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCAc                                                                                                                                          >  2:478794/1‑131 (MQ=255)
                      gCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCAc                                                                                                                                          <  1:478794/131‑1 (MQ=255)
                                              tGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAAAAGATTGCACCttt                                                                                                                                      <  4:960157/111‑1 (MQ=255)
                                              tGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAAAAGATTGCACCttt                                                                                                                                      >  1:1242182/1‑111 (MQ=39)
                                               gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACa                                                                                               >  3:1210471/1‑149 (MQ=255)
                                                  cGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTg                                                                                                                                             >  3:1094989/1‑100 (MQ=255)
                                                  cGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTg                                                                                                                                             <  4:1094989/100‑1 (MQ=255)
                                                      ctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTaa                                                                                                   >  4:687483/1‑138 (MQ=255)
                                                      ctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTaa                                                                                                   <  3:687483/138‑1 (MQ=255)
                                                       tcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAgg                                                                                                                                >  2:455574/1‑108 (MQ=255)
                                                       tcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAgg                                                                                                                                <  1:455574/108‑1 (MQ=255)
                                                               cTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCAt                                                                                                           >  3:625607/1‑121 (MQ=255)
                                                               cTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCAt                                                                                                           <  4:625607/121‑1 (MQ=255)
                                                                    aGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCAt                                                                                                           <  1:524934/116‑1 (MQ=255)
                                                                    aGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGGGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCAt                                                                                                           >  2:524934/1‑116 (MQ=255)
                                                                           cTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtttt                                                                                                                                              >  1:1292727/1‑74 (MQ=255)
                                                                           cTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtttt                                                                                                                                              <  2:1292727/74‑1 (MQ=255)
                                                                            tAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCttt                                                                                                                                               >  4:1072714/1‑72 (MQ=255)
                                                                            tAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCttt                                                                                                                                               <  3:1072714/72‑1 (MQ=255)
                                                                                            gTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtttt                                                                                                                                              <  1:1282690/57‑1 (MQ=255)
                                                                                            gTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtttt                                                                                                                                              >  2:1282690/1‑57 (MQ=255)
                                                                                                  cGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttc                                                      <  1:1180810/139‑1 (MQ=255)
                                                                                                  cGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttc                                                      >  2:1180810/1‑139 (MQ=255)
                                                                                                   gAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTg                                                                                                                                             <  1:1186942/51‑1 (MQ=255)
                                                                                                   gAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTg                                                                                                                                             >  2:1186942/1‑51 (MQ=255)
                                                                                                            gTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCtt                                                                                                                                       <  3:1240342/48‑1 (MQ=255)
                                                                                                            gTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCAt                                                               >  4:1198427/1‑120 (MQ=255)
                                                                                                            gTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCAt                                                               <  3:1198427/120‑1 (MQ=255)
                                                                                                            gTAACCACGCCCTAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCtt                                                                                                                                       >  4:1240342/1‑48 (MQ=255)
                                                                                                                      ccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGTTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAAtt                                                                                       >  3:1008291/1‑86 (MQ=255)
                                                                                                                      ccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAAtt                                                                                       <  4:1008291/86‑1 (MQ=255)
                                                                                                                          ggTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCataata                    >  1:388443/1‑149 (MQ=255)
                                                                                                                             aaTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCtttt                                                                            >  1:1323858/1‑90 (MQ=255)
                                                                                                                             aaTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCtttt                                                                            <  2:1323858/90‑1 (MQ=255)
                                                                                                                              aTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTg                <  2:388443/149‑1 (MQ=255)
                                                                                                                               tCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCAt                                                               <  4:49587/101‑1 (MQ=255)
                                                                                                                               tCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCAt                                                               >  3:49587/1‑101 (MQ=255)
                                                                                                                                cTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCa              >  2:516905/1‑149 (MQ=255)
                                                                                                                                  gAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGt            >  1:853761/1‑149 (MQ=255)
                                                                                                                                  gAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGt            >  3:1123321/1‑149 (MQ=255)
                                                                                                                                    aaaaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTaaa                                                                         >  1:861467/1‑86 (MQ=255)
                                                                                                                                    aaaaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTaaa                                                                         <  2:861467/86‑1 (MQ=255)
                                                                                                                                     aaaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttcttc                                                   <  1:642947/107‑1 (MQ=255)
                                                                                                                                     aaaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttcttc                                                   <  1:123453/107‑1 (MQ=255)
                                                                                                                                     aaaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttcttc                                                   >  2:123453/1‑107 (MQ=255)
                                                                                                                                     aaaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttcttc                                                   >  2:642947/1‑107 (MQ=255)
                                                                                                                                      aaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTaaa                                                                         >  3:263783/1‑84 (MQ=38)
                                                                                                                                      aaGGTGCAAAAGATTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTaaa                                                                         <  4:263783/84‑1 (MQ=38)
                                                                                                                                          tGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCa    >  4:624275/1‑149 (MQ=255)
                                                                                                                                           gCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGAc                                                <  3:311753/104‑1 (MQ=255)
                                                                                                                                           gCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGAc                                                >  4:311753/1‑104 (MQ=255)
                                                                                                                                           gCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCAt   >  4:939311/1‑149 (MQ=255)
                                                                                                                                            cAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCAtg  >  1:1285244/1‑148 (MQ=255)
                                                                                                                                             aaTCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAc                                           <  1:913526/107‑1 (MQ=255)
                                                                                                                                             aaTCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAc                                           >  2:913526/1‑107 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAATCTTCATA  >  NC_000913/4281534‑4281822

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