Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,176,353 (GCA)4→3 coding (479‑481/828 nt) cpdA ← 3',5' cAMP phosphodiesterase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3176344 3176346 3 29 [0] [0] 31 cpdA 3',5' cAMP phosphodiesterase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3176347 =0 (0.000)29 (0.860) 24/242 0.6 100% coding (487/828 nt) cpdA 3',5' cAMP phosphodiesterase
?NC_000913 = 3176343 0 (0.000)coding (491/828 nt) cpdA 3',5' cAMP phosphodiesterase

AGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGC  >  NC_000913/3176205‑3176490
                                                                                                                                          |                                                                                                                                                   
agaTCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaa                                                                                                                                           <  4:125289‑M2/149‑11 (MQ=255)
                 aaTATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggc                                                                                                                           >  1:106916‑M2/1‑122 (MQ=255)
                                      tttGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtt                                                                                                              >  8:21095‑M2/1‑101 (MQ=255)
                                      tttGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtt                                                                                                              <  7:21095‑M2/140‑40 (MQ=255)
                                                        cGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaa                                                                                     >  3:125289‑M2/1‑83 (MQ=255)
                                                               aCGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcag                                                                            <  2:87222‑M2/149‑74 (MQ=255)
                                                                    aTCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttcctgcca                                                                                                 >  6:122110‑M2/1‑71 (MQ=255)
                                                                    aTCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagcca                                                                                                 <  5:122110‑M2/123‑53 (MQ=255)
                                                                              cccGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggc                                                                                                                           <  6:84485‑M2/87‑27 (MQ=255)
                                                                              cccGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggc                                                                                                                           >  5:84485‑M2/1‑61 (MQ=255)
                                                                                 gcgTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgc                                                           <  5:35619‑M2/148‑91 (MQ=255)
                                                                                         cAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacct                                                  >  1:87222‑M2/1‑50 (MQ=255)
                                                                                          aGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaac                                                      >  8:115621‑M2/1‑49 (MQ=255)
                                                                                          aGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaac                                                      <  7:115621‑M2/144‑96 (MQ=255)
                                                                                           gACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctgg                                                >  6:35619‑M2/1‑48 (MQ=255)
                                                                                                    gATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctg                                                 <  5:163604‑M2/139‑101 (MQ=255)
                                                                                                    gATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctg                                                 >  6:163604‑M2/1‑39 (MQ=255)
                                                                                                           ccAACTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctggctatccagcaacagga                                <  2:121121‑M2/149‑118 (MQ=255)
                                                                                                             aaCTACAGCCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctggctatccagcaacaggatt                              <  2:91005‑M2/149‑120 (MQ=255)
                                                                                                             aaCTACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctggctatccagcaacaggatt                              <  8:74818‑M2/149‑120 (MQ=255)
                                                                                                                tACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctgg                                                >  5:42920‑M2/1‑27 (MQ=255)
                                                                                                                tACAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctgg                                                <  6:42920‑M2/128‑102 (MQ=255)
                                                                                                                  cAACCCGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctggctatccagcaacaggatttgc                           >  1:121121‑M2/1‑25 (MQ=255)
                                                                                                                     cccGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctggctatccagcaacaggatttgcca                         >  2:25016‑M2/1‑22 (MQ=255)
                                                                                                                     cccGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctggctatccagcaacaggatttgcca                         <  1:25016‑M2/146‑125 (MQ=255)
                                                                                                                     cccGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccaatt                                                                                                              >  4:66069‑M2/1‑22 (MQ=255)
                                                                                                                     cccGCAGGTAGCGGATGATGATgcagcagcaacgtatggcgttctggcgcatcggccaatt                                                                                                              <  3:66069‑M2/61‑40 (MQ=255)
                                                                                                                                   atgatgatgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctggctatccagcaacaggatttgccattgctcaccaataaaca        >  1:91005‑M2/1‑8 (MQ=255)
                                                                                                                                         atgcagcagcaacgtatggcgttctggcgcatcggccagtttacgttccagccactcaagctgaaactcgctcagctcaccgtgcggcacgccaaacacctgcctatccagcaacaggatttgccattgctcaccaataaacacgcgct  <  2:191683‑M2/149‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
AGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGATGATGCAGCAGCAGCAACGTATGGCGTTCTGGCGCATCGGCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGC  >  NC_000913/3176205‑3176490

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