Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,862,722 Δ1 bp 8.0% coding (3151/3489 nt) smc ← chromosome partition protein Smc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,862,7220C. 8.0% 151.8 / 7.0 50coding (3151/3489 nt)smcchromosome partition protein Smc
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base . (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

AGCGGCGCATCGACCTCGTCGAGCATGCAGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCATTTATCTGATCAAAGGTATCCTTGAAGCGGTTGCGGGTTTCCTTGTCG  >  NC_002947/4862580‑4862870
                                                                                                                                              |                                                                                                                                                     
aGCGGCGCATCGACCTCGTCGAGCATGCAGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGAt                                                                                                                                                 >  3:38893/1‑149 (MQ=255)
 gcggcgCATCGACCTCGTCGAGCATGCAGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATc                                                                                                                                                >  6:282576/1‑149 (MQ=255)
    gcgcATCGACCTCGTCGAGCATGCAGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTc                                                                                                                                             <  3:152254/149‑1 (MQ=255)
        aTCGACCTCGTCGAGCATGCAGAACGGTGCGGGGTTCAACTTGAAGATAGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCAccc                                                                                                                                         <  3:27499/149‑1 (MQ=255)
         tCGACCTCGTCGAGCATGCAGAACGGTGCGGGGTTCAACGTGAAGATGCCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCAcccc                                                                                                                                        <  5:178291/149‑1 (MQ=255)
           gACCTCGTCGAGCATGCAGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCt                                                                                                                                       <  5:282576/148‑1 (MQ=255)
                    gAGCATGCAGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGt                                                                                                                              <  4:183546/148‑1 (MQ=255)
                           cAGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGt                                                                                                                                              >  5:6774/1‑125 (MQ=255)
                           cAGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGt                                                                                                                                              <  6:6774/125‑1 (MQ=255)
                            aGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCg                                                                                                                     <  2:153846/149‑1 (MQ=255)
                            aGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCg                                                                                                                     >  4:194571/1‑149 (MQ=255)
                                          ttCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTc                                                                                                          >  7:72298/1‑146 (MQ=255)
                                             aaCTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGAt                                                                                                     >  3:73302/1‑148 (MQ=255)
                                             aaCTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGAt                                                                                                     >  6:203515/1‑148 (MQ=255)
                                                ttGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGa                                                                                                                                                  >  2:151047/1‑100 (MQ=255)
                                                ttGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGa                                                                                                                                                  <  1:151047/100‑1 (MQ=255)
                                                ttGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAgc                                                                                                >  7:74680/1‑149 (MQ=255)
                                                  gAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAAGGGAGGGTGCTGTTCTTCTTGCCCGGTGGGCGCCCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAgcgc                                                                                              >  8:227178/1‑149 (MQ=255)
                                                     gATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGc                                                                                           >  3:111494/1‑149 (MQ=255)
                                                     gATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGc                                                                                           >  7:166802/1‑149 (MQ=255)
                                                       tGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCa                                                                                         <  8:72298/149‑1 (MQ=255)
                                                          cAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCAccc                                                                                      >  2:246503/1‑149 (MQ=255)
                                                           aaaCACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCAcccc                                                                                     >  1:308293/1‑149 (MQ=255)
                                                           aaaCACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCAccc                                                                                      >  6:149911/1‑148 (MQ=255)
                                                           aaaCACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCAccc                                                                                      <  5:149911/148‑1 (MQ=255)
                                                                ccAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCTGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGTTCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGaaa                                                                                <  5:81334/149‑1 (MQ=255)
                                                                 cAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGaaaa                                                                               >  5:213902/1‑149 (MQ=255)
                                                                      ccAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACtttt                                                                          <  8:230771/149‑1 (MQ=255)
                                                                      ccAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACtttt                                                                          <  3:67767/149‑1 (MQ=255)
                                                                       cAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGgcgcgc                                                                                                                                                       >  1:212173/1‑72 (MQ=255)
                                                                       cAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGgcgcgc                                                                                                                                                       <  2:212173/72‑1 (MQ=255)
                                                                             ggTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGTGCGCGACATGATCGTCACCCCTGTATAGAGTAAATCTTCGCCCGTAAGTTCCAGATAAGCGCTGCCACGCCCGAAAACTTTTGGGaaaa                                                                   >  3:199014/1‑149 (MQ=255)
                                                                                   tGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCt                                                             <  4:73302/149‑1 (MQ=255)
                                                                                     ccTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACtt                                                                            <  1:261008/132‑1 (MQ=255)
                                                                                     ccTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACtt                                                                            >  2:261008/1‑132 (MQ=255)
                                                                                     ccTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTg                                                            >  8:204649/1‑148 (MQ=255)
                                                                                            gcggccgGAGATCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTATATCTTCGCCCGTCATTTCCAGATAAGGGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCgg                                                    <  1:248918/146‑1 (MQ=255)
                                                                                            gccgccGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGACACCCCCGAAAACATTTGGGAAAAGTGCCTGTAATCCgg                                                    >  2:248918/1‑149 (MQ=255)
                                                                                              cgccgGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCa                                                  <  5:203515/149‑1 (MQ=255)
                                                                                                  ggACAGCAAATGGATGGTTCTGTTCTTCTTGCCCGGTGGGCGCGCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATcc                                                      <  2:231313/140‑1 (MQ=255)
                                                                                                  ggACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCgagatcgg                                              >  1:231313/1‑141 (MQ=255)
                                                                                                  ggACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCATGATCGGCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATcc                                                      <  2:235297/140‑1 (MQ=255)
                                                                                                  ggACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCATGAGCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCgagatcgg                                              >  1:235297/1‑141 (MQ=255)
                                                                                                        cAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTg                                                                         <  7:214156/116‑1 (MQ=255)
                                                                                                        cAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTg                                                                         >  8:214156/1‑116 (MQ=255)
                                                                                                              gatggTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATACAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCATTTATCTGATCaaa                                    <  1:246503/148‑1 (MQ=255)
                                                                                                                        ttcttcTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTGCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCATTTATCTGATCAAAGGTATCCTTg                          <  8:74680/147‑1 (MQ=255)
                                                                                                                        ttcttcTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCATTTATCTGATCAAAGGTATCCTTg                          >  2:219776/1‑147 (MQ=255)
                                                                                                                             cttGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCATTTATCTGATCAAAGGTATCCTTGAAGCGGt                   >  7:25692/1‑149 (MQ=255)
                                                                                                                                gCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCATTTATCTGATCAAAGGTATCCTTGAAGCGGtt                  <  6:213902/147‑1 (MQ=255)
                                                                                                                                     gTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCATTTATCTGATCAAAGGTATCCTTGAAGCGGTTGCGGGtt           <  4:38893/149‑1 (MQ=255)
                                                                                                                                      tGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCATTTATCTGATCAAAGGTATCCTTGAAGCGGTTGCGGGttt          <  8:25692/149‑1 (MQ=255)
                                                                                                                                              ccATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCATTTATCTGATCAAAGGTATCCTTGAAGCGGTTGCGGGTTTCCTTGTCg  <  8:166802/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
AGCGGCGCATCGACCTCGTCGAGCATGCAGAACGGTGCGGGGTTCAACTTGAAGATGGCAAACACCAGCGCCAAGGCGGTCAGTGCCTTCTCGCCGCCGGACAGCAAATGGATGGTGCTGTTCTTCTTGCCCGGTGGGCGCGCCATGATCGTCACCCCTGTATCGAGTAAATCTTCGCCCGTCAGTTCCAGATAAGCGCTGCCACCCCCGAAAACTTTTGGGAAAAGTGCCTGTAATCCGGCATTTATCTGATCAAAGGTATCCTTGAAGCGGTTGCGGGTTTCCTTGTCG  >  NC_002947/4862580‑4862870

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