Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 970,407 Δ1 bp 10.5% coding (633/1323 nt) mukF → chromosome condensin MukBEF, kleisin‑like subunit, binds calcium

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273970,4070T.10.5% 83.7 / 7.3 38coding (633/1323 nt)mukFchromosome condensin MukBEF, kleisin‑like subunit, binds calcium
Reads supporting (aligned to +/- strand):  ref base T (17/17);  new base . (2/2);  total (19/19)
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.76e-01

TGGCACCGTAATGTCTATGCGCCACTGAAATATTCGGTAGCAGAAATTTTCGACAGTATCGACCTGACGCAACGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCA  >  CP009273/970228‑970670
                                                                                                                                                                                   |                                                                                                                                                                                                                                                                       
tGGCACCGTAATGTCTATGCGCCACTGAAATATTCGGTAGCAGAAATTTTCGACAGTATCGACCTGACGCAACGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTg                                                                                                                                                                                                                   >  1:113516/1‑234 (MQ=255)
    aCCGTAATGTCTATGCGCCACTGAAATATTCGGTAGCAGAAATTTTCGACAGTATCGACCTGACGCAACGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAgg                                                                                                                                                                                                               <  2:113515/234‑1 (MQ=255)
                                                                  aCGCAACGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGTTATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGNGCATTAAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCg                                                                                                                     >  1:499683/1‑262 (MQ=255)
                                                                  aCGCAACGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTGGCTGTCGGAAAATTCCGGAACGCTGCGCGAATTGCAGGGTAAGCTGGGAGCGGCAAGCGACAAATTGCAGGCTAATCCGGTGCGCATTCAGGGTCCGACGAGGAACCATGAGGATCCGACTTCAGTAGACAGGCGGGGGTGCGATCTGCAGAGCAAACTCgtcattagtttaa                                                                                                         >  2:238695/1‑262 (MQ=255)
                                                                      anCGGCTGATGGAAGACCAGCAGCAGCAGGTGAAGGACGATATCGCCCATTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCg                                                                                                                 <  2:499674/262‑1 (MQ=255)
                                                                                     aaCAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCa                                                                                                                                                                                            >  1:411354/1‑172 (MQ=255)
                                                                                         gcagcagcagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGAt                                                                                                                                                                                        <  2:411349/172‑1 (MQ=255)
                                                                                          cagcagcagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGTCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTg                                                                                                                                  <  1:359549/225‑1 (MQ=255)
                                                                                            gcagcagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGAAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTgcg                                                                                                                                                                                                  >  2:311128/1‑159 (MQ=255)
                                                                                                cagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAg                                                                                                                                                                                           >  1:36875/1‑162 (MQ=255)
                                                                                                cagGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGAAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCAtt                                                                                                                                                                                              <  1:311133/159‑1 (MQ=255)
                                                                                                   gTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACACATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTAACAGTTGGGGCCAGCAATACATCGACtatcacatat                                                                        >  1:73288/1‑265 (MQ=255)
                                                                                                    tGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATg                                                                                                                                                                                       <  2:36875/162‑1 (MQ=255)
                                                                                                    ggAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATGTACAGCTGTGAATTGTTGCTTTCGGACACTTCCgg                                                                                                                                                                                                                                                                    >  2:200285/2‑85 (MQ=255)
                                                                                                        ggACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACg                                                                                                                                                                                                                                                                <  1:200289/85‑1 (MQ=255)
                                                                                                              tatCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGTAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGCCAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGtt                                                                                                                                         >  1:64853/1‑198 (MQ=255)
                                                                                                                  gCCCAGTTGCTGAACAAAGACTGGCGGGTGGCGATTGCCAGCTGTGAAGTGTTGCTTTCGGAAACTTCCGTAACGCTGCGTGAAGTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGAt                                                                                                                                     <  2:64853/198‑1 (MQ=255)
                                                                                                                       gTTTCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACaa                                                                                                                                                                                                                       >  1:553997/1‑111 (MQ=255)
                                                                                                                           cTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTg                                                                                                                                                                                                                   <  2:553986/111‑1 (MQ=255)
                                                                                                                             gAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCgg                                                                                                                                                                                                                                                                    >  2:72975/1‑60 (MQ=255)
                                                                                                                                 aaaGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACg                                                                                                                                                                                                                                                                <  1:72975/60‑1 (MQ=255)
                                                                                                                                  aaGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCg                                                                                                                                                                                                                                                           >  2:143854/1‑64 (MQ=255)
                                                                                                                                   aGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCg                                                                                                                                                                                                                                  >  2:424720/1‑88 (MQ=255)
                                                                                                                                      cTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGaa                                                                                                                                                                                                                                                       <  1:143856/64‑1 (MQ=255)
                                                                                                                                       tGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAg                                                                                                                                                                                                                              <  1:424726/88‑1 (MQ=255)
                                                                                                                                        ggcgggtgGCGATTTCCAGCTGTGAATTATTGCTTTCGGAAACTGCCGGAACg                                                                                                                                                                                                                                                                >  1:506173/1‑53 (MQ=38)
                                                                                                                                        ggcgggcgGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGc                                                                                                                                                                                                                                                   >  2:255767/1‑66 (MQ=255)
                                                                                                                                        ggcgggcgGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTg                                                                                                                                                                                                                   >  2:216156/1‑98 (MQ=255)
                                                                                                                                        ggcgggcgGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGAt                                                                                                                                                                   >  2:391245/1‑146 (MQ=255)
                                                                                                                                            ggcggcGATTTCCAGCTGTGATTTGTTGCTTTCGGAAACTTCCGGACCGCTGc                                                                                                                                                                                                                                                            <  2:506163/53‑1 (MQ=255)
                                                                                                                                            ggcggcGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGa                                                                                                                                                                                                                                               <  1:255771/66‑1 (MQ=255)
                                                                                                                                            ggcggcGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAgg                                                                                                                                                                                                               <  1:216160/98‑1 (MQ=255)
                                                                                                                                            ggcggcGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGc                                                                                                                                                               <  1:391250/146‑1 (MQ=255)
                                                                                                                                                      tCCAGCTGTGAATTGTTGCTTTCGGAAACTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGc                                                                                                                                                                                                              >  1:490737/1‑88 (MQ=255)
                                                                                                                                                      tCCAGCTGTGAATTGTTGCTTTCGGAAACTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGc                                                                                                                                                                                                              >  1:265431/1‑88 (MQ=255)
                                                                                                                                                          gCTGTGAATTGTTGCTTTCGGAAACTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAAt                                                                                                                                                                                                          <  2:265427/88‑1 (MQ=255)
                                                                                                                                                          gCTGTGAATTGTTGCTTTCGGAAACTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAAt                                                                                                                                                                                                          <  2:490729/88‑1 (MQ=255)
                                                                                                                                                                         ttttgggtaCTTCCGGAACGCTGCGTGAATTTCAGGGTACGCTGGAAGCGTCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCGACGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCa  <  1:238699/266‑1 (MQ=255)
                                                                                                                                                                                   |                                                                                                                                                                                                                                                                       
TGGCACCGTAATGTCTATGCGCCACTGAAATATTCGGTAGCAGAAATTTTCGACAGTATCGACCTGACGCAACGTCTGATGGACGAACAGCAGCAGCAGGTGAAGGACGATATCGCCCAGTTGCTGAACAAAGACTGGCGGGCGGCGATTTCCAGCTGTGAATTGTTGCTTTCGGAAACTTCCGGAACGCTGCGTGAATTGCAGGATACGCTGGAAGCGGCAGGCGACAAATTGCAGGCTAATCTGTTGCGCATTCAGGATGCGACGATGACCCATGACGATCTGCATTTCGTCGATCGTCTGGTGTTCGATCTGCAGAGCAAACTCGATCGTATTATCAGTTGGGGCCAGCAATCCATCGACTTGTGGATTGGCTACGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCA  >  CP009273/970228‑970670

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: