Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,474,143 T→C 16.2% intergenic (‑48/+132) 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,1430TC16.2% 61.7 / 6.9 32intergenic (‑48/+132)PP_2167/talPhoD‑like phosphatase/transaldolase
Reads supporting (aligned to +/- strand):  ref base T (7/18);  new base C (3/2);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

GTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTA  >  NC_002947/2473999‑2474291
                                                                                                                                                |                                                                                                                                                     
gTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACa                                                                                                                                                   <  4:75774/149‑1 (MQ=255)
     gttgTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGGGGCGCAGTAGCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGGCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGAt                                                                                                                                             <  4:221074/146‑1 (MQ=255)
     gTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGAt                                                                                                                                             <  6:34160/149‑1 (MQ=255)
     gTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCGATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGAt                                                                                                                                             <  6:225390/149‑1 (MQ=255)
        gTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCCCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCa                                                                                                                                           <  2:288147/148‑1 (MQ=255)
             gAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAg                                                                                                                                     <  7:43744/149‑1 (MQ=255)
                 cagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGc                                                                                                                                 <  1:294043/149‑1 (MQ=255)
                  agacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCg                                                                                                                                <  4:378701/149‑1 (MQ=37)
                   gacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCgg                                                                                                                               <  3:126041/149‑1 (MQ=255)
                     cggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCATCGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAgggcgggc                                                                                                                             <  1:83511/149‑1 (MQ=37)
                                                             ggCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAg                                                                                     <  1:345480/149‑1 (MQ=255)
                                                             ggCGAGAACGAGGGGCAACGGGTCAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCTATCAGTGTAGCAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAg                                                                                     <  4:362800/149‑1 (MQ=255)
                                                                         gggCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAAGCCCGCCCGTTCCTCGATGAGTGTAGGCGCGGGCTTGGCCGCGGACAGGCCCGCGAATTGAATACACAACTAAGGCGGGGgaaac                                                                          >  1:109178/1‑148 (MQ=255)
                                                                              aCGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACg                                                                    <  3:188199/149‑1 (MQ=255)
                                                                                gggTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACggg                                                                  <  1:20165/149‑1 (MQ=255)
                                                                                           ggTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACt                                                       <  1:309360/149‑1 (MQ=255)
                                                                                                      gCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTccc                                            <  6:323823/149‑1 (MQ=255)
                                                                                                       cTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATAATGAACACTTcccg                                           >  5:193981/1‑149 (MQ=255)
                                                                                                             tGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTc                                     <  7:103302/149‑1 (MQ=255)
                                                                                                                    ggcggcCGGTTCGCGGGCACGCCCCCCATAACATCTAACAGTGTAGCCCCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCaaa                              <  1:180739/149‑1 (MQ=255)
                                                                                                                      cggcCTGTTCGCGGGCAAGCCCGCCCGTACATCGCTCAGTGTAGGAGCGGGCTTGCCCGCGAACCGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAACCAAAACGGGCGGGAAAGACTGACCACTTCCCGCCCGGCAGGCAAaca                            >  5:384290/1‑149 (MQ=38)
                                                                                                                        gcCTGTTCGCGGGCAAGCCCGCCGTTACAGAGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAAGACTTCCAGCCCGTCAGGCAAaaac                           >  6:317729/1‑148 (MQ=38)
                                                                                                                          cTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGaa                        <  6:5215/149‑1 (MQ=255)
                                                                                                                          cTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGGGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGa                         >  5:5215/1‑148 (MQ=255)
                                                                                                                           tGTTCGCGGGCAAGCCCGCCCTTACATCGAGTCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGaa                        >  1:233111/1‑149 (MQ=255)
                                                                                                                             ttCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGTAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACt                      <  7:218647/148‑1 (MQ=255)
                                                                                                                              tCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCATACACGAACTGt                    >  4:253487/1‑149 (MQ=255)
                                                                                                                                  gggCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACATACCCGAGAATTCAATACACAACTAAGGCGGGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTAt                  >  1:296640/1‑147 (MQ=255)
                                                                                                                                      aaGCCCGCCCTTACATCGATCAGTGTAGGCGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAACCGGGCGGAAAATACTGCACACTTCCCGCCCGTCTGGCACACCCGAACTGTATCAGGcc            >  4:4497/1‑149 (MQ=255)
                                                                                                                                        gcccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtt          >  2:209458/1‑149 (MQ=37)
                                                                                                                                          ccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtttt        >  1:246719/1‑149 (MQ=37)
                                                                                                                                          ccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCttt         >  3:213482/1‑148 (MQ=37)
                                                                                                                                                ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTa  >  7:352824/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GTTGGGTAGTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTA  >  NC_002947/2473999‑2474291

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: