Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC CP009273 2,809,385 Δ148 bp 100% coding (267‑414/1557 nt) gshA ← glutamate‑cysteine ligase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ CP009273 2809385 2809532 148 25 [0] [0] 25 gshA glutamate‑cysteine ligase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 = 28093840 (0.000)21 (0.840) 21/206 NT 100% coding (415/1557 nt) gshA glutamate‑cysteine ligase
?CP009273 2809533 = 0 (0.000)coding (266/1557 nt) gshA glutamate‑cysteine ligase

GCAGGTAATACATACCGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATACAGCGTTTTAAAGCGTCCGGTGTTAGAAGTGCCGTACTGTGCCAGTTCGATGTCCTGACCTTCTGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATT  >  CP009273/2809114‑2809599
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                       
tcaggtactcanaaccgcacgcgggttgCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAAGAATGACCCAACCGAAACGATAGTAAGTGCGGAGAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaac                                                                                                                                                                                                                     <  2:70893‑M1/262‑5 (MQ=255)
                 cACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGCTGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTTCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATAAGCGCCCGAGATATCACCGCACTTCGTTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgaaaa                                                                                                                                                                                                    >  1:383894‑M1/1‑254 (MQ=255)
                                      aGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaata                                                                                                                                                                               >  1:27255‑M1/1‑233 (MQ=255)
                                                 tgttccTTGAGGGTAAGAAGAACAATACGCCGGAGAAGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatcaccatccac                                                                                                                                                                    <  2:27255‑M1/273‑54 (MQ=255)
                                                     ccTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacaatgcagatcgcgcataaaggtcagcataggttcaatataantatccaatggg                                                                                                                                                                >  2:61509‑M1/1‑218 (MQ=255)
                                                         gCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatcaccatcca                                                                                                                                                                     <  1:61509‑M1/266‑53 (MQ=255)
                                                                           tCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTTGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatcaccatccactggtgtaatgaattccagcaatgctt                                                                                                                                          >  1:282361‑M1/1‑196 (MQ=255)
                                                                                                 aTAAGGAATGACCCAACCGAAACGATAGGAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcngatcgcgcataaaggtcagcatatgttcaatatcaccatccactggtgtaatgaattccagcaatgcttccgcaaaatcggtagtaatcaa                                                                                                                    >  1:521084‑M1/1‑174 (MQ=255)
                                                                                                                                                                    tctTTGGCATGAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatcaccatccactggtgtaatgaattccagcaatgcttc                                                                                                                                         >  1:650925‑M1/1‑107 (MQ=255)
                                                                                                                                                                        tGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatcaccatccactggtgtaatgaattccagcaatgcttccgca                                                                                                                                     <  2:650914‑M1/187‑85 (MQ=255)
                                                                                                                                                                                         atatCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtaagcatatgttcaat                                                                                                                                                                                >  2:11253‑M1/1‑86 (MQ=255)
                                                                                                                                                                                           atCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatat                                                                                                                                                                              >  2:443514‑M1/1‑84 (MQ=255)
                                                                                                                                                                                             cACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtaagcatatgttcaatatca                                                                                                                                                                            <  1:11253‑M1/127‑46 (MQ=255)
                                                                                                                                                                                               ccGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatcacc                                                                                                                                                                          <  1:443520‑M1/127‑48 (MQ=255)
                                                                                                                                                                                                  cACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatccccatccactggtgtaatgaattccagcaatgcttccgcaaaatcggtagtaatccatttgtgcgtcagtgcggaacctaatgcttcaggatgacctgttgttgccagtgtgccatcagcattaacacgcaaagtttcgcgctccagcc                        >  1:65956‑M1/1‑77 (MQ=255)
                                                                                                                                                                                                   cctgttattnacagaaggcctttGGCAAAGAGAAATTGTATTGAACGCCGGAAATGGTTGGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatcaccatacactggtgtaatgaattccagaaatgcttccgcaaaatcggtagtaatacatttgtgcgtcagtgcggaacctaatgcttcatgatgacctgttgttgccagtgtgccatcagcattaacacgcaaagtttcgcgctccagcccacg                    <  2:65956‑M1/262‑198 (MQ=255)
                                                                                                                                                                                                   aCTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggt                                                                                                                                                                                               >  2:156144‑M1/1‑76 (MQ=255)
                                                                                                                                                                                                       cGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagc                                                                                                                                                                                           <  1:156147‑M1/102‑31 (MQ=255)
                                                                                                                                                                                                                   tGCCATTGTAAAAAAGAAATTGTAGGGCACCCCCGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcaggtcgcgcataaaggtcagcatatgttcaatatcaccatccactggtgtaatgaattccagcaatgcttccgcaaaatcggtagtaatccatttgtgcgtcagtgcggaacctaatgcttcaggatgacctgttgttgccagtgtgccatcagcattaacacgcaaagtttcgcgctccagcccacgctgtatcccctttaacgc  <  2:282356‑M1/275‑216 (MQ=255)
                                                                                                                                                                                                                        ttGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttca                                                                                                                                                                                  >  1:587403‑M1/1‑55 (MQ=255)
                                                                                                                                                                                                                            cAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatat                                                                                                                                                                              <  2:587392‑M1/94‑44 (MQ=255)
                                                                                                                                                                                                                                                             ttGCATCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatcaccatccactggtgtaatgaattccagcaatgcttccgcaaaatcggtagtaatccatttgtgcgtcagtgcggaa                                                                                                 >  1:583541‑M1/1‑18 (MQ=255)
                                                                                                                                                                                                                                                                 aTCAGCGCGCCGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatcaccatccactggtgtaatgaattccagcaatgcttccgcaaaatcggtagtaatccatttgtgcgtcagtgcggaaccta                                                                                             <  2:583530‑M1/138‑125 (MQ=255)
                                                                                                                                                                                                                                                                          ccGTAtaacgatgcagatcgcgcataaaggtcagcatatgttcaatatcaccatccactggtgtaatgaattccagcaatgcttccgcaaaatcggtagtaatccatttgtgcgtcagtgcggaacctaatgctt                                                                                       >  2:455275‑M1/1‑5 (MQ=255)
                                                                                                                                                                                                                                                                              ataacgatgcagatcgcgcataaaggtcagcatatgttcaatatcaccatccactggtgtaatgaattccagcaatgcttccgcaaaatcggtagtaatccatttgtgcgtcagtgcggaacctaatgcttcagg                                                                                   <  1:455281‑M1/135‑135 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                       
GCAGGTAATACATACCGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATACAGCGTTTTAAAGCGTCCGGTGTTAGAAGTGCCGTACTGTGCCAGTTCGATGTCCTGACCTTCTGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATT  >  CP009273/2809114‑2809599

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