Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,888,432 Δ18,167 bp tnaCbglG 17 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3888432 3906598 18167 71 [0] [0] 71 tnaC–bglG tnaC,tnaA,tnaB,mdtL,yidZ,yieE,yieF,adeP,yieH,cbrB,cbrC,yieK,yieL,bglH,bglB,bglF,bglG

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 38884310 (0.000)66 (0.930) 46/96 0.0 100% intergenic (+239/‑4) mnmE/tnaC 5‑carboxymethylaminomethyluridine‑tRNA synthase GTPase subunit/tnaAB operon leader peptide
?NC_000913 3906599 = 0 (0.000)intergenic (‑32/+254) bglG/phoU transcriptional antiterminator BglG/negative regulator of the pho regulon

AAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTG  >  NC_000913/3888382‑3888480
                                                 |                                                 
aaaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGC                                                   >  1:1785563‑M1/1‑50 (MQ=255)
aaaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGC                                                   >  1:5533457‑M1/1‑50 (MQ=255)
aaaaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGC                                                   >  1:2228156‑M1/1‑50 (MQ=255)
  aaaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCc                                                  >  1:4744198‑M1/1‑48 (MQ=255)
   aaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcga                                                <  1:4769902‑M1/50‑4 (MQ=255)
   aaCTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcga                                                >  1:3620034‑M1/1‑47 (MQ=255)
     ctctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgact                                              >  1:4726828‑M1/1‑45 (MQ=255)
     ctctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgact                                              <  1:470935‑M1/50‑6 (MQ=255)
     ctctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgact                                              >  1:3753373‑M1/1‑45 (MQ=255)
     ctctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgact                                              >  1:187380‑M1/1‑45 (MQ=255)
       ctGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgactt                                             >  1:5571388‑M1/1‑43 (MQ=255)
        tGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgactt                                             <  1:5288793‑M1/48‑7 (MQ=255)
        tGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgactt                                             <  1:908313‑M1/48‑7 (MQ=255)
         gAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcac                                          >  1:1044574‑M1/1‑41 (MQ=255)
         gAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcac                                          >  1:5398480‑M1/1‑41 (MQ=255)
         gAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcac                                          >  1:3009331‑M1/1‑41 (MQ=255)
          aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcac                                          <  1:1181382‑M1/49‑10 (MQ=255)
          aGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttc                                            <  1:3855846‑M1/47‑8 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagta                                     >  1:4350594‑M1/1‑36 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagta                                     <  1:3575902‑M1/50‑15 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagt                                      <  1:3883764‑M1/49‑14 (MQ=255)
              taataaTGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagt                                      <  1:2919136‑M1/49‑14 (MQ=255)
                ataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtatt                                   <  1:4059188‑M1/50‑17 (MQ=255)
                ataatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtat                                    <  1:1809914‑M1/49‑16 (MQ=255)
                 taatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattc                                  >  1:3961806‑M1/1‑33 (MQ=255)
                 taatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattc                                  >  1:3851653‑M1/1‑33 (MQ=255)
                 taatGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattc                                  <  1:153431‑M1/50‑18 (MQ=255)
                   atGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctc                                <  1:5532455‑M1/50‑20 (MQ=255)
                    tGTAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctct                               >  1:1955616‑M1/1‑30 (MQ=255)
                      tAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctgg                             >  1:32646‑M1/1‑28 (MQ=255)
                      tAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctgg                             <  1:1090259‑M1/50‑23 (MQ=255)
                      tAGCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctg                              <  1:4326116‑M1/49‑22 (MQ=255)
                        gCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggtt                           <  1:1611539‑M1/50‑25 (MQ=255)
                        gCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggt                            >  1:642351‑M1/1‑26 (MQ=255)
                        gCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggt                            <  1:621469‑M1/49‑24 (MQ=255)
                        gCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggt                            >  1:5337153‑M1/1‑26 (MQ=255)
                        gCCTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggt                            >  1:4079100‑M1/1‑26 (MQ=255)
                         ccTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggtta                          <  1:3711993‑M1/50‑26 (MQ=255)
                          cTCGTGTCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttat                         >  1:468887‑M1/1‑24 (MQ=255)
                             gtgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgt                       >  1:4098449‑M1/1‑21 (MQ=255)
                             gtgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgt                       <  1:1360333‑M1/49‑29 (MQ=255)
                             gtgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgt                       >  1:175874‑M1/1‑21 (MQ=255)
                             gtgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgt                       >  1:5243282‑M1/1‑21 (MQ=255)
                              tgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtca                     <  1:4830386‑M1/50‑31 (MQ=255)
                              tgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgt                       >  1:4590381‑M1/1‑20 (MQ=255)
                              tgtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgt                       >  1:1464732‑M1/1‑20 (MQ=255)
                               gtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcag                    >  1:2345239‑M1/1‑19 (MQ=255)
                               gtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcag                    >  1:5008130‑M1/1‑19 (MQ=255)
                               gtCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcag                    >  1:3586559‑M1/1‑19 (MQ=255)
                                tCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcag                    <  1:2184084‑M1/49‑32 (MQ=255)
                                tCTTGCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcag                    <  1:5166655‑M1/49‑32 (MQ=255)
                                    gCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggtttt               >  1:2511261‑M1/1‑14 (MQ=255)
                                    gCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttt                <  1:1914023‑M1/49‑36 (MQ=255)
                                    gCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttt                <  1:4296766‑M1/49‑36 (MQ=255)
                                    gCGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggtt                 <  1:4937017‑M1/48‑35 (MQ=255)
                                     cGAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggtttt               <  1:5586388‑M1/49‑37 (MQ=255)
                                      gAGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgc             >  1:2497026‑M1/1‑12 (MQ=255)
                                       aGGATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcc            >  1:2461008‑M1/1‑11 (MQ=255)
                                         gATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctg          <  1:532674‑M1/50‑42 (MQ=255)
                                         gATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctg          <  1:106283‑M1/50‑42 (MQ=255)
                                         gATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctg          <  1:4589082‑M1/50‑42 (MQ=255)
                                         gATAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctg          <  1:1509673‑M1/50‑42 (MQ=255)
                                          aTAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgc         <  1:2265455‑M1/50‑43 (MQ=255)
                                          aTAAGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgc         >  1:2114825‑M1/1‑8 (MQ=255)
                                             aGTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcga       <  1:943282‑M1/49‑45 (MQ=255)
                                              gTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcgaat     >  1:2621866‑M1/1‑4 (MQ=255)
                                              gTGCcgacttcaccagtattctctggttatgtcaggttttgcctgcgaat     <  1:1766999‑M1/50‑47 (MQ=255)
                                               tGCcgacttcaccagtattctctggttatgtcaggttttgcctgcgaatg    <  1:3051597‑M1/50‑48 (MQ=255)
                                               tGCcgacttcaccagtattctctggttatgtcaggttttgcctgcgaatg    <  1:1793301‑M1/50‑48 (MQ=255)
                                                 ccgacttcaccagtattctctggttatgtcaggttttgcctgcgaatgca  <  1:4671891‑M1/50‑50 (MQ=255)
                                                 ccgacttcaccagtattctctggttatgtcaggttttgcctgcgaatgca  <  1:3143919‑M1/50‑50 (MQ=255)
                                                 |                                                 
AAAAACTCTGAGTGTAATAATGTAGCCTCGTGTCTTGCGAGGATAAGTGCATTATGAATATCTTACATATATGTGTGACCTCAAAATGGTTCAATATTG  >  NC_000913/3888382‑3888480

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