Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,355,029 T→G 13.6% T151P (ACA→CCA)  U712_11925 ← Stage III sporulation protein AA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,355,0290TG13.6% 82.7 / 2.9 48T151P (ACA→CCA) U712_11925Stage III sporulation protein AA
Reads supporting (aligned to +/- strand):  ref base T (13/26);  new base G (5/2);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 9.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTCTGACCAAATTGATGCTGCGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGATATAAGTATGGAAGCAGCGGCTCCGCAATCCCCAGCTTTTGCCGGGCAATTCGTATATTAAAAGATGCGATGTCACGGAG  >  CP006881/2354892‑2355164
                                                                                                                                         |                                                                                                                                       
tctcTGACCAAATTGATGCTGCGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCCCAATGTTTTTTTTGCCGGTGGTGGAGAGCCGCCCAAGTGCCCGGAGAAGggtt                                                                                                                                        >  1:46787/1‑139 (MQ=255)
    tGACCAAATTGATGCTGCGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGttt                                                                                                                                    >  8:92113/1‑139 (MQ=255)
    tGACCAAATTGATGCTGCGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGGTGttt                                                                                                                                    >  1:313576/1‑139 (MQ=255)
      aCCCAATTTATGCTGCGGAATACCGCGAAGACAGCCGGCATTCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGGATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTg                                                                                                                                  <  4:193263/139‑1 (MQ=255)
      aCCAAATTGATGCTGCGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTg                                                                                                                                  <  4:219493/139‑1 (MQ=255)
      aCCAAATTGATGCTGCGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTg                                                                                                                                  <  8:113018/139‑1 (MQ=255)
      aCCAAATTGATGCTGCGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGGTGTTTTg                                                                                                                                  >  5:238165/1‑139 (MQ=255)
          aaatGATGCTGCGGAATACCGCGAAGACCTCCGGCAAACACTGACCGTTCATCGACGATCCAAGTATTTACGGGAAGCATGTTTTTTTTGCCGGTGCTTGAGCGCATCTCGAGGTACCCGAGAAGTGTTTTGTTGcggt                                                                                                                              >  7:224217/4‑136 (MQ=255)
                  cgccGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCCTCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTgggg                                                                                                                       <  6:75669/135‑1 (MQ=255)
                   gtcGGTATTCCGCGAATACCTCCGGCCATGTCTGACCGGTGATCGTCGATCCCTGTGTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCg                                                                                                                     <  2:46787/137‑1 (MQ=255)
                     cGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGTCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGGCGGTTTGggggggg                                                                                                                   >  3:193322/1‑138 (MQ=255)
                     cGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGt                                                                                                                   >  3:74909/1‑139 (MQ=255)
                              gcgAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATc                                                                                                          <  8:105778/139‑1 (MQ=255)
                              gcgAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATc                                                                                                          <  6:119241/139‑1 (MQ=255)
                              gcgAAGACATCCCGGAATCTCAGACCGTGCATCGACGCTCCCTGTCTCTCCGGCCCGCATGTTTATTGTGCCGGTGCCTGAGAGCCCCGCAAGGCCCCGGACAAGTGATGCTTTGCCGGTTTGGGGCGGTCCGATAATc                                                                                                          <  2:286273/139‑1 (MQ=255)
                               cgAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGTCAGCATGTTTTTTTTGCCGGTGATTGAGAGCCTCGCAAGGTCCCGTATAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCa                                                                                                         <  7:83870/139‑1 (MQ=255)
                               cgAAGACATCCGGCAATCTCTGACCGTTCATCGACGAGCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCa                                                                                                         <  8:259994/139‑1 (MQ=255)
                                gAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACCGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAg                                                                                                        <  8:88897/139‑1 (MQ=255)
                                             aaTCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTcagccag                                                                                           <  7:92113/139‑1 (MQ=255)
                                              aTCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGGCCGATAATCAGCGTATTcagccag                                                                                           <  2:285881/138‑1 (MQ=255)
                                              aTCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGGTGTTTTGCCGGTTTGGGGCGGGCCGATAATCAGCGTATTcagccag                                                                                           >  1:285881/1‑138 (MQ=255)
                                                         ttgaTCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGatat                                                                               <  5:357210/136‑1 (MQ=255)
                                                         ttaaTCGTCGCTCCCTGTCTTTACGGGCATCATGTTTTTTTTGCCTGAGCTTGAGAGCGTCGCTAGGTCGCGGAGAAGGGGTGTTTTGTCGGTTTGTGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGatat                                                                               <  2:336670/136‑1 (MQ=255)
                                                         ttCATGGACGATCCCTGTCTTGACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCTAGGTCCCGGAGACGTGTTGTTTTGCCGGTTGGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGatat                                                                               <  5:290560/139‑1 (MQ=255)
                                                         ttCATCGACGATCCTTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGGTGTTTTGGCGGTTTTGGGCGGGCCGAGAATCAGCGTATTAAGCCAGCTGTTTTGata                                                                                >  7:269805/1‑138 (MQ=255)
                                                         ttCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTTTTTTGCCGCTTTGGGGCGGTCCGAGAATCAGCGTATTCAGCCACCTCTTTTtata                                                                                >  5:349285/1‑138 (MQ=255)
                                                         ttCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGGCGGGTTTGGGGGGGCAGCTAATCCGCGTATTAAGCCCGCTGgtttttttt                                                                               >  5:298874/1‑134 (MQ=255)
                                                         ttCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGatat                                                                               <  4:74909/139‑1 (MQ=255)
                                                         ttCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGatat                                                                               <  2:62747/139‑1 (MQ=255)
                                                         ttCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGatat                                                                               <  6:238165/139‑1 (MQ=255)
                                                         ttCATCGACGATCCCTGTCTATACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGatat                                                                               <  5:6269/139‑1 (MQ=255)
                                                         ttCATCGACGATCCCTGGCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGatat                                                                               <  6:197831/139‑1 (MQ=255)
                                                               gacgaTCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGATATAAGTa                                                                          <  6:156611/138‑1 (MQ=255)
                                                                       cTGTCTTTCCGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCCCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGGCCGATAATCAGCGTATTCAGCCAGCTGTTTTGATATAAGTATGGAagcag                                                                 <  1:182936/139‑1 (MQ=255)
                                                                              tACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTTCCGGTTTGGGGGGGGGCCATAATAACGGTAATCAGCCAGCGGGTTTGATATTAGTATTGAAACAGCGGCAcc                                                          >  4:5588/1‑139 (MQ=255)
                                                                              tACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATaa                                                                                                            >  3:331144/1‑89 (MQ=255)
                                                                              tACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATaa                                                                                                            <  4:331144/89‑1 (MQ=255)
                                                                                 gggCTGCATGTTTTTTTTGCCGGTGCTAGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCTTATTCAGCCAGCTGTTTTGATATAAGTATGGAAGCAGCGGCTCCGCa                                                       <  4:214842/139‑1 (MQ=255)
                                                                                 gggCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTTATATTATTATGGAAGCAGCGGCTCCgaa                                                       <  5:307409/139‑3 (MQ=255)
                                                                                    cagcaTGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTTCCGGTTTTGGGCGGGCCGATAATCAGCGCATTTAGCCAGCTGTGTTTATATAAGTATTTAACCAGCGGcgccgccaac                                                    >  3:29508/1‑135 (MQ=255)
                                                                                    cagcaTGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGGCGGTTTGGGGCGGGCCGATAATCAGCGTATTCAGCCCGCTGTTTTGTTATAAGTATGGAAGCAGCGGCTCCGCAATc                                                    >  2:59887/1‑139 (MQ=255)
                                                                                    cagcaTG‑TTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGc                                                                                                                                 <  3:84841/61‑1 (MQ=255)
                                                                                    cagcaTG‑TTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGc                                                                                                                                 >  4:84841/1‑61 (MQ=255)
                                                                                             ttttttGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGATATAAGTATGGAAGCAGCGGCTCCGCAATCCCCAGCttt                                           >  5:147600/1‑139 (MQ=255)
                                                                                                tttGCCGGTGCTTGAGCGCCTCGCAAGGGCCCGGAGACGGGATGGATTGCCGGTTGGGGGCGGGCCGATAATCCGCGTATGCAGCCAGCTGTTTTGATATAAGTATGGAAGCAGCGGCTCCGCAATCCCCAGCTTTTg                                         <  2:78324/138‑1 (MQ=255)
                                                                                                    ccGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGGTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGATATAAGTATGGAAGCAGCGGCTCCGCAATCCCCAGCTTTTGCCggg                                    <  2:313576/139‑1 (MQ=255)
                                                                                                                  ccTCGCAAGGTCCCGGAGAAGTGGTGTTTTGGCGGGTTGGGGCGGGGCGATAATAAGCGGATTCCGCGCGCGGTTTTGATTTTAGTATGGCAGCAGCGGGTCCGCAATCCCCAGCTTTTGCCGGGCCATTCggttttt                       >  6:202486/1‑132 (MQ=255)
                                                                                                                     cGCAAGGTCCCGGAGAAGTGTTGTTTTGGCGGTTTGGGGCGGTGCGGTAATAAGCGTATTAAGCCAGCGGTTTTTCTATAAGTATGGAAGCAGCGGCTCCCCAAAACCCACCGTTTTGCGGGGAATTAGTATTTTaaa                    >  6:8102/1‑138 (MQ=255)
                                                                                                                           ggcccGGAGAAGTGTTGGTTTCCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCAAGCTGTTTTTATATAAGTATGGAAGCAGCGGCTCCGCAATCCCCAGCTTTTGCCGGGCAATTCGTATATTAAAAGATGCg             <  6:241952/137‑1 (MQ=255)
                                                                                                                                       tgttgtTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGATATAAGTATGGAAGCAGCGGCTCCGCAATCCCCAGCTTTTGCCGGGCAATTCGTATATTAAAAGATGCGATGTCACGGAg  <  3:5588/138‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TCTCTGACCAAATTGATGCTGCGGAATACCGCGAAGACATCCGGCAATCTCTGACCGTTCATCGACGATCCCTGTCTTTACGGGCAGCATGTTTTTTTTGCCGGTGCTTGAGAGCCTCGCAAGGTCCCGGAGAAGTGTTGTTTTGCCGGTTTGGGGCGGTCCGATAATCAGCGTATTCAGCCAGCTGTTTTGATATAAGTATGGAAGCAGCGGCTCCGCAATCCCCAGCTTTTGCCGGGCAATTCGTATATTAAAAGATGCGATGTCACGGAG  >  CP006881/2354892‑2355164

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