Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,397,239 Δ10,478 bp IS5‑mediated ynaJ[dgcM] ynaJ, uspE, fnr, ogt, abgT, abgB, abgA, abgR, mcaS, smrA, [dgcM]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1396205–1397243 1407716 10474–11512 23 [22] [0] 79 [insH4]–[dgcM] [insH4],ynaJ,uspE,fnr,ogt,abgT,abgB,abgA,abgR,mcaS,smrA,[dgcM]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1397238NA (NA)79 (0.900) 51/598 0.5 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 1407717 = 0 (0.000)coding (79/1233 nt) dgcM diguanylate cyclase DgcM

AAACAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCTGAGCGATGGCG  >  NC_000913/1395905‑1396337
                                                                                                                                                                                                                                                                                                           |                                                                                                                                     
aaaCAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGAtt                                                                                                                                       <  1:61734/300‑1 (MQ=255)
 aaCAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTc                                                                                                                                      <  1:410129/300‑1 (MQ=255)
     gTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGCCAATTCATTCAGATCATCAATATTCAGGGAAGAAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCAATAGAACAGGGTTCATCATGAGTCATCAACTTCCCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGACTGGAGCAGATTCTGcc                                                                                                                                  >  2:278274/1‑300 (MQ=255)
     gTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGcc                                                                                                                                  <  2:77882/300‑1 (MQ=255)
               cATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATNTGAGATCATGTTTGTCATCTGGAGCCNTAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAAc                                                                                                                        <  2:462386/300‑1 (MQ=255)
                   gTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACAtgg                                                                                                                    >  1:378011/1‑300 (MQ=255)
                                cctccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGa                                                                                                       >  1:640852/1‑300 (MQ=255)
                                          aTTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATTATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCCGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTAATCGAGCCGTTTTAc                                                                                             >  2:215432/1‑300 (MQ=255)
                                            tGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTAccc                                                                                           <  1:370851/300‑1 (MQ=255)
                                                            aCCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCACGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGc                                                                           <  1:36803/300‑1 (MQ=255)
                                                                           aTTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCg                                                            >  1:365810/1‑300 (MQ=255)
                                                                                   gcgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAAc                                                    <  1:426876/300‑1 (MQ=255)
                                                                                        aaTTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGc                                               <  2:373093/300‑1 (MQ=255)
                                                                                                 aGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATTAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCATAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGTTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTc                                      <  2:266795/300‑1 (MQ=255)
                                                                                                     catcaATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACtg                                  >  1:680187/1‑300 (MQ=255)
                                                                                                             gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATgcagc                          <  1:546052/300‑1 (MQ=39)
                                                                                                                 gggAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTg                      >  1:89312/1‑300 (MQ=38)
                                                                                                                     aggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTNCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTAc                  <  1:480459/300‑1 (MQ=25)
                                                                                                                         agTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAAcc              <  1:445970/300‑1 (MQ=34)
                                                                                                                          gTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCt             <  1:374237/300‑1 (MQ=34)
                                                                                                                                aCATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGACTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCTGAGCGa       >  1:513941/1‑300 (MQ=17)
                                                                                                                                 cATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCTGAGCGAt      >  1:549261/1‑300 (MQ=31)
                                                                                                                                     aGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCTGAGCGATGgcg  >  1:344070/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                     
AAACAGTTCAGCGATCATTGTTTTATCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCTGAGCGATGGCG  >  NC_000913/1395905‑1396337

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