Predicted mutation
evidence seq id position mutation freq annotation gene description
JC CP009273 2,656,840:1 (GCGCGCA)2→3 81.2% coding (40/804 nt) suhB → inositol monophosphatase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 = 26568263 (0.140)15 (0.670)
+GCGCGCA
16/228 NT 81.2% coding (26/804 nt) suhB inositol monophosphatase
?CP009273 2656827 = 4 (0.170)coding (27/804 nt) suhB inositol monophosphatase

TTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAA  >  CP009273/2656687‑2656970
                                                                                                                                            |                                                                                                                                                 
tttcttCATCTAAAGATTATTCAGGCATGTTATCGTAAAACGAAGACAGATGCAGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagc                                                                                                                                          >  2:239432‑M1/1‑140 (MQ=255)
              gATTATTCACGCATCTTATCATAAAACGAAGACAGATGCGGATCTCGCTGCTATACTCTGCGCCGTTTTCCTGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGATCATCGCCGTGCGCGCAGCGCGCAAgg                                                                                                                                  >  1:102379/1‑143 (MQ=255)
              gATTATTCACGCATCTTATCAGAAAACGAAGACAGATGCCGATCTCGCTGCGATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGta                                                                                                                            <  2:120146/149‑1 (MQ=255)
                  attCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGCTGATGACACACCGGGTgcgcgcagcgcgcagcgcgcaaggcgg                                                                                                                        <  2:212039‑M1/149‑28 (MQ=255)
                                  aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaagccgggaaatttcaatgccaa                                                                                                        >  2:110710‑M1/1‑106 (MQ=255)
                                  aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGGTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggaactttacttgacaa                                                                                                        >  2:301344‑M1/1‑106 (MQ=255)
                                   tAAAACGAAGACAGATGCCGATCTCGCTGCTATACTTTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaatt                                                                                                             <  1:309513‑M1/143‑39 (MQ=255)
                                    aaaaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattg                                                                                                            >  1:46742‑M1/1‑104 (MQ=255)
                                         gAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgagcagcgccccaggcgggtcatttcattgcccaa                                                                                                       >  1:180550‑M1/1‑99 (MQ=255)
                                                  tGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaa                                                                                              <  1:123134‑M1/143‑54 (MQ=255)
                                                   gCCGATCTCGCTGCTATACTCTGCGCCGTTTACCCGTTCTTTAACATCCAGTGAGCGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcggagcgcgcagcgcgcaaggcgggtaagttaattgccaaaaaacttgtaa                                                                                             >  1:55445‑M1/1‑89 (MQ=255)
                                                        tctcGCTGCTATACCCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgct                                                                                  <  2:29685‑M1/149‑66 (MQ=255)
                                                                    aCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgctgtagaa                                                                            <  1:84221‑M1/143‑72 (MQ=255)
                                                                    aCTCTGCGCCGTTTTCCCGTTCTTTAACAGCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgctgtagaa                                                                            <  1:125935‑M1/143‑72 (MQ=255)
                                                                                     cGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGCAACCCCGGACGCTGTAGAAGCGAGCCAGCAAGGCAGTAACGATTTCGTg                                                    >  2:285115/1‑149 (MQ=255)
                                                                                                        agagggaCCGATGCATCCGAGGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaatggccaaaaactatgaaaccccggacgctgtagaagcgagccagaaaggcagtaacgatttcgtgaccaacgtagat                                  <  2:60246‑M1/149‑114 (MQ=255)
                                                                                                        agagagaCCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgctgtagaagcgagccagaaaggcagtaacgatttcgtgaccaacgtagat                                  >  2:32746‑M1/1‑36 (MQ=255)
                                                                                                                      aTCCGAGGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgctgtagaagcgagccagaaaggcagtaacgatttcgtgaccaacgtagataaagctgccgaagc                    <  2:53174‑M1/149‑128 (MQ=255)
                                                                                                                         cGATGCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgctgtagaagcgagccagaaaggcagtaacgatttcgtgaccaacgtagataaagctgccgaagcggt                 >  2:202348‑M1/1‑19 (MQ=255)
                                                                                                                             gCTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgctgtagaagcgagccagaaaggcagtaacgatttcgtgaccaacgtagataaagctgccgaagcggtgatt             >  2:197622‑M1/1‑15 (MQ=255)
                                                                                                                              cTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgctgtagaagcgagccagaaaggcagtaacgatttcgtgacccacgtcgatacagctgccgaagcggtgatta            >  2:30291‑M1/1‑14 (MQ=255)
                                                                                                                              cTGAACATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgctgtagaagcgagccagaaaggcagtaacgaattcgtgaccaacgtagataaagctgccgagtggg                  >  1:241873‑M1/1‑14 (MQ=255)
                                                                                                                               tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGCAGCGGTGATTATc                >  1:50099/1‑143 (MQ=255)
                                                                                                                                 aaCATCGCCGTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgctgtagaagcgagtcagaaaggcagtaacgatttcgtggccaacgtagataaagctgccgaagcggtgattatcg         >  2:196584‑M1/1‑11 (MQ=255)
                                                                                                                                        gCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTaa  >  2:375613/1‑149 (MQ=255)
                                                                                                                                          cGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACACCGGCCGACGTAGAAGCGAGCCAGAAAGGAAGCAACGGGCGGCCTGCCGACGTAAACAAATCGGCCGAAGCGGTGATTATCGACACGAtcc      >  1:374330/1‑141 (MQ=255)
                                                                                                                                           gTgcgcgcagcgcgcagcgcgcaaggcgggtaatttaattgccaaaaactatgaaaccccggacgctgtagaagcgagccagaaaggcagtaacgatttcgtgaccaacgtagataaagctgccgaagcggtgattatcgaca      >  1:268160‑M1/1‑2 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
TTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAA  >  CP009273/2656687‑2656970

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.