Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,474,162 A→G 26.3% intergenic (‑67/+113) PP_2167 ← / ← tal PhoD‑like phosphatase/transaldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,474,1620AG26.3% 58.9 / 20.8 39intergenic (‑67/+113)PP_2167/talPhoD‑like phosphatase/transaldolase
Reads supporting (aligned to +/- strand):  ref base A (17/11);  new base G (5/5);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

GACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCTTTT  >  NC_002947/2474018‑2474309
                                                                                                                                                |                                                                                                                                                   
gacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCgg                                                                                                                                                 <  6:339207/149‑1 (MQ=255)
  cggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGTTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCggg                                                                                                                                                <  4:197276/148‑1 (MQ=37)
      cagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGTGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCTCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGc                                                                                                                                           <  5:133697/149‑1 (MQ=255)
      cagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCGGGCTTGc                                                                                                                                           <  6:349612/149‑1 (MQ=37)
        gacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGccc                                                                                                                                         <  7:101034/149‑1 (MQ=255)
             gggggTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTGGCCCGCGaa                                                                                                                                    <  5:321776/148‑1 (MQ=255)
              ggggTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCg                                                                                                                                                  >  3:151697/1‑134 (MQ=255)
              ggggTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCg                                                                                                                                                  <  4:151697/134‑1 (MQ=255)
                  ttCCAGGCGGCGCAGAACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGCTGGCTATGGCGGCCTATTCGCGGCCAAGCCAGCCCTTACATCGATCAGTGTAGGACCGGGCTTGCCCGCGAACAGGc                                                                                                                               >  2:281546/1‑149 (MQ=255)
                         cggcgAAGTACCGGGCCTGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCGGGCTTGCCcg                                                                                                                                        <  8:350776/133‑1 (MQ=21)
                          ggcgCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCtt                                                                                                                                             <  3:212624/127‑1 (MQ=255)
                          ggcgCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCtt                                                                                                                                             >  4:212624/1‑127 (MQ=255)
                                               gaACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCGGGCTTGc                                                                                                                                           <  1:259522/108‑1 (MQ=18)
                                               gaACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCGGGCTTGc                                                                                                                                           >  2:259522/1‑108 (MQ=18)
                                                     ggggCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAAtt                                                                                                                      <  3:233932/123‑1 (MQ=255)
                                                     ggggCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCAGTTAGCGGGCCAGCCCCCCCATACATCGATCACTTTAGGAGCGGGCGTGGCCCCGAAACGGCCCCACAAtt                                                                                                                      >  4:233932/1‑123 (MQ=255)
                                                                                 cAGCTGGGGTGGGTTTGGCGGCCTGTGCGCGGGCATGCCCGCCCTTACATCGATCAGTTTAGCTGGGGGCTTGCCCGCGAACAGGCCCGAGATTTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTc                                                                <  4:379436/149‑1 (MQ=255)
                                                                                                         gTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACt                                        >  5:174146/1‑149 (MQ=255)
                                                                                                             gcgGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTa                                     >  3:377581/1‑148 (MQ=255)
                                                                                                                   aaGCCCGCTCCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGcc                              >  2:2524/1‑149 (MQ=37)
                                                                                                                     gcccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtt                            >  1:159419/1‑149 (MQ=37)
                                                                                                                     gcccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtt                            >  1:153701/1‑149 (MQ=37)
                                                                                                                      cccgcccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCttt                           >  5:260410/1‑149 (MQ=255)
                                                                                                                       ccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGACAGGCAAACACGAACTGTATCAGGCCttt                           >  7:89724/1‑148 (MQ=37)
                                                                                                                       ccgcccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTa                                     >  6:382472/1‑138 (MQ=255)
                                                                                                                       ccgcccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTa                                     <  5:382472/138‑1 (MQ=255)
                                                                                                                       ccgcccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCttt                           >  4:228020/1‑148 (MQ=255)
                                                                                                                         gcccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtttt                          >  5:46425/1‑147 (MQ=255)
                                                                                                                             ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTa                    >  5:182110/1‑149 (MQ=255)
                                                                                                                             ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTa                    >  8:159181/1‑149 (MQ=255)
                                                                                                                             ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTa                    >  5:294231/1‑149 (MQ=255)
                                                                                                                                cATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCTGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATc                  <  4:377581/148‑1 (MQ=255)
                                                                                                                                cATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATc                  <  3:228020/148‑1 (MQ=255)
                                                                                                                                cATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATc                  >  7:326876/1‑148 (MQ=255)
                                                                                                                                cATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATc                  >  7:77931/1‑148 (MQ=255)
                                                                                                                                 aTCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCaa                >  1:306717/1‑149 (MQ=255)
                                                                                                                                  tCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAAc               <  6:294231/149‑1 (MQ=255)
                                                                                                                                        gatgtAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCa         <  8:89724/147‑1 (MQ=255)
                                                                                                                                               gAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCGTTTCAGCTATCAACTGCGCAAGCtttt  >  8:333248/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCTTTT  >  NC_002947/2474018‑2474309

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