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

GCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATACAGCGTTTTAAAGCGTCCGGTGTTAGAAGTGCCGTACTGTGCCAGTTCGATGTCCTGACCTTCTGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGC  >  CP009273/2809262‑2809747
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                      
tcaggtactcanaaccgcacgcgggttgctcaaacggcagcgacgttggttttccttgcaggaaagaagaacaaatcgccggagatgcaccaaacagataaagaatgacccaaccgaaacgatagtaagtgcggagaacgcggaaatagcccgcagaaattttctctttggcatcagcgcccgagatatcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAAc                                                                                                                                                                                                                     <  2:70893‑M2/4‑1 (MQ=255)
                 cactcggttttctcaaacggcagcgacgttggttttccttgcaggaaagaagaacaaatcgccggagctgcaccaaacagataaggaatgacccaaccgaaacgatagtaatttcggataacgcggaaatagcccgcagaaattttctctttggcataagcgcccgagatatcaccgcacttcgtttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCagatcgcgaaaa                                                                                                                                                                                                    >  1:383894‑M2/255‑271 (MQ=255)
                                      agcgacgttggttttccttgcaggaaagaagaacaaatcgccggagatgcaccaaacagataaggaatgacccaaccgaaacgatagtaattgcggataacgcggaaatagcccgcagaaattttctctttggcatcagcgcccgagatatcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAata                                                                                                                                                                               >  1:27255‑M2/234‑275 (MQ=255)
                                                 tgttccttgagggtaagaagaacaatacgccggagaagcaccaaacagataaggaatgacccaaccgaaacgatagtaattgcggataacgcggaaatagcccgcagaaattttctctttggcatcagcgcccgagatatcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCAc                                                                                                                                                                    <  2:27255‑M2/53‑1 (MQ=255)
                                                     ccttgcaggaaagaagaacaaatcgccggagatgcaccaaacagataaggaatgacccaaccgaaacgatagtaattgcggataacgcggaaatagcccgcagaaattttctctttggcatcagcgcccgagatatcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACAATGCAGATCGCGCATAAAGGTCAGCATAGGTTCAATATAANTATCCAAtggg                                                                                                                                                                >  2:61509‑M2/219‑274 (MQ=255)
                                                         gcaggaaagaagaacaaatcgccggagatgcaccaaacagataaggaatgacccaaccgaaacgatagtaattgcggataacgcggaaatagcccgcagaaattttctctttggcatcagcgcccgagatatcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCa                                                                                                                                                                     <  1:61509‑M2/52‑1 (MQ=255)
                                                                           tcgccggagatgcaccaaacagataaggaatgacccaaccgaaacgatagtaattgcggataacgcggaaatagcccgcagaaattttctcttttgcatcagcgcccgagatatcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCtt                                                                                                                                          >  1:282361‑M2/197‑275 (MQ=255)
                                                                                                 ataaggaatgacccaaccgaaacgataggaattgcggataacgcggaaatagcccgcagaaattttctctttggcatcagcgcccgagatatcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCNGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAatcaa                                                                                                                    >  1:521084‑M2/175‑273 (MQ=255)
                                                                                                                                                                    tctttggcatgagcgcccgagatatcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTc                                                                                                                                         >  1:650925‑M2/108‑187 (MQ=255)
                                                                                                                                                                        tggcatcagcgcccgagatatcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCa                                                                                                                                     <  2:650914‑M2/84‑1 (MQ=255)
                                                                                                                                                                                         atatcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTAAGCATATGTTCAat                                                                                                                                                                                >  2:11253‑M2/87‑127 (MQ=255)
                                                                                                                                                                                           atcaccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAatat                                                                                                                                                                              >  2:443514‑M2/85‑127 (MQ=255)
                                                                                                                                                                                             caccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTAAGCATATGTTCAATATCa                                                                                                                                                                            <  1:11253‑M2/45‑1 (MQ=255)
                                                                                                                                                                                               ccgcacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCAcc                                                                                                                                                                          <  1:443520‑M2/47‑1 (MQ=255)
                                                                                                                                                                                                  cacttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCCCCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGcc                        >  1:65956‑M2/78‑270 (MQ=255)
                                                                                                                                                                                                   cctgttattnacagaaggcctttggcaaagagaaattgtattgaacgccggaaatggttggcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATACACTGGTGTAATGAATTCCAGAAATGCTTCCGCAAAATCGGTAGTAATACATTTGTGCGTCAGTGCGGAACCTAATGCTTCATGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACg                    <  2:65956‑M2/197‑1 (MQ=255)
                                                                                                                                                                                                   acttcgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGt                                                                                                                                                                                               >  2:156144‑M2/77‑102 (MQ=255)
                                                                                                                                                                                                       cgcttgccagaatgccattggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGc                                                                                                                                                                                           <  1:156147‑M2/30‑1 (MQ=255)
                                                                                                                                                                                                                   tgccattgtaaaaaagaaattgtagggcacccccgaaatggtttgcatcagcgcgccgtaTAACGATGCAGGTCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGc  <  2:282356‑M2/215‑1 (MQ=255)
                                                                                                                                                                                                                        ttggcaaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCa                                                                                                                                                                                  >  1:587403‑M2/56‑94 (MQ=255)
                                                                                                                                                                                                                            caaagagaaattgtagtgcacgccggaaatggtttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAatat                                                                                                                                                                              <  2:587392‑M2/43‑1 (MQ=255)
                                                                                                                                                                                                                                                             ttgcatcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGaa                                                                                                 >  1:583541‑M2/19‑138 (MQ=255)
                                                                                                                                                                                                                                                                 atcagcgcgccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTa                                                                                             <  2:583530‑M2/124‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ccgtaTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCtt                                                                                       >  2:455275‑M2/6‑135 (MQ=255)
                                                                                                                                                                                                                                                                              aTAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAgg                                                                                   <  1:455281‑M2/134‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                      
GCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATACAGCGTTTTAAAGCGTCCGGTGTTAGAAGTGCCGTACTGTGCCAGTTCGATGTCCTGACCTTCTGCGATGTAGCATGGCATACTTAACGGCCACATCCGCTCATCGCCCATATTGCGCGCCGTATAACGATGCAGATCGCGCATAAAGGTCAGCATATGTTCAATATCACCATCCACTGGTGTAATGAATTCCAGCAATGCTTCCGCAAAATCGGTAGTAATCCATTTGTGCGTCAGTGCGGAACCTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGC  >  CP009273/2809262‑2809747

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: