Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,157,589 C→T 7.9% D343D (GAC→GAT)  PP_5400 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,157,5890CT 7.9% 148.8 / 2.4 55D343D (GAC→GAT) PP_5400hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base T (2/2);  total (28/27)
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.23e-01

TGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCTTTTTACACGAAGT  >  NC_002947/6157442‑6157729
                                                                                                                                                   |                                                                                                                                            
tgtgGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGaca                                                                                                                                             <  2:152128/149‑1 (MQ=255)
  tgGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGacac                                                                                                                                            <  8:139720/148‑1 (MQ=255)
    gTACTGAATTTTGCCAAGCACAATGATGCCTTGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGTGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCAt                                                                                                                                         <  3:114842/149‑1 (MQ=255)
     tACTGATTTTTGCCAAGAACAATGATGCGTCGGCGCTAGGAATGGCGCCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTCGGCCGGATGTTCTCCCTGCCGACACAATc                                                                                                                                        >  6:283493/1‑149 (MQ=255)
       cTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGa                                                                                                                                      >  5:91471/1‑149 (MQ=255)
        tGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTTACGACACCATCGAt                                                                                                                                     >  3:87439/1‑149 (MQ=255)
        tGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGAt                                                                                                                                     >  5:320591/1‑149 (MQ=255)
        tGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGAt                                                                                                                                     >  1:279110/1‑149 (MQ=255)
                  cAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCTTAGTCAACGATCTGATCTGGCGTTGTATTATGCCGGATGTTCTCACTGAAGAAACCATCGATCTCCTCTAcc                                                                                                                           >  3:225399/1‑149 (MQ=255)
                  cAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCATGACGCCACCAACGATCACCGCTAcc                                                                                                                           >  7:228100/1‑149 (MQ=255)
                        caATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTATTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTgg                                                                                                                     <  3:122779/149‑1 (MQ=255)
                        caATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTgg                                                                                                                     >  6:339878/1‑149 (MQ=255)
                           tgatgCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCgtgt                                                                                                                       >  3:279802/1‑144 (MQ=255)
                           tgatgCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCgtgt                                                                                                                       <  4:279802/144‑1 (MQ=255)
                             ctgtcTGGGAGGTAGGTTTGGAGGTGATGAAACTTCTACAGGCGGGTGATAATAGAGTTGTGATCTTAGTGAACGATATGCTATTGGGTTGTTTTTGGGAGGATGTTGTCCCTGACGAGTCCATGGTTGTGCTTTACCGTGTGGGTcta                                                                                                                <  8:228100/145‑3 (MQ=255)
                                ccTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGcc                                                                                                             >  4:126450/1‑149 (MQ=255)
                                     gcgcTAGGAATTGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATa                                                                                                        <  5:339878/149‑1 (MQ=255)
                                     gcgcTAGGAATGTCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCTTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGTATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATa                                                                                                        <  5:136842/149‑1 (MQ=255)
                                      cgcTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAg                                                                                                       <  5:40106/149‑1 (MQ=255)
                                          aGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGAtt                                                                                                   >  4:69523/1‑149 (MQ=255)
                                                       aTCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTg                                                                                      <  3:126450/149‑1 (MQ=255)
                                                           aaCTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAAt                                                                                  >  7:265904/1‑149 (MQ=255)
                                                                cTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGt                                                                               <  7:12472/147‑1 (MQ=255)
                                                                   cAGGCGTGTGAAAAAACATCCGGCATCGCAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGTATGTTCTCCCTTTCGAGACCATCGATCTGCTCTACCGTTTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATc                                                                          <  5:112217/149‑1 (MQ=255)
                                                                           tGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGgcg                                                                                         <  2:82898/126‑1 (MQ=255)
                                                                               aaaCCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGt                                                               >  7:313612/1‑148 (MQ=255)
                                                                                  ccAGCCGGCAGCGTAGTCAACGAGCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTAGCGTGTGGCTCCAGCAGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATc                                                           <  8:138960/149‑1 (MQ=255)
                                                                                   cAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATcc                                                          <  8:265904/149‑1 (MQ=255)
                                                                                     gCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCAt                                                                                                                                         <  4:89765/68‑1 (MQ=255)
                                                                                     gCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCAt                                                                                                                                         >  3:89765/1‑68 (MQ=255)
                                                                                       cGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTg                                                        >  4:205694/1‑147 (MQ=255)
                                                                                              ggaGTCAACGATCTGCTCTGGCGTGGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTTCCGTGGGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATCACGACGTGATCCTGAAGCTGGcc                                               <  4:87439/147‑1 (MQ=255)
                                                                                                 gTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAg                                            >  8:108090/1‑149 (MQ=255)
                                                                                                    aaCGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGa                                                                                     >  5:193463/1‑105 (MQ=255)
                                                                                                    aaCGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGa                                                                                     <  6:193463/105‑1 (MQ=255)
                                                                                                           tGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGc                                  <  3:69523/149‑1 (MQ=255)
                                                                                                             ctctGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTg                                <  2:279110/149‑1 (MQ=255)
                                                                                                              tctGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTg                                >  7:4496/1‑148 (MQ=255)
                                                                                                                  gCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGa                              <  6:91471/146‑1 (MQ=255)
                                                                                                                     ttGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCt                        >  4:240520/1‑149 (MQ=255)
                                                                                                                        tATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTc                      >  1:344742/1‑148 (MQ=255)
                                                                                                                          ttAGGCCGGATGTTCTCCCTGACGACACCATCGATctcc                                                                                                                                 <  3:330416/39‑1 (MQ=255)
                                                                                                                          ttAGGCCGGATGTTCTCCCTGACGACACCATCGATctcc                                                                                                                                 >  4:330416/1‑39 (MQ=255)
                                                                                                                            aGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTTAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCaaa                  <  3:240520/148‑1 (MQ=255)
                                                                                                                              gCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCt               >  7:16068/1‑149 (MQ=255)
                                                                                                                               ccGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCtt              <  6:320591/149‑1 (MQ=255)
                                                                                                                                         tcCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCTTTTTACACg      >  6:66961/1‑147 (MQ=255)
                                                                                                                                         tcCCTGACGACACCATCGATCTCCTCTACCAAGTGGCACCAGCCGAATAGGATTTCGCATCGCGCAGAAATTTTACATAGATAACGACGAGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCAAAATCATATTACACCaa    >  4:103516/1‑149 (MQ=255)
                                                                                                                                           ccTGACGATACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCttt             <  5:201221/138‑1 (MQ=255)
                                                                                                                                           ccTGACGATACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCttt             >  6:205983/1‑138 (MQ=255)
                                                                                                                                           ccTGACGATACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCttt             >  6:201221/1‑138 (MQ=255)
                                                                                                                                           ccTGACGATACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCttt             <  5:205983/138‑1 (MQ=255)
                                                                                                                                           ccTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGa                              <  7:70538/121‑1 (MQ=255)
                                                                                                                                           ccTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGa                              >  8:70538/1‑121 (MQ=255)
                                                                                                                                           ccTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCTTTTTACACGAAg   <  3:205694/148‑1 (MQ=255)
                                                                                                                                           ccTGAAGACACCATCGATCTCCTCTACCGTGTTGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACTACGTTATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCTTTTTACACGAAGt  <  3:103516/149‑1 (MQ=255)
                                                                                                                                            cTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCTTTTTACAGGAAGt  >  6:80909/1‑148 (MQ=255)
                                                                                                                                                 gacACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATc                                                           <  6:314060/86‑1 (MQ=255)
                                                                                                                                                 gacACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATc                                                           >  5:314060/1‑86 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
TGTGGTACTGAATTTTGCCAAGCACAATGATGCCTGGGCGCTAGGAATGGCGGCGATCAAACTTCTACAGGCGGGTGAAAAACCAGCCGGCATCGTAGTCAACGATCTGCTCTGGCGTTGTATTAGGCCGGATGTTCTCCCTGACGACACCATCGATCTCCTCTACCGTGTGGCTCCAGCCGAATAGGATTTCGCATGGCGCTGAAATTGTTCATCGATAACGACGTGATCCTGAAGCTGGCCCAGTACGGAATGCTGGAGCATCTTCCAAAGCTTTTTACACGAAGT  >  NC_002947/6157442‑6157729

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: ‑