Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,474,154 A→G 14.7% intergenic (‑59/+121) 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,1540AG14.7% 72.6 / 6.1 34intergenic (‑59/+121)PP_2167/talPhoD‑like phosphatase/transaldolase
Reads supporting (aligned to +/- strand):  ref base A (13/16);  new base G (3/2);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

GTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCA  >  NC_002947/2474007‑2474302
                                                                                                                                                    |                                                                                                                                                    
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)
                                                                                                                                            aTCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtttt                   <  1:298703/139‑1 (MQ=255)
                                                                                                                                            aTCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATc           <  7:347070/147‑1 (MQ=255)
                                                                                                                                            aTCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATc           >  8:347070/1‑147 (MQ=255)
                                                                                                                                            aTCGATCAGGGGAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAACCAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCGAACACGAACTGTATCAGGCCtgtg                   >  2:298703/1‑136 (MQ=255)
                                                                                                                                             tCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAAc        >  8:310638/1‑149 (MQ=255)
                                                                                                                                             atcgtccAGGGTCGGAGCGGCCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTACGGCGTGGCAAACACACCGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTCAGCTATCAAc        <  6:384290/143‑1 (MQ=255)
                                                                                                                                                    aGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCa  <  5:390731/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCA  >  NC_002947/2474007‑2474302

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: